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TRIP- TRAJECTORY INTEGRATION PROGRAM 


Bernard Kaufman 


ABSTRACT 


This document describes a numerical integration program (TRIP) based on 
Encke's method of integrating only the perturbations in the cartesian coordi- 
nates. The mathematical model can include Venus, Earth, Mars, Jupiter, the 

ol 


ness potentials are included for the Earth and Moon and Jj only for Mars. Mod- 
ification to expand the oblateness models for the other planets is simple. Solar 
radiation pressure may also be included as an option in the input. The program 
uses the latest exported JPL double precision planetary ephemerides (DE :^.9) 


TRIP has the option of calculating shadow times using penumbra and umbra; 
and printout may be obtained at specified times or, in the case of a planetary 
orbiter, at apoapsis or periapsis. Input data is very simple and there is a 
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TRIP- TRAJECTORY INTEGRATION PROGRAM 


INTRODUCTION 

TRIP was developed in 1968 primarily to provide an accurate and fast trajectory 
program to replace the n-body numerical integration programs then in use for 
mission analysis at Goddard. These were primarily large and very slow pro- 
grams such as the old JPL Space Program. A second and equally important 
objective was to make the input format simple enough to be used without extensive 
study of a manual. 

TRIP was developed along modular lines so that it could be easily modified. It 
uses a fourth-order predictor-corrector numerical integration subroutine 
(described later) which because of the modularity may be easily replaced by 
any integration package. Many options have been added to TRIP since it was 
first developed; however, since Encke*s method is used which requires rectifica- 
tion of the reference orbit, it was not considered advantageous to incorporate 
atmospheric drag. 

Throughout the development of TRIP, its speed and simplicity have been main- 
tained, making it a very useful as well as accurate, tool. It was most recently 
used for the launch window analysis for Imp-L 

TRIP is currently running on the IBM 360/95 and is available in the Goddard 
program library (number D00140). 


DESCRIPTION OF SUBROUTINES 

The functions of the subroutines used in TRIP are described briefly below. 
Many of them came from the Quick Look Mission Analysis Program developed 
by Ptsilco-Ford (reference 1). A complete listing of TRIP appears in Appendix 
A. 


MAIN 

ADOT 


The driver for TBIR Its primary function is to set up certain 
keys from the input^ta. 


(Reference 1) Obtains the an^e between two vectors. 


ABES - (Reference 2) Calculates the matrix necessary to transform co- 
ordinates from Earth mean equinox and equator of date to Mars 
mean equinox and equator of date. 
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ARKTNS 

CROSS 

DERIV 

DOT 
FNOL II 

FNORM 

GETACC 

GETTAP 

GHA 

GOTOR 

HARMON 

INTRl 

ITER 

KONSTK 

Ml 


(Reference 1) Obtains the arctangent of y/x. 

(Reference 1) Finds the cross product of two 3 -dimensional vectors 

Calculates the disturbing accelerations due to the presence of 
other bodies and* when required, calls other subroutines to calcu- 
late accelerations d^ o to oblateness and solar radiation pressure. 
Performs rectification test. 

(Reference 1) Finds the dot product of two vectors. 

(Reference 3) Used for the numerical integration of a system of 
ordinary differential equations. A fourth-order Adams-Moulton 
predictor- corrector with a Runge-Kutta starter, FNOL II has 
the option of automatically varying the step size. 

(Reference 1) Finds the magnitude of a 3 -dimensional vector. 

(Reference 4) Evaluates the disturbing acceleration due to 
oblateness of the Moon or Earth. 

(Reference 5) Used to position the ephemeris data set and read 
the correct data record, 

(Reference 1) Calculates the Greenwich hour angle. 

(Reference 1) Solves Kepler's equation for incremental eccentric 
anomaly on a conic section given the incremental mean anomaly 
(time) . 


Computes acceleration due to oblateness of Mars and is the driver 
for GETACC 

(Reference 5) Sets up the call to READE to find position and velocity 
of the other bodies with respect to the central body. 

Calculates times of entrances and exits from penumbra and umbra. 

Defines all constants needed by TRIP. 

(Reference 1) Multiplies a 3-dimensional vector by a 3 x 3 matrix. 
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MNA 

MULT 

NUTAIT 

ORB 

ORB2X 

OUT 

OUTX 

OVRLAY 

READE 

ROTEQ 

RVOUT 

SHADOW 

SHIFTP 

SOLAR 

SPEH 


(Reference 1) Pro' ies the matrix to transform moon-centered 
coordinates in Earth's true equator and equinox to Moon-centered 
coordinates m tlie moon’s true equator. 

(Reference 1) Forms the product of two 3X3 matrices. 

(Reference 1) Calculates the matrix to relate coordinates in 
Earth's mean equator and equinox to the true equator and equinox. 

(Reference 1) Transforms cartesian state vector to keplerian 
elements. 

(Reference 1) Transforms keplerian elements to cartesian state 
vector. 

Prints out various types of output as determined by keys set from 
input. 


(Reference 1) Prints out cartesian and spherical coordinates. 
Called from subroutine OUT. 

(Reference 1) Reads input data. 

(Reference 5) Heads s interpolates and translates ephemeris data 
from the JPL double-precision ephemeris system (DE-19) . 

(Reference 1) Provides the matrix to transform from mean 
equator and equinox of 1950.0 to mean equator and equinox of 
date. 


(Reference 1) Calculates spherical coordinates of position and 
velocity. 


Driver for subroutine ITER. Sets up keys to be used in calculating 
shadow times. 


(Reference 1) Provides the cartesian state vector of a probe 
relative to all the bodies in the ephemeris. 

Calculates the perturbing acceleration due to solar radiation 
pressure. 

(Reference 1) Converts cartesian position to spherical coordinates. 
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STEPD - (Reference 1) Obtains the state vector of a probe which is sepa- 
rated from an initial state by a time or angle increment, assuming 
only two -body motion. 

TCONIC - (Reference 1) Calculates time from periapsis corresponding to a 
given true anomaly on a conic section. 


TERM 


Contains termination logic for the numerical integrator. Deter- 
mines when the proper conditions for termination have occurred. 


TFRAC - (Reference 1) Provides the integral and fractional parts of the 
sum of two numbers. 

TIMEC - (Reference 1) Computes whole and fractional days from January Ij 
1950 at 0*"". 


TRAJ 


- The driver for FNOL II. 


INPUT 

All constants and variables needed for the program are contained in subroutine 
KONSTK and built-in values are supplied thf re for everything except the state 
vector and the epoch times. The state vectors option keys and epoch times are 
input through subroutine OVRLAY» which also provides for overlaying any of 
the built-in constants or variables. A complete description of the input necessary 
to run TRIP is given below. 

The input format is 4 (13 » D12.8). The first 3 digits are fixed-point input (13) 
and represent the C -array location of the variable in subroutine KONSTK. The 
C” array location is followed by the value desired. 


-array Name (if any) 

Variable 

750 

X (km) or a (km) 

751 

Y (km) or e 

752 

Z (km) or ^ (true anomaly-deg.) 

753 

X (km/ sec) or U (deg) 

754 

Y (km/sec) or 1 (deg) 

755 

Z (km/sec) or cd (deg) 

756 

Epoch time (year, month, day in form of 
yymm.dd) 

757 

Epoch time (hour, minutes, seconds in form 
of hhmm.sec) 
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Reference coordinate system of input 

1. Earth mean equator and equinox of 1950, 0^, 

Jan, 1, 

2. True ecliptic and equinox of date, 

3. Selenographic-true moon equator — rotating , 

4. Mars me^ equator of date, 

5. True Earth equator and equinox of date. 

Stop time (year, month, day as above) 

Stop time (hours, minutes, seconds as above) 

TINC - total days from start to stop time if 
only 1st and last point, or printout at peri- 
apsis or apoapsis is wanted, or printout at 
approximate steps of At is wanted (nearest 
time step to At). 

TINC = A t if printout at exact intervals of At is 
wanted. 

Built-in value is 180. (days) 

1, If printout is wanted at the 1st and last points 
only, or if printout at periapsis or apoapsis 
is wanted, or if printout is wanted at ap- 
proximate steps of At, 

total days If printout is wanted at exact steps 
M of At. 

Built-in value is 1. 

For printout at exact steps, an additional print 
occurs if TOTAL * At ^ stop time. 

IQO 000. If printout at 1st aixi last points is 
wanted, or if printout is wanted at periap- 
sis, apoapsis or exact intervals of At. 


0. If printout is wanted at approximate 

steps of At 


Built-in value is 0. 


0. E print it is wanted at first and last points 
only; or if printout is wanted at periapsis , apo 
apsis or exact steps of At. 

At (seconds) E printout is wanted at approxi- 
mate steps of At. 














AH£A 


Cross sectional area of spacecraft in cm^. 
Built-in value is 22225. Used only if c(769) = 1. 


771 


772 


XMASS 


Eeflectivity 
of spacecraft 


Total mass of spacecraft in gms. Built-in 
value is 113000. Used only if ci 


^Values from 0. to l.j inclusive. Used only if 
c(769) = 1. 

1, Perfect reflector 
0. Perfect absorber 
Built-in value is 0.2 


773 


SPEC. REFLECTION 
of spacecraft 


Fraction of particles specularly reflected. 

1 . 


0 . 0 % 

Built-in value is 1. May be 0, to 1., inclusive. 
Used only if c{769) = 1. 


Solar pressure at 1 A.U. 

774 CONSTANT Mlt-in value is 4.7 x 10”® dynes/cm 

Used only if c(769) = 1. 


775 INPUT TYPE 


0. Cartesian. 

1. Orbital elements. 

Built-in value is 0. Actual values are input 
throu^ c(750) — c(766). 


776 APPER 


777 OUTPUT = 


Apoapsis and periapsis key. 

0. Do not print at apoapsis or periapsis. 

1. Print only at apoapsis. 

-1. Print only at periapsis. 

Built-in value is 0. 

Output options in form of abcdefgh. 

1 a; Central body meem Earth equator of 1950 
b: Stenographic 
c; Blars mean equator of date 
d; Central body true ecliptic and equinox of date 
et Earth-centered true ecliptic and equinox of 
date 
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f: Sun-centered true eclipt’ ? and equinox of 
date 

gi Central body true Earth equator and equinox 
of date 

ht Earth-fixed latitude and longitude (only if 
^ Earth is central body) 

0 Do not print this form. 

1 Print this form. 

Built-in value is 1000(J011. 

Cases may foe stacked. When this is desired, cards for the 1st case are fol- 
lowed by a blank card and then the cards for the second case. There is no limit 
to the number of cases that may be stacked. Again, since an overlay method is 
used on input, only those variables different from those for the preceding case 
‘ need be input. A blank card at the end of the last case is necessary to end the 
program. A sample of the ii^ut is shown in Appendix B. 


OUTPUT 

The output options are listed above in the description of array c(777). TRIP has 
the capability of automatically switching body centers when the probe leaves the 
sphere of influence of one body and enters that of another body. This switching 
does not interfere with the printout options. 

Printout for the central body consists of both the cartesian state vector and 
instantaneous osculating elements; for other bodies, only the state vector. A 
sample of the output with an e^lanation is shown in Appendix B. 


CONSTANTS 

Ail constants needed by TRIP are listed in subroutine KONSTK and are either 
described there by a comment card or are described above in the input section. 
The constants in c(31) throu^ c(38) are no longer used by the program but were 
not taken out. 

Ai^ or all of these constants may be changed for a particular run by including 
the proper array number and the new value in the input data. 

KNOWN RESTRICTIONS AND ANOMAUES 

Ihere are two restrictions to be noted in using TRIPs 
1.) Only one type of printout can be used In a run. 


2. The printout interval must be greater than two minutes which is the initial 
time step in the integrator* If a request for a smaOer interval is made, the 
program automatically sets the printout interval to be slightly greater than 
2 minutes and prints out a message to this effect. 

There are also two known anomalies in the program, which cause inconvenience 
but do not produce incorrect results: 

1. When periapsis printout is requested, occasionally the program will inter- 
sperse a printout at apoapsis and then continue with periapsis printouts. 
Similarly, a printout at periapsis may occur when printout is requested at 
apoapsis * 

2. When shadow times are requested, if the satellite is already in shadow at the 
time the shadow search is initiated, the printout of the shadow time is incor - 
rect for that shadow interval only . All subsequent printout are correct. This 
bug occurs because the search procedure used requires a point before shadow 
occurs in order to iterate to the correct time. These incorrect shadow times 
will be obvious to the user (see ease 6, Appendix B). 

A final note of caution — when exact printout is asked for, the numerical inte- 
grator must restart after each printout time in order to continue. This procedure 
may result in a loss of accuracy due to numerical problems? that is if 2 cases 
are stacked using the same trajectory, one asking for exact printout and the 
other asking for printout only at the last point, the final state vectors at the stop 
time may disagree somewhat. To a certain extent, the numerical errors may be 
minimized by increasing the size of the error col rol in array c(767> or by de- 
creasing the rectification limit in array c(720). H jwever this is a problem in- 
herent to a certain extent in any numerical integration routine and probably can 
not be completely eliminated. 
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C ENCKE PROGRAM MAIM , « o o o » » 100 

IMPLICIT REAL* 8 { A-H, 0 “Z? S) . 200 

DIMEMSION K( 3 I »V (3 ) yH 1( I » VI ) 300 

DIMENSION OUM(3 93 »9nATEFC<3v3JynATECT( 393)yXX(3l9)(V(3) AOO 

COMMON Cl lOOOIylCI^OtZO) ,VV(3) 9M(JR9INP SOO 

DIMENSION XN0T(3I AOO 

DIMENSION ELKA) 700 

COMMON/APPER/ FPRnLD 90 L 09 KMAX flOO 

COMMQN/JTERM/MEARTHffMMOnNjMMARS 900 

DATA BLK/ 6 H / » XNOT 1 /3HNUT/ 1000 

EQUIVALENCE (C (401 ) 9 R 11 1» 1 ) ) 9 t C( 425 1 9 V I f 1 p 1 M 9 1100 

1 (C1750I,xn ) ) 9 (C( 753 ) ,V( 1 U 1200 

DIMENSION A(393)9EN(393I»EMN(393),DATE(3v3)9SDATF(393)9 1300 

1DATET( 39 3» 9SDATET(393) 9 AOATEI 3p 3 1 9 ADATET ( 3p 3 ) 9 AR ( 3 , 3 ) 14 00 

COMMON/KOT/ date, SDATE 9 DATET 9 SOATFT 9 AOATF 9 ADATFT 9 OATFEC pDATECT 9 FN 15 00 

CnMMON/XLIH/ SG 9 TA 9 XMLIE IfrOO 

COMMON/TOP/TLU 1700 

COMMON /SHAD/SHADK IPOO 

COMMON/ PR I NT/ IC 5 O 9 I SElpMARDaTp ICFNECp I FAR EC, ISUNFC p ICOATE 9 lEARFX 19 00 

CALL KONSTK 2000 

1 CALL QVRLAY 2100 

C(462»=lo00 2200 

C(122)=loD0 2300 

FPROLUsOpOO 2400 

OLD = O.nO 2500 

IF(C(775)oEOoO»00) GO TO 2 ?f»00 

NTFMP=C(766I 2700 

DU 5 1=1 96 ?R00 

5 FLHII = C(I+749) 2900 

CALL 0RB2X(X9V9ELS9C(NTEMPH 3000 

2 CONTINUE 3100 

IF( FNORM(X I .EOsOoOOl GO TO 20 3200 

SHAOK«OoODO 3300 

TS=C(756J 3400 

0F=C(757) 3500 

INPsC(75ai 3600 

ST^C(759) 3700 

STF=C(760) 3800 

TINC=C(761) 3900 

TOTAL»Ct762) 4000 

C( 1000)=C(765I 4100 

XNhR=C(766> 4200 

NEsC(767) 4300 

MEARTHaC(76BI/100»DO 4AO0 

TEMP=Ct76BJ-MEARTW*100 4500 

MMOONsTEMP/lOoDO 4600 

TEHP»TEMP-MMDON*10 4700 

MMARSaTEMP 4800 

IC50=C(777»/10000000900 4900 

TEMP“C(777I-IC50*10000000 5000 

ISEL*TEMP/1000000«00 5100 

TF«P=TEMP-ISEL*1000000 5200 

MARDAT«TEMP/ 100000*00 5300 

TEMP«TEMP»MAR0AT*100000 5400 

ICENEC*TEMP/10000,D0 5500 

TEMP«TEMP-ICENEC*10000 5e>00 

gEAREC*TEMP/1000*DO 5700 

TEMP«TEMP-!EAREC*1000 5fiOO 

ISUNEC»TEMP/100«00 5900 

7EMP-TEMP-ISUNECPI00 6000 

!COATE»T£MP/I0*f)0 6100 

TEMP»TEMP-IC0ATE*10 6200 

IEARFX«TEMP 6300 

NTOTAL*TOTAL 6400 

NOR*XNOR 6500 

CALL TEMEUTSfOFpTWpTFl 6600 

CALL TIMEC 8 ST 9 STP 9 ST 1 »STF 1 » 6700 

A-1 



IFCCdOOOn 60950,50 

6R00 

60 

NR1TE(6930) 

6900 


GO TO 70 

7000 

50 

WRITE(6940J 

7100 

40 

FORMAUlHl 9 l9HSHAnOH IS REOUESTFO) 

7200 

30 

■FORMAT (1H1 923HSHA00W IS WOT hEOIIFSTEnP 

7300 

70 

COWTINUE 

74 00 


XWOm )=BLK 

7500 


XN0T(2)=BLK 

7600 


XN0T(3)=BLK 

7700 


IF(MEAKTH,LT,n XWOTm = KWUTi 

7R00 


IFIMMOOWoLTd » KNOT! 2)*KN0T1 

7900 


IFIMMARS.LTol ) XNOT I 3 ) sXNOT 1 

8000 


WRITEt69l00f Cf61l jKNOTf IP 

8100 


WRITE( 69 lOOJ C(62»9XfM0H2) 

8200 


WRITE( 69 lOO) C(65) 9XN0T(3P 

8300 

100 

FORMATUHO, 13HHARM0NICS OF 9 A 694 H IS 9 A4,RH IWCLODFDJ 

8400 


XWOTf 1 )*BLK 

8500 


IF(C4 769J oEOcOqODO) XW0T( D^XWOTl 

8600 


wRiTE( 69 io 2 P xworm 

8700 

102 

FORMAT! oQ SOLAR RAOIATIOW PRESSllRF IS S A4, “ INCLUOED * ) 

8800 


TIME=TW+TF 

8900 


TLUi=TlME 

9000 


CALL R0TE0ITIME9A) 

9100 


CALL WUTAIT!TIME9nM9CR9 0T9EN9ePSm 

9200 


CALL MWA{T1ME90M,CR90T9EPSIL9K0909GP9WW9FMM) 

9300 


XML IB=180.1)0*,01745329252DO * CR-OM-^TA-SG 

9400 


CALL MULTIEN9A90ATE) 

9500 


OUM! I 9 U-I .000 

9600 


DUM!l92I*=0.00 

9700 


OtJMi l,3P=0»00 

9800 


0UM!29Ua0.00 

9900 


DlJM(2 92 J=OCnS(EPSlL P 

10000 


D0M(293P«0SIW!EPSILP 

10100 


OOM (3911*0. 00 

10200 


DUM(3p2P*-DS IN (BPS ID 

10300 


DUM(393P*OCnS(EPSIL P 

10400 


CALL MULTtmJMpOATEgOATEFCP 

10500 


CALL MULT!EMN9DATF,S0ATe' 

10600 


call ARESaSpDFpARP 

10700 


CALL MULT! ARgAp AOATEP 

10800 


00 iOi 1*193 

10900 


00 101 J*U3 

11000 


OATECTiJg 1 P*DATEEC! I, JP 

lUOO 


DATETUpn-OATEnp JP 

11200 


SOATET!Jpl P-SDATEU pj) 

11300 

101 

AOATET !J 9 1 P * ADATE ! 1 9 J P 

11400 


WRlTE!6p6P 

11500 

6 

FORMAT! IHOpBlHlNPOT IN MEAN EO ANO EO OF 1950P 

11600 


CALL OATOUTITW 9 TF 90 AYM*FI)ATE 90 P 

11700 


DO 600 1*193 

11800 


XX! n«xn p 

11900 

600 

XV! I P*V! n 

12000 


GO TO!200«30094009SOOp700P«INP 

12100 

300 

CALL MItXfDATECTtXXP 

12200 


CALL MHVpDATECTtXVP 

12300 


GO TO 200 

12400 

400 

CALL MUXpSOATETtXXI 

12500 


CALL MUVpSDATETpXVB 

12600 


GO TO 200 

12700 

500 

CALL Ki!X»ADATETpXX» 

12800 


CALL MUVpiM^ATET«XVP 

12900 


GO TO 200 

13000 

700 

CALL NHX90ATET«XXP 

13100 


CALL MHVpOATETfXVI 

• 13200 

200 

CONTINUE 

13300 


00 11 I=lf3 13A00 

C^!■^499)=^XSr I) 13500 

11 C( |+502J=KV(I » I3ft00 

CALL 0UTK(K)f*KV»NnUT) 1370U 

C(592)-TW l3ftno 

C(593)=TF 139O0 

TFF=TF 1*000 

KKsl 14100 

13 TFFaTFF+TINC 14200 

CALL TFRAC(TW9TFF«TWFffTFF( 14300 

16 KMAX=0,DO 14400 

CALL TRAJ(TWF»TFF#XXbXV9N6 J 14500 

IFrC«462).EO*O.O0)GO TO 1 14600 

TW=TWF 14700 

DO 10 JJ=i»3 14R00 

C( JJ+499)*XX ( JJ) 14900 

10 CIJJ+502)=XV< JJ» iSOOO 

ST0T=(Sn-C(592f + STFl-C(593l >«86400o00 15100 

IF( ST0T-120,D0,LE»0oD0J GO TO 1 15200 

XMAX-OoOO 15300 

!F( caOOOJ oGT,0,D0) Cl 10O0»=0oD0 15400 

14 1F1C1462J .eOoOoOOf GO TO 1 15500 

IFlFPROLOoEO. 0*001 GO TO 12 15600 

FPR0LO»0*O0 15700 

OLD =0*00 15BOO 

STOTPRs* 1 T WF-Cl 592 l+TFF-C 15931 l*R6400o00 15900 

IF! ST0TPft*LEol20*D0) GO TO 12 16000 

IF IKK .LE» MTOTALl GO TO 16 16100 

GO TO 15 16200 

12 KKaKK+i 16300 

15 IF(KK*LE,NTOTAL 1 GO TO 13 16400 

CALL TRAJ(ST1,STF1»XX9XV,WE) 16500 

IF(FPROLD*NE*OoOO) GO TO 14 16600 

GO TD 1 16700 

20 CONTINUE 16R00 

STOP 16900 

end 17000 


o d 


R»=AL FUNCTIOW AD0T«8(XpV> 17100 
IMPLICIT R6AL*flf A-HpO-Z,S) 17R00 
DIMEWSION X(3)p Y(3) 17300 
FUNCTION RETURNS ANGLE IN OEGREFS bETWEFW TWO INPUT VECTORS 17400 
RANGE IS RESTRICTED TO (0»180) 17500 
A=FNORN(X) 17ft 00 
BsFNIJRMIVI 17700 
ANG=OOT( XpVf/A/P 17R00 
AnOT = ARKTNS( 18 Op ANGpOSORTI 1, 000- ANG*ANG I )«5 7* 295 779D0 17900 
RETURN IftOOO 
END 18100 
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SUBROUTIPJR ARES (T 1 »T2fl AR ) 

IMPLICIT RFAL«8( A-H»0-Z 9* ) 

DIMENSION AR(3v3)flEE(3l9nO),FXO(3) -FOni 

dimension TEMI3I 

DUMaOoOO 

DTR=0,017<V5329252D0 

IYR=n/100oD0 

YR=1YR 

DT=YR-50oD0 

AnP= 3 1 7c 793A16667D0+0 oA 5 208330-2* DT 
OOP=5^o6575DO+Oo35D-2*DT 
BYR-YR*l00oD0+lc01D0 
CALL TIMEC (BYRjOIJMjTW, TF ) 

CALL TIMECITIpTZ^TWItTFII 

TAU= < (TWl-TWJ+tTFl-TF) 1/365.2500 

A0= AOP-0.OOlO13DO*TAU 

D0= DOP-0. 00063 ID0*TAU 

TE= ( TW1+TF1+IR262.5DOI/36525.DO 

nM= .7 86441 6700+0. 77 099 16 700* TF-Oc 1388Bft89D«5*TF.*TE 

X 1= 1. 85033333300-0. 6750-3*TE+0. 1261 11 lin-4*TE*TF 
E= 23.4522944400-0, 1301250- 1*TE“0. l63R88fl9D-5*TF*TE 
1 +0.5027777«0-6*TE*TE*TE 

AOR=An*DTR 

OORa 00*OTR ' 

OMRs OM*DTR 
XIR= XI*DTR 
ER=E*DTK 
SADsDSINIAOR ) 

CAO=DCOSI AOR ) 

SOO=DSIN(DGR I 
COO=DCOS(DOR) 

SOMsDSINIOMR > 

CnMs=DC(JS(QMR) 

SI=OSIN(X IR) 

CI=DC0S(XIR) 

SE=OSIN(ER) 

CE*OCQS(ER» 

C the following is KAUFMAN® S METHOD +♦««•* 

eei n*cAo*CDo 
EE(2)sSA0*C00 
EE(3I=SD0 
01 1)=SI*S0M 

012)= -SI*COM*CE-CI*SE 
0( 3)«-Si*C0M*SE+Cl*CE 
call CROSSlOfEE^EXQ) 

EOm=-SAO 

E0I2)*CA0 

EQ(3)-0.D0 

CALL CROSSIEO^EXQfTEM) 

SOMEG*DSORT(TEM( n**2+TEM(2I**2+TEMI3»**2) 

COMEG* OOT(BO»EXO) 

0MG»ARKTNS1180i^C0MEGvS0MEG) 

SOMEG*DS!NCOMG) 

COMEG»OCOS(OMG) 

AOM - AOR+1 .57079632700 
XAI*la570796327DO-OOR 
SAOM»OS1N< AOM) 

CAOM«OCOS(AOM) 

SAI«DS!N(XAn 
CA1-0C0S(XAI ) 

ARI I V U-COMEG*C AOM-SOMEG«SAOH*CAI 
A R U 9 2 ) *C ONE G*S AOM+ 50MEG*C AOM* CAI 
AM le3)>S0MEG*SAI 

AR I 2« U— SOMEG*C A0M-C0MEG*SA0M*CA1 
AR4 292)«°50MEG*SAOM+COMEG*CAUM*CAI 
ARS2»3)-COMEG*SA! 

ARC 391)bSAOH*SA| 
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18200 

18300 

18400 

18500 

18600 

18700 

18800 

18900 

19000 

19100 

19200 

19300 

19400 

19500 

19600 

19700 

19800 

19900 

20000 

20100 

20200 

20300 

20400 

20500 

20600 

20700 

20800 

20900 

21000 

21100 

21200 

21300 

21400 

21500 

21600 

21700 

21800 

2J900 

22000 

22100 

22200 

22300 

22400 

22500 

22600 

22700 

22800 

22900 

23000 

23100 

23200 

23300 

23400 

23500 

23600 

23700 

23800 

23900 

24000 

24100 

24200 

24300 

24400 

24500 

24600 

24700 

24800 


ARO,2)a-CAOM*SAI 2^900 

AR13p3)«CAI 25000 

RETURN 25100' 

END 25200 
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REAL FUNCTION ARKTNS^S (N, K 9 Y ) 25300 

IMPLICIT REAL*a(A-H9D*Zp4) 25400 

COMPUTES 4-UUADRANT ARCTANGENT OF Y/K IN RADIANS 25500 

Nb360 angle lies IN RANGE (O 936 O) DFG 25600 

N=180 ANGLE LIES IN RANGE (-lflO»lBOI DEG 25700 

TPI * 6»2831853000D0 25AOO 

»A - DABS(») 25900 

YA = DABSIV) 26000 

IFUA-YA) UI 92 26100 

1 Z»X/VA 26200 

GO TO 3 26300 

2 Z=V/XA 26400 

YA*)CA 26500 

3 D=0S0RTUo0D0+Z*Z ) 26600 

VAsYAOD-^X 26700 

IFIYA»4,4,5 26800 

4 ARKTNS«=TP!/2*0D0 26900 

GO TO 6 27000 

5 ARKTNS*2*ODO*OATAN<y/YA) 27100 

6 IF(N-180)7,997 27200 

7 IFIARKTNS )8»999 27300 

0 ARKTNS=ARKTNS+TPI 27400 

9 RETURN 27500 

END 27600 


A-7 


SUBROUTlwe CROSS(A,B*C) 27700 
IMPLICIT REAL*8(A-H,0-Z*$} 27flOO 
DIWeMSlOM AI 3) ,B(3 I »C(3) 27900 
DIMENSlOfy A(3)9RI3)*C( 3) 28000 
Cm-A(2)*B(3)-A(3)>fB<2) 28100 
CI2)sA(3»#Bm-Am'»'B{ 3) 28200 
Cl )«B( U 28300 
RETURN 28A00 
END 28500 
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SUBROUTINE DATOUT ( T I » T2» OATE p FDATE» KHATE ) 2fi600 

IMPLICIT REAL#8< A-HpO-ZpS)- 28700 

DATE RETURNED IN FORMATIVEARS FROM 1900 M MONTH ),( DAY S > 2R800 

FOATE RETURNED IN FORMAT ( HOUR SM M I N) » ( SFCONDS) 2R900 

K SET 0 FOR PRINTOUT? SET - OR + FOR NO OUTPUT 29000 

GIVES CALENDER DATE FROM TKWHOLF DAYS FROM 19S0) 29100 

ANDT2(FRACT OF DAY) 29200 

DIMENSION C(12) 29T00 

C DIMENSION C<12) 29400 

DATA 0000HL/6H JAN*/ 29500 

C( 1 )=0000HL 29600 

DATA 0001HL/6H FEB,/ 29700 

C(2 J^OOOIHL 29R00 

DATA 0002HL/6H MARCH/ 29900 

C(3 )”0002HL 30000 

DATA 0003HL/6H APRIL/ 30100 

Cl 4 )=0003HL 30200 

DATA 0004HL/6H MAY/ 30300 

C(5 )=0004HL 30400 

DATA 0005HL/6H JUNE/ 30500 

U6 ) = 0005HL 30600 

DATA 0006HL/6H JULY/ 30700 

Cl 7 )=0006HL 30800 

DATA 0007HL/6H AUG*/ 30900 

C18)»0007HL 31000 

DATA 0008HL/6H SEPT./ 31100 

Cl 9)s0008HL 31200 

DATA 0009HL/6H OCT,/ 31300 

Cfl0)«0009HL 31400 

DATA 0010HL/6H NOV,/ 31500 

C111)=Q010HL 31600 

DATA 0011HL/6H DEC,/ 31700 

C112)=0011HL 31R00 

MD»0 31900 

K-Tl+365c0D0 32000 

IFlTl)6fl6p7 32100 

6 WRITE 1 6pl00) 32200 

8 RETURN 32300 

100 FDRMAT123M DATE IS 1950 OR BEFORE) 32400 

7 K1*K/I46l 32500 

K2*K“K141461 32600 

KWN=K2/365 32700 

KKsK1*4+1949+KWN 32800 

N*1 32900 

IFlKHN-3)U2plllp 111 33000 

111 N«0 33100 

112 KWN«K2-KWN*365 33200 

IFIKWN) II f 11 »13 33300 

11 MONTH* 1 33400 

M0«1 33500 

GO TO 25 33600 

13 CONTINUE 33700 

KNW*0 33800 

JJ*0 33900 

IBS 34000 

IFa2'°JJ)23923pl37 34100 

C JeF vM9A«MpJ,J9A»S*O«N«M0NTHS FOR GO TO 34200 

137 GO TO lL49l59l49l6«l4»16»l49i4»16»l4pI69l4)pJJ 34300 

14 CONTINUE 3V400 

KNW«KNW43i 34500 

GO TO 17 34600 

15 CONTINUE 34700 

IFlNU0«189l9 34800 

18 KNW-KNU+l 34900 

19 KNW*KNW428 35000 

GO TO 17 35100 

16 KNW*KNW-«'30 35200 
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17 CQNTINUc 35300 

It^iKWW-KNW 120*21,22 36400 

22 CONTINUE 35500 

«D«KNW 35600 

GO TO 135 35700 

20 CONTINUE 35800 

MDNTH^*JJ 35900 

MD®KWN-MD+l 36000 

GO TO 25 36100 

21 CONTINUE 36200 

MONTHsJJ-U 36300 

MD«1 36400 

GO TO 25 36500 

23 MONTH =12 36600 

MD=KWN-MD*H 36700 

25 continue 36800 

" H=T2424oODO 36900 

THP«IDINT(TH1 37000 

FDATE=THP*100.0DO 37100 

NHOUR=THP 37200 

THP=ITH-THP )460o000 37300 

TH=IDINTITHPJ 37400 

FDATB*PDATE+TH 37500 

NMIN=*TH 37600 

THP=(THP-TH»460o0r)0 37700 

NSEC=THP 37800 

F0ATG*F0ATE+THP/100,0D0 37900 

TSEC*NSEC 38000 

0ATE=IKK-1900)’*'100+M0NTH 3R100 

SP=MO 38200 

DATE=DATE+SP/10Oo0D0 38300 

IF(KOATE» 8,50,8 38400 

50 CONTINUE 38500 

THP=(THP-TSEC 1*1000*000 38600 

NFSEC»THP 38700 

TT»!0INT{T2+,5O0) S^HOO 

TP2=T2+*5D0-TT 38900 

TPl*TL'>'TT-s-2433282,0D0 39000 

WRITE I 6,l01)C<M0NTH),M0,KK,NHnUR,NMlN,NSEC ,NPSEC,TP1 ,TP2 39100 
101 FORMAT! IH , A6, 1 3, IH, , IK , 14 , IH , IK, I 2,4HHRS, , I3,4HMIN, , 13,1H* , 13, 39200 

13HSBC,50X,12HJULIAN DATE , F8. 0,T1 1 1 ,F9*8 1 39300 

GO TO 8 39400 

END 39500 
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SUBROUTINE DERIVIXpYsD) 39600 

IMPLICIT REAL*8( A-HsO-Z9St 39700 

DIMENSION V(20»9R1(3»«)9VK3»8I ,0(20) vCB( 20l9NP«5) » 39A00 

lDACC(3J9Zl(3l9VVl(3)pRE(3),CACC{3) sRZOpB) 39900 

DIMENSION DACCI3I 40000 

DIMENSION SACC(3» 40100 

COMMON Ca000),lC(50»9Z(3)pW(3i9NOR 40200 

COMMON/JTERM/MEARTHpMMOONpMMARS 40300 

COMMON/DISTRB/DACCpCACCpOACC 40400 

COMMON /SHAD/SHADK 40500 

EQUIVALENCE ( C ( 401 ) pR U I * I H » ( C ( 425 ) * VI ( 1 9 L ) P 40600 

C MUST CALCULATE POSITION OF PROBE FROM INCREMENTS 40700 

OT=X 40800 

U=:C(NOR) 40900 

IFIXoNEpOoDOI go TO 321 41000 

DO 322 i=lp3 41100 

ZI n*C( i-s-499) 41200 

322 VVn )=Cn+502) 41300 

GO TO 323 41400 

321 CONTINUE 41500 

MN=l 41600 

CALL STEPD(MNpDTpTApC<500) pC(503»pUpZp VV p I pCB) 41700 

323 CONTINUE 41800 

DO 10 I = U3 41900 

znn=zm+vm 42000 

10 vvim=vv(n+vn+3i 42100 

TW=Ct592) 42200 

TF=C(593I+tC(719)+C«506I+X) /86400 p00 42300 

CALL TFRAC(TWpTFpTWpTF) 42400 

CALL SHIFTP(NORpZlpVVl,TWpTFpRlpVU 42500 

NP4 1)=2 42600 

NPt2J»3 42700 

NPI3>»4 42800 

NP(4)«5 42900 

NP|5J«6 43000 

GO TOt Ip2p3p90p80p 100) pNOR 43100 

2 NP( U«l 43200 

GO TO 1 43300 

3 NP(2)*1 43400 

GO TO i 43500 

90 NP43)«1 43600 

GO TO I 43700 

80 NPI4)«1 43800 

GO TO 1 43900 

100 NP(5I»1 44000 

GO TO 1 44100 

I CONTINUE 44200 

DO 4 l-lpS 44300 

NPl«NPm 44400 

RN{ n«*FNQRHiRmpNPm«*3 44500 

00 5 J«l»3 44600 

5 R21 J^I^'BMRl CJpNORI-RlIJpNPU 44700 

4 RMim«FN0RM{R3(|pl43H**3 44B00 

NPl«NPm 44900 

NP2*>NPS2) 45000 

NP3«NPO) 45100 

NP4»NPS4I 45200 

NP5«NPI5S 45300 

DO 6 l»l«3 45400 

6 0ACCm->aNPU4iRH!*NPUMMI U4R2Ut4»/RNU in-CiNP2)4IRllUNP2S 45500 

IMNI2KR2U«5I/RNU2H 45600 

2-UNP3l«Uin«NP3»/R*H^4(12n»6mMH3n-€«NP4l«4RiH«NP4l 45700 

3/RNI4I 4R2nvTI /RNU4BBH:«NP5£4«faU«NP5»/RN4 5l«R2HvflI/RMUSBI 45B00 

DO 7 I«i»3 45900 

7 REm«zn»4.5D04Vn S 46000 

0« D0TIRE»Y}/SFN0RHUI4429 46100 

SUM -I.DO 46200 
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TERMaKOO ^6300 

DO 8 K=2v50 ^6A00 

XK=K 46500 

TERMaTERM4((.5D0-XK)*2oD0>«‘0/(XK»l#D0) \ 46600 

SUMaSUM+TERM 46700 

IF( DABS( TERM) oLToOoOOOOlOO^SUM) GO TO 9 46800 

8 CONTINUE 46900 

WRITE(6»20) 47000 

20 FORMATnHl,34HO SERIES IN ENCKE 010 NOT CONVERGE) 47100 

STOP 47200 

9 FO=lcDO-SUM 47300 

H?CCN0R)/( FNORM(Z )**3) 47400 

00 11 1=193 47500 

11 CACC( i )=-H4Vn )+H*FO*Rin sNOkI 47600 

00 30 1=193 47700 

30 0ACCn)*0*D0 47800 

GO TO (4095092l92l960921)9NOR 47900 

40 IFIMEARTH.NE*!) GO TO 21 48000 

’ GO TO 70 48100 

50 IF(MMOONoNEol) GO TO 21 48200 

GO TO 70 4B300 

60 IFIHMARS.NEa ) 60 TO 21 48400 

70 CONTINUE 48500 

TWTE=C(592) 48600 

TFTExCI593)+(C(506)+X)/8640O.D0 48700 

CALL TFRAC(TWTE9TFTE9TWTE9TFTE) 40800 

TIMUP=TWTE+TFTE 48900 

CALL HARMONIRivVltOACCiiTlMUP) 49000 

21 CONTINUE 49100 

DO 13 1=1 j3 49200 

13 DAccm* OACcm+oAccm 49300 

RAT = FNORM«Y(U) / FNORMtZd)) 49400 

IF(RAT#GEoC(720)*0R.09GEoC« 720)) CU22) = 0o0DC 49500 

00 ir 1=1,3 49600 

17 0Accm«DAccn )+cAccm 49700 

IF(C<769)9E0.0.0D0«0R.SHA0KeE09 1«000) GO TO 551 49800 

4 CALCULATE SOLAR PRESSURE ♦*♦#*# 49900 

CALL S0LAR{Rm»3),SACC) 50000 

00 1000 I«l,3 50100 

1000 DACCm-OACCU ) * SACCm 50200 

551 C O N T I N U E 50300 

DO 16 J-1,3 50400 

Dm-YIHS) 50500 

16 Oi |4-3)»DACCn ) 50600 

RETURN 50700 

END 50800 
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REAL FUNCTION D0T*8(K«Y) 

50900 

IMPLICIT REAL*8( A-HjO-Z,*) 

51000 

DIMENSION 

51100 

DIMENSION )^(3I»V(3I 

51200 

OOT ' xm»i«Ym X(2)*Vt2I + X(3)*Y(3) 

51300 

RETURN 

51400 

END 

51500 
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SUBROUTINE FN0L2(J9N9G9L9H«NE 97,0,009 17,00) 

51600 


IMPLICIT REAL*8( A“H,0-Zfl*) 

51700 


COMMON BERNIE(IOOO) 

51800 

51900 

52000 


DOUBLE PRECISION )(0,VD,VA, VC, YP,Yl 

52100 


DIMENSION Y(50)9D(50),YB(3096),GI2430)9GI3(30) pGl4(30) ,EF(30), 

52200 


DIMENSION VtaOOO) 

52300 


DIMENSION ZZ(5),0(50},VB43096),GI2(30),G!3(30> ,G 14(30) ,EF(30), 

52400 

lEFl ( 30) «EF2 ( 30},EF3 f 30) 9 VI { 30) vERRDRt 30» «HA( 301 9 YA« 50), DA (50) « 

52500 

2VC( 30) 9YPI30) 9YD(50) 

52600 


EC*Y(N+3) 

52700 

9876 

FORMAT (IHO, 17 ) 

52800 

1 

HbG 

52900 

2 

HZ*H 

53000 

3 

LN-N-»-MAX0{L9 3) 

53100 


SUBROUTINE FN0L24J,N9G,L,M,NE ,X »V,D*0ERIV9TERM,0UT) 

53200 

6876 

F0RMATUH0,17) 

53300 

A 

NA«0 

53400 

5 

NB«1 

53500 

6 

NF*0 

53600 

7 

NG-0 

53700 

8 

F«0«00 

53800 

9 

FA>0.D0 

53900 

)0 

FB^OoOO 

54000 

11 

FC>*0«00 

54100 

12 

CONTINUE 

54200 

13 

ENE»NE 

54300 


DO 200 I-lfiLN 

54400 

; 200 

VD4 n»0BLE(V4I )) 

54500 


X0«D6LE4X) 

54600 


DO 200 !»i,LN 

54700 

200 

Y04 n-vn ) 

54800 


xo>x 

54900 

14 

IFU-3H5,2l9l5 

55000 

IS 

IF4N£)ie«l6,18 

55100 

16 

JA=4 

55200 

17 

GO TO 22 

55300 

18 

RE1«10.00**(-ENE) 

55400 

19 

RE2*10,O04»«-ENE“3*0D0S 

55500 

20 

REM«i0.004>*4-ENE>1.500) 

55600 

21 

JA*1 

55700 

22 

00 25 I«i,N 

55800 

23 

00 24 KC«1,5 

55900 

24 

YB4 1,10-0.00 

56000 

25 

ERROR! I I-O.OO 

56100 

26 

CAU DERIV4X,V,D) 

56200 


CALL TER«HX«VvO,F» 

56300 


|F4DABS4F)-I.OO"9) 731,731, 5209 

56400 

§209 

CONTINUE 

56500 


DO 300 1-1, N 

56600 


GI2U )-om 

56700 


Giam-om 

56800 


Gi4m-om 

56900 


EFin-OlU 

57000 

: 27 

CALL OUTPUT! X, V,0, ERROR, N,L,H) 

57100 

27 

CALL OUT«X,VtO,ERROR0N,L,H) 

57200 


IF!S£RN1£8§06».NE*0.00» GO TO 29 

57300 

2« 

F0-VIN4U 

97400 

29 

lFM-21 30,129,30 

57500 

30 

GO T0«3l, 37,39,371, JA 

97600 

31 

00 33 

57700 

32 

vAn»-vo«n 

57800 

33 

oA«i»-o«n 

57900 

34 

GO TO 37 

58000 

3§ 

Mft-N 

98100 

36 

H-2*00«H 

98200 
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37 HD2 * •5nO*H 5R300 

DO 39 I-IbW 58400 

38 VBU9NB)BD(n 58500 

XL « 0(1) * HD2 58600 

C 39 V( n»SNGUVOn f ’^XL t 58700 

39 Y(I)»VO(n+XL 58000 

C X«SNGL(XD+HD2) 58900 

X**X0+H02 59000 

40 CALL DERIV (X9Y9GI2) 59100 

41 DO 42 I»L«N 59200 

XL * GI2( I )*HD2 59300 

C 42 Y( n>SNGLl YD( n-t-XL) 59400 

42 Y(n*VD(n+XL 59500 

43 CALL DERIV (X^Y^GIS) 59600 

44 DO 45 IsljN 59700 

XL-GlBUJfH 59800 

C 45 YU )mSNGL(YD( n^XL ) 59900 

45 Ym«YDm■^XL 60000 

C X*SWGL(XD+H) 60100 

XeXD*H 60200 

46 CALL 0ERIV(X«Y9GI4) 60300 

47 HD6 “H/6*D0 60400 

GO T0( 48» 55 96 O 966 t « JA 60500 

48 DO 52 8*l9N 60600 

XL*(D(n ■»- 2.00*(GI2(I) G13(m -»-GI4U n»H06 60700 

49 VC«n-YDtn+XL 60800 

51 YOU )*tYA(l) 60900 

52 ERRORn)«0«D0 61000 

53 JA«3 61100 

54 GO TO 35 61200 

55 DO 57 I-UN 61300 

XL»(Dm * 2«D0*(G12m * GI3Un +GI4m>4tHD6 61400 

56 VD(| IwYOn )*#*XL 61500 

C 57 ERRORS! )*SKGLfY0m-Y9>n n/15. 61600 

57 ERRORm»( YD(E )-YP(I H/|5«D0 61700 

58 JA^I 61800 

59 GO TO 681 61900 

60 DO 62 62000 

61 YOm»YCm 62100 

XL«SOm + 2.004(GI2m GIBU)) ^GI4m)«MD6 62200 

62 YPU KYAtl I-i-XL 62300 

63 62400 

64 JA-2 62500 

65 GO TO 661 62600 

66 DO 66 I«lvN 62700 

XL«SOm *■ 2«D0*(GI2m * GI3U H •*'014(1 H«HD6 62800 

67 Y0U)«Y0U)4>XL 62900 

66 ERRORU )»0<»DO 63000 

66i 00 69 I«lfN 63100 

C 69 Ym-SNGUYOUn 63200 

69 Yn£«YOm 63300 

XD-XIHH 63400 

C XaSNGUXDl 63500 

X«XO 63600 

70 CALL DERlVUfV^D) 63700 

71 FC»F 63600 

72 call TERMIXffVfOfFl 63900 

73 IFIOA6SaFS<-’l*OOo'9 ITSlvTMvTSB 64000 

731 NF-9 64100 

732 GO TO 124 64200 

733 IF«F 874, 124*76 64300 

74 FA-1.00 64400 

75 GO TO 77 64600 

76 FO-1.00 64600 

77 1F4FA*FBB83,78,63 64700 

78 WF-WF-H 64800 
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79 

JA=4 


64900 

80 

NB=l 


65000 

81 

HsH*F/<FC-Ft 


65100 

82 

IF(NF-4)37,37,124 


65200 

83 

IF(NE S84v 117,84 


65300 

84 

IF( JA-mi7,B5, 117 


65400 

85 

IF< J»3>06,U7,86 


65500 

86 

DO 95 1*1, N 


65600 


IFtVn H886, 885,886 


65700 

885 

HA( n*iOOO.DO 


65800 


GO TO 95 


65900 

B86 

!F4EC»B80,890,87 


66000 

07 

IFfOABSUn n-EC) 880,880,890 


66100 

880 

iFf DABS! ERROR n H-RE2} 882,94, B81 


66200 

881 

!F(DA8S(ERR0R< I H-REU94,94,882 


66300 

882 

HAn)*Hl«tREM/JOABSlERRORn)) + ,0000000001DO) >*♦1,200) 


66400 

883 

GO TO 95 


66500 

890 

1F( DABS ( ERROR m/V (I ) ) •RE2)892, 94,89 1 


66600 

391 

IF(DABS«ERROR| I) /V 41 )) -RE 1 >94,94,892 


66700 

892 

HA( n-H*(REM/C0ABS(ERR0R4 U/Y4 1 >>+* 000000000100) >** ( , 

200) 

66800 

893 

GO TO 95 


66900 

94 

HA4 U»H 


67000 

95 

CONTINUE 


67100 

96 

HBbOA8S4HA(N) ) 


67200 

97 

00 98 1*1, N 


67300 

C 90 

HB*AM1N1 tABStHAU )>,HB) 


67400 

98 

HB*0H1NUDABS4KA4 1 >),HB) 


67500 

99 

IF!DABSlH)-HB)100,n7, 101 


67600 

100 

|F(OABS4H2)-DABS4H>) 101,101,116 


67700 

101 

00 103 1*1, LN 


67800 

102 

VD4n»VA< S > 


67900 

C 

Y( n*5NGL( V04 n } 


68000 


V4 n-von } 


68100 

103 

mU^QMU 


68200 

104 

IF4NB-6) 107,105,105 


68300 

105 

XO*XO-»M 


68400 

106 

GO TO 109 


68500 

107 

X0*XD-2.00*H 


68600 

108 

H£*H 


60700 

109 

H*OSlGNf HQ,H> 


68800 


X*XD 


68900 

C 

X*SNGUXO) 


69000 


CALL DER3VtX,V,0) 


69100 


CALL TERMU,Y,0,F) 


69200 

110 

NB»1 


69300 

111 

XAB$*OABSf «000001D0*X > 


69400 

m 

IF40A8$<M)-XABS} 113, 113,117 


69500 

113 

NG»NG^1 


69600 

U4 

H*0SIGNUABS«H> 


69700 

115 

IFiNG - 10)124,150,150 


69800 

150 

WRIT£(6,126U H 


69900 

1261 

FORMAT UH1,107HEXECUT1QN TERMINATED BECAUSE INTERVAL 

OF INTEGRATi 

70000 

ION LESS THAN l.OE -6 TIMES INOEPENDENT VARIABLE UK 

X »,1PD15,7> 

70100 


STOP 


70200 

116 

H£»H 


70300 

U7 

iFiMma,iie,i2i 


70400 

U8 

IFI0ABS!VIN4l)-F0)-YiN4>2> >29,119,119 


70500 

119 

F0*F04YIN*2> 


70600 

120 

GO TO 124 


70700 

121 

NA«NA<«>1 


70800 

122 

IF4M-NA>123tl23,29 


70900 

123 

NA»0 


71000 

C126 

CALL OUTPUT U,Y,0,£RR0R,N,LfH> 


71100 

124 

CALL OUTU,Y,0,ERROR,N,L,H» 


71200 

125 

IFINF'>4)29,29,126 


71300 

126 

CONTINUE 


71400 





o n 


128 RETURN 71500 

129 NB=NB+1 71600 

130 IF(NB»6130fl I3lfl 136 71700 

13lD0l34Issl9N 718 00 

132 BF3n 1 = YB<U31 71900 

133 EF2U 1 = YB1 I ) 72000 

13A EFin l*YBn p51 72100 

135 GO TO 137 72200 

136 NB=10 72300 

137 HD24 =H/ 24,00 72400 

DO 138 I»UN 72500 

KL =(55oD0*0( n -59,00*EFHn +37,00*FF2m -9 , D0»FF3 ( II 1 *H024 72600 

YPl n=VOl I )+KL 72700 

138 Y( n=SNGL( YPII n 72800 

K = SNGL fXn-9-H) 72900 

138 vm*YPm 73000 

X=XO+H 73100 

139 CALL OfiRiVIXpYDEFl 73200 

140 DO 142 IslsLN 73300 

141 YAI ! ) = YDn 1 73400 

142 DAUlsOm 73500 

143 DO 148 1=1 jN 73600 

XL = (9,00*EFm +19oD0^0m "5onO*EFim +EF2mi*H024 73700 

1'4 YD( I )=YD( IKXL . 73R00 

C 145 FRRORm = -SNGL(YOm-YPn n/14, 73900 

145 ERRORm=“(YDl I )-YPtI M/14,D0 74000 

145 EF3U fsEF2U 1 74100 

147 EF2(I )=EF1 t I ) 74200 

148 EFK I t=l)t I ) 74300 

999 CONTINUE 74400 

149 GO TO 681 74500 

END 74600 



REAL FUNCTION FN0RN«8(X) 

74700 


IMPLICIT REAL*8<A-M^O“2jtI 

74R00 


DIMENSION Xt3I 

74900 


DIMENSION X(3) 

7*5 000 

1 

FNORM*DSORT(K« 1 J f 21 *>C2+X I 3 J 2) 

75100 

3 

RETURN 

75200 


END 

75300 


SiJBROUTINE GETACC ( LCT s 1 CT? CNM* SNM* ECN, PVA? CTT 9 CTB 9 UMT »VMT jNTE ) 75400 

IMfLICIT REAL* 8 ( A-H^O-ZjU 75500 

DIMENSION LCTlDfflCn U*CNM( U, SNMt UsECWm . 75AOO 

DIMENSION PVAUKCTTmflCTbt nfOMTI 1) 9 VMTt U 75700 

EOUIVALENCE I XP 9 PX > 9 ( VP? PY) 9 1 ZP »PZ ) 75ROO 

ER0=ECN(1I 75900 

XMU*eCNI2I 7ft000 

ALFbECNO) 74100 

PXs PVA( n 76200 

PV= PVAI2) 76300 

PZ* PVAI3i 76400 

CGBsOCOSi ALF) 76500 

SGB°0SIN( ALFI 76600 

RMTaKOOO/(PX*PX+PV*PV+PZ*PZ ) 76700 

TMAaRMT*ERD 76B00 

CLT*TMA*I XP«l'CGB+YP*SGBI 76900 

SLT*»TMA*( YP*CGB-XP*SGB) 77000 

RZR«TMA*ZP 77100 

RMRaTWA*ERO 77200 

TWR=0.5D0/ER0 77300 

RMOaDSORTIRMTI 77400 

LMX»0 77500 

LNC=l 77600 

LA*tCTl2t 77700 

IF (NTE-III 59 I 592 O 77R00 

15 LA»0 77900 

20 LBaLCim 78000 

LC*LA 78100 

NB»NTE+1 78200 

DO 75 I»l9NB 78300 

KA»LA-l 78400 

KB-LB 78500 

KC-LC+1 78600 

IF (KB-KA)30»35«35 78700 

30 KB«KA 78800 

35 IF ILCI 5 O 95 O 94 O 78900 

40 IF (KB“KCI45j50»50 79000 

45 KB*KC 79100 

50 LNC»LNC+KB+1 79200 

ICTm*lNC 79300 

LA*LB 79400 

L8*LC 79500 

(.C-LCm+21 79600 

IF aMX-LAJ60*65»65 79700 

60 LMX-LA 79800 

65 IF n*NTE + U75?70s70 79900 

70 LC«0 80000 

75 CONTINUE 80100 

THA>2,0D0 B0200 

TMB-Z.ODTj 80300 

DO 90 l>lfLNX 80400 

CTTm— TMA 805 00 

CTBm> TMB 80600 

TNA«TMA-«-2«OnO 30700 

TH&«THB'4>TNA 80800 

90 CONTINUE 80900 

UMT! U«XNU*RNO 01000 

vNTm»o«ooo siloo 

LA-l 01200 

t&«l 81300 

NA*2 01400 

K6«2 81500 

00 160 1-1 »NB 01600 

KC-ICim-1 01700 

KA-!«1 01800 

TNB-KA+KA 01900 


A-19 


■ ! 

i 


i 



’i 

'"I 


TMAsTMB+UODO 

TMC«1.0DO 

IF (‘^C-KBUAOfl 100» 100 
100 DO 130 LC«KB»KC 

TM0»R2R^UMT(LB)^TMA 
TMEaRZR*VMT(LB)*TMA 
IF (KB-LCI 110» 120^120 
110 TMD=TMD''RMR*UMTaR"U*THfl 
TM6 = TMe-RMR*VMT (LB-U*TMB 
120 UMT(NA)*TMD/TMC 
VHT(NAI=TME/TMC 
T?1AeTMA+2Q0D0 
TMb'* TMB+ 1 • 000 
TMC«TMC<-loODO 
NA»NA-«-l 
130 LB*LB+1 

140 IF n~NB)1509 I60e 160 
150 TMC«KA+KA+1 

UMHNA>sTMC*(CLT*UMT( LAJ-SLT*V«T( LA) ) 

vmt(na)stwc*<clt^vmt(la)>slt*umt« lam 

NAbNA-^1 

KB»KC+2 

LB-LB4-1 

LAsLB 

160 CONTINUE 
TMA*0«0D0 
TMB»OoODO 
TMC=OoODO 
TMDa0o5D0 
KB^l 

LA* ICT(I) 

LB*2 

LCaLA 

00 260 I*UNTE 
KC>»LCTn ) 

IF (KC)260»2609210 
210 IF n-'2)230»2159220 
215 LA=3 

GO TO 225 
220 LA»ICTU-2)+2 
225 LB«ICTU-n+l 
LC*ICTm 

230 DO', 250 NA«1«KC 
TME*CTB«NAJ 
TWF*TME 

IF { 1-1)235*235,240 
235 TMFtaTMO^TWe 
THE»NA 
TM0*7M0*TM6 
TME«NA+2 
TMD»TMD/TME 
TWe«“TMF 

240 RMR-T!ME4UMT(LA) 

R2R«TMF«VMTILA) 

CLT-UMTJLC) 

SlT«VMmC) 

TME-RMR-CLT 

TNF«R2R+SLT 

RMR*RMR-*>CLT 

R2R»R2R-SLT 

CLT*CNMU6) 

SLTaSNMCRB) 

TMA> TMA<»>CL T T4>R 2 R 

TM6«THtt-CLT*T«F+SlT*RHR 
TMC*TMC^CTTINA)i*JCLT*U«T4iBS'»>SLT*VMT*LeJ S 

LA«tA4-l 


82000 
82100 
02200 
82300 
82400 
82500 
82600 
82700 
82R00 
8 2*? 00 
83000 
83100 
83200 
83300 
83400 
83500 
83600 
83700 
83800 
83900 
84000 
84100 
84 200 
84300 
84400 
84500 
84600 
84700 
84800 
84900 
85000 
85100 
85200 
85300 
85400 
85500 
85600 
85700 
85800 
85900 
86000 
86100 
86200 
86300 
86400 
86500 
86600 
86700 
86800 
86900 
07000 
87100 
87200 
87300 
07400 
87500 
87600 
87700 
07800 
87900 
88000 
88100 
86200 
88300 
88400 
88500 




LB«LB+1 

8H600 

LC-LC+1 

80700 

KB=KB+1 

08800 

COWTINUE 

88900 

tma«tma*twr 

89000 

TMB=TMB#TWR 

89100 

PVA(7t=(CGB«TMA-SGB*TMB) 

89 200 

PVA(8 ) = <SGB*TMA+Cr,B*TMB) 

89300 

PVA(9J=TMC*TWR 

89400 

RETURN 

89500 

END 

89600 


A-Jl 



SUBRnUTINE GETTAP 89700 

DEFINE FILE 104 95ft p6636v L» I V ( 89ft00 

DOUBLE PRECISION TABS? JDI p JDF* JDI, TDAY, JOIF p JDP 89900 

REAL NUT AT 90000 

REAL^e TA831 90100 

COMMON/CETBL2/ICW,ICENTpIRE04 13 P 90200 

C0MM0N/CETBL3/TAB34829J,NUTAT420^P pCKSUW 90300 

C0MM0N/CETBL9/JDlpTDAVp JDIFpIERRl 90400 

DATA lV/0/ 90500 

TATA IWPftV/0/ 90600 

JDPaJDl+TOAV 90700 

IF( J0P-LT,2438760,500.0RoJDPeGTo2446424«5n0)G0 TO 100 90ft00 

IERR1»0 90900 

IrnCWaNEal) GO TO 205 91000 

DO 333 I«l9204 91100 

333 NUTATIDbO. 91200 

ICM=2 91300 

IW»1 91400 

READ! 10° IWITAB3 91500 

TAB31=TAB3m 91600 

205 IW=ID|NTnj0P-TA63U/8<.D0J+lo00 91700 

IFI IW.EOoI WPRV tGO TO 220 91B00 

IWPRV»IW 91900 

200 READ( 10® IWJTAB3 92000 

IF4 JDP-TAB3tU)2109220p230 92100 

210 IW=IW-2 92200 

GO TO 200 92300 

230 IW*IW+1 92400 

I F { ODP- 4 TAB3 4 1 ) +8 . DO ) ) 220, 200p 200 92500 

220 JDIFsJDP-TAB3m 92600 

RETURN 92700 

101 F0RMAT4® EPHEMERIS ERROR — JULIAN DATF °pD25.7p'IS BEFORE JAW 1965 92A00 

lOR AFTER DEC 1985® J 92900 

100 MR;T 8(69 iOUJDP 93000 

STOP 93100 

END 93200 


A-n 


SUBRDUTIWE aHA(TSFC9D*GHAMtDA*UM6GA) 93300 

IMPLICIT REAL*B< A-HjO-Z I 93^00 

0MEGA*o00Al7807^200/< 210-lB^D) 93500 

DD*0*DO 93600 

00*0 93700 

DD»DD*( »9856^735D0/360*n01 93P00 

0«=1D!NT<DDI 93900 

D0**00-0« 94000 

DF=DD 94100 

TEM 1«* 100 *07554D0-»'360oD0*nF+2o90 150-13*040 94200 

140MEGA*TSEC 94300 

1 IFITEMl)2»3s3 94400 

2 T6MlaTEMl+360oO0 . 94500 

GO TO 1 94600 

3 IF( TEMl-360oD0»5,4,4 94700 

4 TEMUTEMl-360,00 94R00 

GO TO 3 94900 

5 GHANaTEMH-OA457o2957795nO 95000 

RETURN 95100 

END 95200 


A-Xl 



o o n n o 


C 


C 


SUBROUTIwe GOTORiK,VMoC»F»En 
IMPLICJT REAL*8U-H^0“Z9S> 
DIMENSIOM C(2)pF(4» 

DIMENSION Cr2)pE(4> 


NMAK«20 

N=0 

GO TO (192), K 

1 CONTINUE 

GUESS IS OBTAINED FOR ELLIPTICAL CASE 
ft CONTINUE 

IF( El-lo0D0J30e31p31 

30 CONTINUE 
02»E1*E1 

F( 1 IsEl^Oa*! o l6666667D+00-D2«( 0 833333333D-02-D2* ( ol9R4 126980-03-02 
l*( 9275573192D-05-02«. 2505210830-07 )M ) 

F( 2)«D2*l *5D0-02'<'Io416666670-01-02'M o 138ft8«fl9D-02-D2#( <,248015870-0 
14-D2»*27557319D“06 > )) ) 

F( 3>*E1-F( 1) 

F(4)=1bOOO-F(2) 

GO TO 3 

31 CONTINUE 

F( 3)sOSIN( EU 
F(4)*0CQS(E1 ) 

F( 1)=E1-F(3) 

F(2)=1o000”F(4) 

GO TO 3 

2 CONTINUE 

FIRS? GUESS IS OBTAINED FORHYPE RROL ! C CASE 
9 CONTINUE 

IF( E1tK000)32933933 

32 CONTINUE 
D2*E1«61 

F( 1) ■ei»D24( o l66666670+00-t-D2* (.83333333D-02+02* to 19841 aaOeo-OS-frOa 
1*( *275573l92D-05-«-D2*.250521083D-07n » ) 

F(2)«02*t0.5D0+02>M .4l6666670-01+D2«t*138888890-02+D2^( *248015870- 
104+02*. 275573190-06)) ) ) 

F(3)-El+Fm 
FI4)»1.0DO+Ft 2) 

GO TO 3 

33 CONTINUE 
£X«.5D0*0EXPtEl) 

0X». 2500/EX 

F(3)«EX-0X 

Ft4»*EX+0X 

Ft i)«Ft3)“61 
F(2)*Ft4)-l.000 . 

3 CONTINUE 

CM-FtU+Cti )*F(3)+Ct2)*Ft2B 
DM«Ft2)+CU)*Ft4)+Ct2I*FC3» 

0E*tVM-C«) /DM 

ERRORbOE 

AB«DABStDE) 

1F(AB-1.000UO«10« U 
U DE-0E/A6 
10 E2-EI+0E 

IF HE2+EU.E0* 0.00) GO TO 4 
IFt0ABSUE2-EU/(E2+Ein-3.0-8)494«5 
5 CONTINUE 

IFtN-NMAXB6v7«7 
7 CONTINUE 
GO TO 4 


95300 
95400 
95500 
95600 
95700 
95800 
95900 
96000 
96100 
96200 
96300 
96400 
965 00 
96600 
96700 
96B00 
96900 
97000 
97100 
97200 
97300 
97400 
97500 
97600 
97700 
97ft00 
97900 
9 ft 000 
98100 
98200 
9R300 
9ft400 
9«500 
98600 
98700 
98800 
98900 
99000 
99100 
99200 
99300 
99400 
99500 
99600 
99700 
99800 
99900 
lOCOOO 
100100 
100200 
100300 
100400 
100500 
100600 
100700 
100800 
100900 
101000 
lonoo 
101200 
101300 
101400 
101500 
101600 
101700 


A-M 



ooezoT 

00220 ? 

00I20X 

00020? 

006101 

OOHTOl 


0M3 

30NUNO0 *f 
Oi. 09 
23«ig 
T ■f>N»N 

gnwiiNoo 9 



SUBROUTINE HARMON < R 1a V lp ACC» T I MF ) 

IMPLICIT REAL#fl <A-HfO-Z,AI 
DIMENSION DATEEC(393)ffDATECT( 3*3) 

DIMENSION RU3*B)pVH39«)*RR(3) * ACC( 3)* TFMP( 3) *OATEt3*3) * 

1 SDATE(3*3) 9DATET(3*3)9SDATET(3* 3 )* ADATE(3* 31* ADATET(3,3) 
DIMENSION Al3*3)*EN(3*3)*EMN(3f 3)*AR(3*3)9LCT(fl)* ICTI9) * 

1 CNM( 36)*SNMt3A) *ECNI3) 9PVA(9)*CTT(R J ,CTRI R) *UMT( 45) *VMTi45 ) * 

2 LCTMl 5)* ICTMI 6 ) *CNMM< 15)*$NMMI 15) » CTTM| 5 ) * CTRMI 5)*«)MTM(2l ) * 

3 VMTMI21) 

COMMON/ROT/ DATEeSDATE90ATET*SDATFT* ADATF* ADATET*DATEeC*DATECT*EN 

COMMON/XL IB/ SG*TA*XMHB 

COMMON CaOOO)* tC(50) *Z( 3) *VV(3)*NQR 

COMM ON/ TUP/ TLU 

DATA LCT/0 9 7* 6 9 5 *4 9 3 * 2 * I / * LCTM/5*4f 3 * 2* 1/ 

EQUIVALENCE I C ( 200 ) *CNMm ) * ( C ( 236 ) * SNM ( 1 )) 9 ( U 272 ) *CWMM ( U ) 

EQUIVALENCE I C I 287 ) * SNMM U H 

KEV*1 

1F4 DABS<TLU-TIME) 0 LT 9 C(901) ) GO TO lOU 

TLU=TIME 

KEY=2 

CALL R0TEQ4TIME* A) 

CALL NUTAmTlME,0M*CR9nT*EN*EPSlU 
CALL MULT4EN*A*0ATE) 

DO 101 1=1*3 
DU 101 J = l*3 

101 DATETtJ,!) = 0ATE(1*J) 

100 CONTINUE 

GO TO 120*4*80*80*5)* NOR 
20 CONTINUE 
NTE =8 

.ECNIi) = C(21) 

ECN12) = cm 
ITIME = TIME 

T 2 » (TIME-ITIME )*86400o0n0 
T1 = ITIME 
DA « EN12*l) 

CALL GHA(T2*Tl*GHAN»0A»nMEGA) 

ECN13) * GHAN*f.01745329252D0 
CALL HI(Rm«n«DATE*RR) 

DO 1 1=1*3 
1 PVAl n = RRII ) 

CALL GETACCtLCT*ICT*CNM,SNM*ECN*PVA*CTT 9 CTB 9 UMTpVMT*NTE) 

CALL Ml 4PVA17)*DATET»ACC) 

RETURN 
4 CONTINUE 

GO TO 1200*202)* KEY 

202 CALL MNAmME 90 M*CR* 0 T*EPSIL«R 0 *G* 6 P*WW«EMN) 

CALL MULTI EMN*0 ATE* SOATE) 

DO 201 1=1*3 
DO 201 J=l*3 

201 SDATETUfl) = SOATEIUJ) 

200 CONTINUE 
DP-TIME 

TP=0P/36525.D0 

T2P»TP*TP 

T3P=T2P*TP 

0M-12.1I279D0-.O52953922D040P 0020795DO+TP 

1 ♦•0020eiO0*T2P+,2D-5*T3P 

CR=64*375452D0 ♦IS* ITGB^TDO^OP -.OOllSlSTHDO^TP 
X “■.00H3015OO*T2P+ • 190-5»T3P 

EPSIL=23*4457587D0».01309404D0»TP -•880-6*T2P ♦*50*6#T3P 
0N»ON*C175) 

CR=CR«C175T 

6PSIL-EPS1L4C875) 

XNLlfi » leo.oo*. 0174532925200 ♦ CR'OM^TA^SG 
NTE«5 


102500 
102600 
102700 
102800 
102900 
103000 
103100 
103200 
103300 
103400 
103500 
103600 
103700 
103800 
103900 
104000 
104100 
104200 
104300 
104400 
104500 
104600 
104700 
104800 
104900 
105000 
105100 
105200 
105300 
105400 
105500 
105600 
105700 
105800 
105900 
106000 
106100 
106200 
106300 
106400 
? 06500 
106600 
106700 
106800 
106900 
K'TOOO 
107100 
107200 
1 07300 
107400 
107500 
107600 
10T700 
107800 
107900 
108000 
108 100 
108200 
10R300 
108400 
1 08500 
108600 
108700 
108800 
108900 
109000 


A-li 



ECN(l) = C(22l 

lOQlOO 


ECW<2» = a 2) 

109200 


ECN(3> « OoOO 

109300 


CALL MURK lB2)«SnATEpRR) 

109400 


no 2 I“U3 

109500 

2 

pvA( n = RR( n 

109600 


CALL GETACC ( LCTMo ICTMpCNMMj SMMM»ECN« PVApCTTM,CTBMp|IMTM»V«TMflWTE) 

109700 


CALL Mil PVA(7)pSnATETpACC) 

109R00 


return 

109900 

5 

CONTINUE 

110000 


KJ2 = C(725I 

110100 


KJ3 = Cl 726) 

110200 


RE = C(25) 

1 10300 


U = C<5) 

110400 


R « FNORMjRl ( lp5 ) ) 

110500 


GO TO 4300»302)p KEV 

110600 

302 

ITIME = TIME 

1 10700 


T2 » TIME-ITIME 

110800 


Tl = ITIME 

110900 


CALL DATOUT « T 1 p T2 ,0 ATE , FOATE , 0 ) 

111000 


CALL ARESlOATEpFOATEf AR) 

1 UlOO 


CALL MULTJ AftpApADATE) 

U1200 


00 301 1-1,3 

1 11300 


DO 301 J*l,3 

111400 

301 

AOATETtJp 1 ) ■ ADATE(IpJ) 

1 11500 

300 

CONTINUE 

111600 


CALL MURl(l,5)p AOATEpRR) 

111700 


R2 s R4*2 

UlflOO 


R5 * ft**5 

1 11900 


R7 » R5*R2 

112000 


RE2 » RE*RE 

112100 


RES = RE2*RE 

112200 


22 B RKI3)’«<RRt3) 

112300 


23 B 22*RRI3) 

112400 


COEF ■ -1 .SDO^U^XJaifREZ/RS 

112500 


DO 10 Is Ip 3 

112600 

10 

PVAn+6) = COEF*t-5,DO'*22*R«( n/R2 + RRIIM 

112700 


PVAt9J B PVA(9) + C0EF*2,D0*RRI 3) 

112800 


!FUJ3.EO,0«D0) GO TO 13 

112900 


COEF * -,5DO*U*KJ3^RE3/R7 

113000 


DO 14 I«lp3 

113100 

14 

TEMPI!) » COEF*(-35oOO*2 3*RRm /R2 + l5,D0*RRm*RRI 3) ) 

113200 


TEMPO) ■ TEMPO) ♦ COEF*j 15, D0*22-3,00'»R2) 

1 133O0 


00 15 1*10 

113400 

15 

PVA(I^6) B PVAU+A) + TEMPI!) 

113500 

13 

CALL MlIPVAI7)pA0ATETpACC) 

113600 


CONTINUE 

113700 

60 

RETURN 

U3800 


END 

113900 


A-S7 



SUBROUTINE INTRlITWvTFtNH^PUS? I niJM, VFL » DDMiJ ) 

IMPLICIT RBAL *8 (A-H,0-Z) 

C THIS VERSION OF THE 360 JPL OnUBlE PRECISION EPKFMFRIS 

C WAS MODIFIED BY JERRY S* LINNEKIW OF COOP 551 ftSFC »GREFNBELT MO, 

C IT HAS BEEN MODIFIED TO HANDLE SFVEN BflOlES INSTEAD OF THE 

C NORMAL ELEVEN, THE BODIES AND ORDER ARE EARTH* MOON , SUN , VENUS 9 MARS » 
C SATURN AND JUPITER, SATURN AND JUPITFR HAVE BEEN INTFRCHANGEO 

C TO CORRESPOND TO THE QUICK LOOK MISSION ANALYSIS PROGRAM, 

C THE READ ROUTINE-GETTAP HAS BEEN MODIFIEO TO USE THE DIRECT ACCESS 

C READ FEATURE OF 360 FORTRANH,IN ADDITION, THE NUTATION VARIABLES 

C HAVE BEEN ELIMINATED FROM THIS VERSION SO THAT AN EPHEMERIS 

C RECORD WILL FIT ON ONE TRACK OF A 2316 DISK PACK, 

C THE DECK AS IT APPEARS HERE IS SET UP AS A DIRECT REPLACEMENT 

C FOR THE SINGLE PRECISION INTRI IN THE QUICK LOOK MISSION 

C ANALYSIS PROGRAM, 

C IN ADDITION TO THE ABOVE CHANGBSpTHE FOLLOWING COMMENT CARDS 

C DESCRIBE THE CHANGES NECESSARY TO USE THIS PACKAGE WITH THE 

C MARK2 PROGRAM ,THE MARK2 COVERAGE PROGRAM AND THE QUICK LOOK 

C COVERAGE PROGRAM 

C **j^**j!t MARK2 PROGRAM 

C SUBROUTINE ANTR I TW,TF ,NB,POS, VFL I MK2 PROGRAM 

C DATA NCENT/3, 11, 10,2,4,5,6/ MARK2 

C NCENTft=NCENTINB I FOR MARK II 

C *«>*♦*♦ MARK2 COVERAGE PROGRAM^*’***^*** 

C SUBROUTINE ANTR 1 (TH,TF ,NB,POSvKK, VEL,OI SI MK2 COVERApE PROGRAM 

C DATA NCENT/3,11,10,2,4,5,6/ MARK2 

C NCENTR=^NCENT(NB ) 

C COVER “QUICK LOOK COVERAGE PROG 

C NO CHANGE NECESSARY DUE TO CARFFUL PLANNING BY PROGRAMMER 

C TO EXTEND THE fc°HEMERIS FOR ALL BODIES, CHANGE THE FOLLOWING 

C P0S(2U TO P0St33i ,VEL(21) TO VELI33) AND NCENT(7» TO NCENTIll) 

C REPLACE THE NCENT DATA CARD WITH 

C .DATA NCENT/3, 11, 10,2,4,5,6,1,7,8,9/ 

C THIS WILL RESULT IN THE FOLLOWING BODIES IN THE ORDER 
C THAT THEY ARE LISTED EARTH,MOQN, SUN, VFNUS, MARS, JUP ITER, SATURN , 

C MERCURY, URANUS, NEPTUNE AND PLUTO, 

Q « « « « « « * ^ « 4i « « * * 4 « * 4t * 4r * 4t 4< 4< 4< >l( « * 41 >(< * « * * )(t 4t * « * 41 4c 4t 4( * )(< « i$i 41 m * « « 4i 41 )9r « 

C ARGUMENT DESCRIPTION 

C TW INPUT, WHOLE DAYS FROM 1950, REAL *8 

C TF INPUT, FRACTIONAL DAYS FROM 1950, RFAL*R 

C NB INPUT, CENTRAL BODY NUMBER U-EARTH, 2-MnON, 3-SUN, 

C 4-VENUS, 5-MARS, 6 -SATURN, 7-JUPITER), 1NTEGFR44 

C POS OUTPUT, PO'^ITION ARRAY, POSt U -POS ( 3 ► , EARTH X,Y,Z 

C VEL OUTPUT, velocity ARRAY, RFALfS 

C NCENTft central BODY NUMBER ( 3-EARTH, U-MOON , 10-SUN, 2-VENUS, 

C AU ASTRONOMICAL UNIT 

C REM EQUATORIAL RADIUS OF EARTH 

C 4-MARS, 6 -SATURN, 5-JUPlTER) 

C EMRAT EARTH MOON MASS RATIO 

REALMS JEO,NUT 
REAL *8 JEOPRV 
REAL «8 JON 
REAL *8 JEZ 
REAL48 JDi 

DIMENSION P0SI2n,VEU2U,NCENT{7l , IREQf 131 * TABOUT I 6 , 121 , NUT (41 

COMMON /CETBLl/ AU,REM,TPO, EMRAT 

COMMON /CET6L2/ ICW,NC£NTR, IREO 

COMMON /CETBL4/ TABOUT,NOT 

DATA NCENT/3, U, 10,2,4,6,5/ 

DATA JE0PRV/0«00/,IERR/0/ 

DATA IFIRST/0/ 

DATA NCENTP/0/ 
rF(IF|KST«N€»0»G0 TO 20 
icw>i 


114000 
114100 
114200 
1 14300 
114400 
1 14500 
114600 
114700 
Il4fl00 
114900 
U5000 
115100 
115200 
1 15300 
115400 
115500 
115600 
115700 
115800 
115900 
U6000 
116100 
116200 
1 16300 
116400 
116500 
116600 
116700 
U6B00 
116900 
117000 
117100 
117200 
1 17300 
U7400 
117500 
117600 
117700 
117800 
117900 
U8000 
118100 
118200 
H8300 
118400 
118500 
118600 
118700 
118800 
U8900 
U9000 
119100 
119200 
119300 
119400 
119500 
119600 
119700 
119800 
1 19900 
120000 
120100 
120200 
1 20300 
120400 
120500 
120600 


A-48 



IF iRSTsl 
AU*149597900oD0 
REM-6378 p 149200 
TPD«=86400,DO 
EMRAT=81o30lD0 
20 DO 10 I = U13 
10 IRE0(n»2 

IF« TF-<,5OOl30,39j39 
3fl SECl=TF«fl6400*D0 
GO TO 40 

39 DT=TW+*5D0 

SEC1 = ( TF-e5D0)»J‘R6400oDn 
GO TO 41 

40 DT*TW 

41 JD1=OT+2433282«5DO 

4 ^ 4c 4^ « « 4 :«i 4c > 01 4i 4 4i * * 4t ^ lit 4 4i 4 4 V 4c 4> <9t 4c Xt 4c ^ ^ 4i ^ i4c !4< 4( « 4 4 4 i4t i4t 

TSEC=SEC1 
JED=J01 
JEZ=JED+TSEC 
NCENTR=:NCENT (N6+1) 

IF( JEZdEO® JFDPRVoflNDol'dCENTRcEOoNCENTP) GO TO 1 

JEDPRVocJEZ 

NCENTPsMCENTR 

CALL READEt JEO 9 TSEC 9 lERR) 

1F( lERR^EOoOJ GO TO 1 
WRITE( 69 l 00 » IERR 9 NCENTR 9 ICW 9 JD3 
100 FORMAT(23HOFPHEMERI S ERRORj I ER R = 9 I 3 9 21 5* ni4o2 ) 

1 CONTINUE 
DO 2 J=i,7 
DO 2 1=1,3 
II=< J-DfS-M 
JJ=NCENTIJ) 

POS( in = TABnUT( I 9 JJ» 

2 VELUI )=TABOUTU+3, JJ) 

RETURN 

END 


120700 

120R00 

120900 

121000 

121100 

121200 

121300 

121400 

121500 

121600 

121700 

121R00 

121900 

122000 

122100 

122200 

122300 

122400 

122500 

122600 

122700 

122R00 

122900 

123000 

123100 

123200 

123300 

123400 

123500 

123600 

123700 

123R00 

123900 

124000 

124100 

124200 

124300 


A-S9 



SUBROUTINE ITER ( I R9R 1 * V I 9F 9 K« KK , OT, R , V I 

IMPLICIT REAL*B( A-H9O-Z 9» I 

COMMON C ( 1 000 )9lC(50)9243)9 VV(3)9 NOR 

common/steps/ STEP 9«0L0 

DIMENSION Ki( 3) 9RK 3t8 IcROvB ) 9 VI (39R )« VI 3v AIvRRi 31 vVVVO ) vCB(20) 

COUWTxO.DO 

ZRS*6c965D5 

XK*l 

F=0*DO 

DTxOoDO 

T® 0 Q DO 

LK*l 

TOL D»OoOO 
FOLOxO.DO 
NN-I 

IF4KK,EOo2) XK»1 

C K>1 MEANS EARTH SHADOW « K«2 MEANS MOON SHADOW 
C K»3 MEANS OTHER PLANET SHADOW 
C KKxl MEANS PENIIMBRA9 KK>2 MEANS UMBRA 
DO 2000 1*198 
DO 2000 J-U3 
RUellsRl' 

2000 V(J9n>Vltu, , 

GO TO n9?9lO)9K 

1 DO 3 1*193 

3 XU n»RU9U-Rn93» 

SEVA* ADOUXlvRIUin 
RFSsFNORM(R(l93n 
RFE* FNDRM(R(l6in 

DELS* DARSIN(6500900/RFE)'f 57 *295779500 

ZRE*6500.00 

GO TO ^ 

2 DO 5 1*193 

5 XUl l*R(I »2>-R(I 93) 

SEVA* AOOTIXURt l92U 
RFS* FNORMIRIUBU 
RFE* FNQRM4RU92H 

DELS*DARSIN( 1738 *D0/RFE )*57o 2957795D0 

ZRE*1738.D0 

GO TO 4 

10 DO 20 1 * I9 3 
20 Xim » Rfl9 NQRI “ RU93I 
SEVA* ADOTlXltRIUNORH 
RFS»FN0RM(R4 USH ' 

RFE- FNORM4RU9NORn 

DELS*DARSlNfC1NaR•^20^/RFEI♦57.29 5 7 7 95n0 
ZRE»CiNaR4>20I 

4 CONTINUE 
ZRE-DSIGNiZREfXK) 

SIG»DARS|N| 4 ZRS’t'ZRE ) /RFS I *^^7 « 2957795D0 
F-SEVA+OELS-a-XK^SlG-lSO.DO 
IF 40ABS4F I «LE. 0.00000 100) RETURN 
IFICOUNT.GT.IOO.DO) GO TO 11 
GO TO «697)«LK 

6 U-2 
OT**-STEP 
GO TO B 

7 OT*ITQLO**F-T*FOLOS /« F-FOLD) 
a FOLO-F 

TOLO*T 

T*OT 

CALL ST£F09NN«07»TA9R18 19 IBlvVH 1» I B) 9 CU 8) 9RR9 VVV 9 1 9C8) 

TIM* C « 593 'fe +« X41L0-K; 4 5 06 1 4-T 1 /86400. DO 
CALL TFRAC4C4592)9TIM9TT»TTF) 

CALL SNEFTFnM,RR*VVV9TTeTTF9R»V) 

COUNT«€OUNT-M.DO 


124400 
124500 
124600 
124700 
124800 
124900 
125000 
125100 
125200 
125300 
125400 
125500 
125600 
125700 
125800 
125900 
126000 
126100 
126200 
126300 
126400 
126500 
126600 
126700 
126800 
126900 
127000 
127100 
127200 
1 27300 
127400 
127500 
127600 
1277O0 
127800 
127900 
128000 
12810U 
128200 
128300 
12 8400 
128500 
12B600 
1ZB700 
12'^800 
1 28900 
129000 
129100 
129200 
129300 
129400 
129500 
129600 
129700 
129800 
129900 
130000 
130100 
130200 
130300 
1304 4M) 
130500 
130600 
130700 
130800 
1 30900 





GO TO Uf2a01 vK 

131000 

11 

MRITE46e 12) 

13J100 

12 

FOftMAT(lHO»22HMAX ITERATIONS IN ITER) 

131200 


STOP 

131300 


END 

131<^00 





SUBROUTINE KONSTK 


IMPLICIT f 

DIMENSIUN 

niMENSION 

DIMENSION 

DIMENSION 


FAL*ft (A-H,0-Z 
C( 1000) ,IC(50) 
SIAOs 1)»SSSI AO 96 ) 
DNAME(2»6) 

0(7,3 ) 


C THE FOLLOWING STATEMENT IS) HAVE BEEN MANUFACTURED BY THE TRANSLATOR TO 
C COMPENSATE FOR THE FACT THAT EO»il VALENCE DOFS NOT REORDER COMMON 

COMMON C 

COMMON IC 9 SSS 

C DIMENSION C( 1000) »IC(50) 

C DIMENSION 5(40, 1),SSS(A096) 

C dimension' DNAME(2,6) 

C DIMENSION 0(7,3) 

C COMMON C, IC 9 SSS 

EOU I VALENCE (C ( 201 ) , S( 1 ? II ) j ( C( 341 ) , D( 1 , 1 1 ) , ( C ( 44ft) ,nNAME(l,l) ) 

C DEFINES CONSTANTS NEEDED BY PROGRAM 

C 


DO 1 1=1,1000 
1 C( I)=0,000 

C planetary gravitational constants (KH3/SEC2) 
cm = 39H603,2D0 
C(2> = 4902^77800 
C(3) « til3271545D12 
C(4) = 324760«5D0 
C( 5)=42915.513D0 
C(6)«.126710f 909 
C(7) = 10000*0L 0 
C( fl)=100nOdOLO 

C RADIUS OF SPHL-RE OF INFLUENCE (KM) 

C( 11 )»925000o0D0 
C( 12)=6o000,000 
C( 13)*1.010 
C( 14) = 616000<,000 
C( 15)=565000*000 
C( 16') = *48De 
C( 17 )» 100000.000 
CU8) = 1000 00.070 

C IMPACT RADIUS (KM) 

Ct21)=6378,165D0 
C(22)*1738.0D0 
C(23)=694000. ODO 
C( 24) >6050.000 
C(25) s 3417.000 
C(26)x71OOD.0O0 
r(27)»lOOOO.OOO 
C ( 28 )x 10000. 000 

C MINIMUM TIME STEP FOR OUTPUT (SFC) 

C( 31) >43200. 000 
C(32)>21600.000 
C(33l>1728000.000 
Cl 34) >43200.000 
C( 35) >432 00. ODO 
Cl 36)>86400.0DO 
Cl 37) >43200. 000 
CI38)>43200.000 

C SIDEREAL ROTATION RATES IHAD/SEC) 

Cl4i)>.7292U52D-4 

C INTERPOLATION INTERVALS FOR EPHEMERIS lOAVS) 

CI5U>40«000 
Cl 54) >40. ODD 
Cl 55 )>40.D00 
Cl 56) >100.000 

C NAMES OF EPHENERI5 BODIES 

DATA Q000HL/6H EARTH/ 

Cl6l)>OOOOHt 


131500 

131600 

131700 

131R00 

131900 

132000 

132100 

132200 

132300 

132400 

132500 

132600 

132700 

132ftOO 

132900 

133000 

133100 

133200 

133300 

133400 

133500 

133600 

133700 

133R00 

133900 

134000 

134100 

134200 

134300 

134400 

134500 

134600 

134700 

134fl00 

134900 

135000 

135100 

135200 

135300 

135400 

135500 

135600 

135700 

135R00 

135900 

136000 

136100 

136200 

136300 

136400 

136500 

136600 

1^6700 

136800 

136900 

137000 

137100 

137200 

137300 

137400 

137500 

137600 

137700 

137800 

137900 

138000 


A>» 


nonooooo 


C 

C 


C 


C 


DATA OOOIHL/6H MOOM/ 138100 
C(62 J:=Q001HL 138200 
DATA 0002HL/6H SUN/ 138300 
C<63)=0002HL 138400 
DATA 0003HL/6H VENUS/ 138500 
C(64)=0003HL 138600 
DATA 0004HL/6H MAKS/ 138700 
Cf65)**0004HL 138800 
DATA 0005HL/6HJUPTeR/ 138900 
C<66t=0005HL 139000 
DATA 0006HL/6HBODV 7/ 139100 
C(67)-Q006HL 1392C0 
DATA 0007HL/ 6HBnDY 8/ 139300 
C{681=0007HL 139400 
ASTRONOMICAL UNIT IN KM 139500 


Cl 701 = * 1495D9 139600 
SECONDS PER MEAN SOLAR DAY AND RECIPROCAL 139700 
C( 71)s86400cODO 139800 
C( 72 l = lol574074r)-5 139900 
SEMI-MAJOR AND SFMI-MINQR AXES OF THE EARTH (KMl lAOOOO 
C( 731=6378,206400 140100 
C( 74)=6356o583800 140200 
DEGREE-RADIAN CONVERSIONS 140300 
CI75)=, 01745329300 140400 
Cl 761 = 57,295779500 140500 
CI77I = 1,5707963200 140600 
CI781 = 3,1415926500 140700 
CI791 = 6,2831053000 140000 
DATA 0008HL/6H K/ 140900 
C(85)=0008HL ‘ 141000 
DATA 0009HL/6HM-KM/S/ 141100 
C(86 »=0009HL 141200 
DATA 0010HL/6HEC-OEG/ 141300 
CI87(=0010HL 141400 
OATA 0011HL/6H U/ 141500 
C188 1=0011HL 141600 
OATA 0012HL/6HN1TS A/ 141700 
C189 I*0012HL 141800 
DATA 0013HL/6HU— OEG/ 141900 
C(90 i=30l3HL 142000 
OATA 00i4HL/6H FOU/ 142100 
CI9U-0014HL 142200 
OATA 0015HL/6HATOR 0/ 142300 
CI92I*0015HL 142400 
DATA 0016HL/6HF DATE/ 142500 
CJ93I»0016ML 142600 
OATA 0017HL/6H ECLI/ 142700 
CI94J*0017HL 142800 
DATA OO10HL/6HPTIC 0/ 142900 
C!95I*Q0i8HL 143000 
CI96}*Q016HL 143100 
OATA 0(U9ML/6M S/ 143200 
CS97»s0019Ht 143300 
OATA 0020HL/6HELEN0ti/ 143400 
CI98}«0020HL 143500 
DATA 0021HL/6HRAPH1C/ 143600 
U99I-0021HL 143700 
U 1221=1,000 143800 


THE fOUOWING ARE EARTH POTENTBAL COEFFICIENTS ♦'»'** 143900 


CI200S 

THRU 

C«2071 

IS 

U 0901 

THRU 

Cg7»0S 

144000 

eg 2091 

THRU 

Cg2141 

IS 

Cl If 11 

THRU 

CI79 U 

144100 

CC2151 

THRU 

CI2201 

IS 

ciZfZi 

THRU 

CI7921 

144200 

CI2211 

THRU 

Ci2251 

IS 

Cg3ir31 

THRU 

€^7*31 

144300 

C g 226 1 

THRU 

CT2291 

IS 

Cl4«4l 

THRU 

Cg7»41 

144400 

Cg2301 

THRU 

CI2321 

IS 

Cf 5^51 

THRU 

C«7»51 

144500 

C«233B 

THRU 

Cg234l 

IS 

Cl 6 96 1 

THRU 

C|7«61 

144600 




nnonoriorio 


C(235) 

IS C(7,7) 

144700 

C(236) 

THRU CC243) IS S(0»0) THRU S(7»n) 

144000 

C(244) 

THRU C(250) IS S(l,l) THRU S(7,l) 

144900 

Cl25i» 

THRU CI256) IS SI292) THRU S(792) 

145000 

C{257> 

THRU C(261) IS S1393) THRU S(7»3) 

145100 

C(262> 

THRU CI265) IS S1494) THRU S(794) 

145200 

C(266) 

THRU C(268J IS S(5o5) THRU SI7»5) 

145300 

C4269) 

THRU CI270) IS 6 ( 696 ) THRU S(7v6) 

145400 

C(271) 

IS S(7,7) 

145500 

C(200> 

= 0 

oODO 

145600 

C(201 ) 

=^0 

cODO 

145700 

C(202) 

= - 

o 1082 645000020-2 

145000 

C(203» 

*0 

.25459999990-5 

145900 

C(204P 

Si 

ol64 0-5 

146000 

C( 205> 

s 

.210 0-6 

146100 

C1206) 

= 

-.6459999999930-6 

146200 

CJ 207> 

s 

.3329999999980-6 

146300 

Ct20B) 

= 

0.00 

146400 

C(209) 

Si 

0.00 

146500 

C(2i0) 

a 

.209111899862 0-5 

146600 

t( 2 U > 

3 

-.542646846405 0-6 

146700 

CJ212I 

= 

-.6765155824770-7 

146000 

C(213) 

S 

-.36979402246 0-7 

146900 

C12UI 

3 

.14418923578 0-6 

147000 

C1215» 

3 

.1535637096760-5 

147100 

C(2L6) 

s 

.2507007287360-6 

147200 

C(217J 


.7379024325770-7 

147300 

C<218) 

Sr 

.1021177075260-6 

14 74 00 

C(219) 


.8583830313850-8 

147500 

U 220 » 

C 

.362552678174 0-7 

147600 

C(22U 

= 

.7822771248130-7 

147700 

C(222) 

= 

.5085697732730-7 

147000 

C{ 2231 

s 

-.1717787428880-7 

147900 

Ct224) 

St 

-.119630040940-8 

140000 

CJ225J 

s 

-.352147606138 0-8 

140100 

C(226) 

s 

-.11190293875 D -8 

140200 

C(227) 

s 

-.2063363282040-0 

140300 

U228) 

3Z 

-.1665608129460-9 

14 06 00 

a 229} 


-.3227765545750-9 

140500 

C(230f 


.3841089461880-9 

14 0600 

C«231) 

K 

-.252611520257 0-9 

140700 

C«232) 


.2689804621530-10 

140000 

C(233> 

■ 

-.9319190365590-11 

148900 

Cf234( 

X 

- 0 1450663738260-10 

16 q 000 

C4235} 

3 

.1020279453390-11 

149100 

C(236) 

3 

0 . DO 

1^9200 

C(237> 

V 

0.00 

169300 

C(238» 

■- 

0.00 

149600 

C(239S 


0,00 

149500 

C«240) 


0.00 

169600 

C4241) 

» 

0.00 

169700 

Cl 242) 

« 

0.00 

169000 

C1243I 

S 

O.DO 

169900 

C1244) 

m 

0.00 

150000 

Cl 245) 

m 

0.00 

150100 

Cl 246 ) 

a 

•2873128376320-6 

150200 

CI247) 

a 

-.4449324667790-6 

150300 

Cl 248) 

a 

-.86203930746 0-7 

150600 

CI249) 

m 

-.21243466397 0-7 

150500 

CI2S0) 

m 

•114180306536 0-6 

150600 

CC25i) 

m 

-.872066750116 0-6 

150700 

CI252) 

a 

-•183761983737 0-6 

150800 

CI253) 

a 

•1478040933150-6 

150900 

U254) 

a 

-.3754565474840-7 

151000 

CI255) 

a 

-.455316216646 0-7 

151100 

U256) 

a 

•162351886106 0-7 

151200 





o o o o o o n o o o o nooonoooooo 


• V- 


C1257) = ,22589R207164 0-6 
C<258) = 1135467178870-7 

Cl 259) = ,231240615425 D-9 
C126U) = ,6427505790650-9 
Cl 261) = ,2535462764220-9 
C1262) = ,4859636964750-8 
C1263) = ,498321698304 0-9 
C1264) » 1960875025120-8 

C1265) = -,2166000563620-9 
C1266) = -,1457644206050-8 
Cl 267) = -,3696360264510-9 
C1268) = ,1911176967930-10 
C1269) = -,3611186266600-10 
C1270) = ,4372814139240-11 
C1271) = ,1780851409540-11 

THE FOLLOWIWG ARB MOON PUTBNTIAL COEFFICIENTS 
CI272) THRU C1276) IS C10,0) THRU C14*0) 

C1277) THRU CC2R0) IS Cll»l) THRU C14,l) 

C1281) THRU C(2B3) IS CC2»2) THRU C(4,2) 

C(204) THRU C1285) IS C(3»3) THRU C14,3) 

C1286) IS C(4v4) 

C1287) THRU C(291) IS S1090) THRU S1490) 

C1292) THRU C1295) IS Sll«l) THRU S(4vl) 

C1296) THRU Cl 298) IS S12»2) THRU S14»2> 

C1299) THRU C1300) IS S13,3) THRU SI4»3) 

C1301) IS S14,4) 

C1272)«0,0D0 
C1273)»O,0DO 
C(274)=-,207108D-3 
C1275)=0, 210-4 
C(279)=0,34D-4 
Cl 281)»0, 207160 -4 
C1284)*0,025830-4 
C1462)* FLAG FOR IMPACT 
Cl 462) a 1,00 
Cl506)=0oD0 

C(719)= REDUCTION FROM UNIVERSIAL TO EPHM TIME 1966 
Cl 719)*40.0D0 

C1720)a RFTIFICATION LIMIT IN ROUTINE OERIV 
C(720)a0. 02500 

Cl725)a J2 OF MARS 
Cl725)a0, 194666 70-2 
C(726)a J3 /F MARS 
C(7Z6)«0;D0 

Cl 758) STANDS FOR TYPE OF INPUT 

«1 EARTH MEAN EOUATOR AND FOIJINOX OF 1950 
■2 TRUE EQUINOX ANO ECLIPTIC OF DATE 
■3 SELENOGRAPHIC TRUE MOON EQUATOR OF DATE 
*4 MARS MEAN EQUATOR OF DATE 
*5 TRUE EARTH EQUATOR AND EQUINOX OF DATE 
Ct 758)«1.D0 
C Cj76U IS DAYS FOR INTEGRATION 

Cl 761 )a|80o0D0 
C1762)al,0D0 
Cl 763) *0,000 
C( 764)«864O0«D0 
Cl 765} >-1,000 
C1766)>1«0D0 
Cl 767) >5, 000 
C1768)»UU00 

C Cl 769) SOLAR PRESS KEY 0,0«00 NOT INCLUDE l,0» INCLUDE 

C(769)>0*000 

C Cl 770) AREA FOR SOLAR PRESSURES ICM,*«2) 

C1770)>22225,00 

C Cl 771) MASS FOR SOLAR PRESS ICM,) 

Cl 771) >113000,00 


151300 

151400 

151500 

151600 

151700 

151800 

151.900 

152000 

152100 

152200 

152300 

152400 

162500 

152600 

152700 

152800 

152900 

153000 

153100 

153200 

153300 

153400 

153500 

153600 

153700 

153800 

153900 

154000 

154100 

154200 

154300 

154400 

154500 

154600 

154700 

154800 

154900 

155000 

155100 

155200 

155300 

155400 

155500 

155600 

155700 

155800 

155900 

156000 

156100 

156200 

156300 

156400 

156500 

156600 

156700 

156800 

156900 

157000 

157100 

157200 

157300 

157400 

157500 

157600 

157700 

157800 


A-SS 



ooorjooo no 


C C<772) «EFLeCTIVITY 

C(772»=Oo2DO 

C C(773I FRACTION OF PARTICLES REFtFCTEO 

C<773)=1«000 

C C(77A» SOLAR PRESS AT ONE AoU* 

C( 77A)*=4*7n-5 

C(775) IS KEY FOR INPUT^’^O. INPUT IS CARTFSIANo^lo INPUT IS 

ORBITAL ELEMENTS 
C( 7751=0,00 

— — C(776) IS KEY FOR STOPPING AT APOGEF OR PERIGEE, 

*=0, DO NOT STOPv =1, STOP AT APOGEE* =-l. STOP AT PERIGEE 

Cl 776)o0,D0 

UV?7f IS THE KEY FOR EXTRA QUTUUT AND IS IN THF FORM APCOEFGH, 

0 MEANS DO NOT PRINT AND 1 MEANS PRINT, 

A= CENTRAL ftOOY MEAN EARTH EQUATOR OF 1950 

B= MOON FIXED 

C= MARS MEAN EQUATOR OF DATE 

G= CENTRAL BODY TRUE EARTH EQUATOR OF DATE 

H=EARTH FIXED LATITUDE AND LONGITUDE 

Cl 7771 = 10000011 , 

C C<901) IS DELTA T FOR UPDATING TRANSFORMATIONS IN HARMON 

C(90U = 15,D0 

C THESE CONSTANTS WERE DEFINED IN MAIN 
IC( 1)=0 
ICI2I=0 
RETURN 
END 


1579D0 

ISflOOO 

158100 

158200 

158300 

158400 

158500 

158400 

158700 

158800 

158900 

159000 

159100 

159200 

169300 

159400 

159500 

159900 

160000 

160100 

160200 

160300 

160400 

160500 

160600 

160700 

160800 





SUBROUTINE Ml a,A#Y) 
IMPLICIT REAL*H{ A-H,n-Z 9 S I • 
C SUBROUTINE GIVES Y^^AX 

DIMENSIUN X (3) ,Y<3 I » A( ^9?S I 
C DIMENSION X(3),Y<3)9 A13*3) 

DO 1 1*1,3 
Yt I JsOoDO 
DO 1 J*l,3 

1 Y( I )*V(I )+AU , J)«4«X( Jl 
RETURN 
END 


160900 

161000 

161100 

161200 

161300 

161^00 

161S00 

161600 

161700 

i6iaoo 

161900 


A-37 
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I. 
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SUBROUTINE MNA(TlMEflOM,CRpDT«ePSIL»Rn,G9GPvWW,FM) 
IMPLICIT REAL*BI A-H,0-ZjS) 

DIMENSION EM(3,3) ,0F13) 

COMMON/XLIB/ SG»TA 


D s TIME 
T - 0/36525»D0 
T2 = T*T 
T3 a T2*T 
Aal3g06^992D0 
DO 6 I = U3 
DDaO^DO 
OOaD 

OOaOD*( A/360,001 
00*:0D-IDINTIDDI 
DF( IlaDO 
GO TO 1^,5,6)»I 

4 As, 985600500 
GO TO 6 

5 A=,1643586D0 

6 CONTINUE 

G=215, 5401300+360, OOfDFt 1) 

Gps358o00906 700+360, DO*nF< 2) 

WW= 196 ,74563200+360, 00 *r>F ( 3 ) 

G a G»if ,01 745329600 
GP = GP#, 01745329600 
WW a WW4,01745329600 
YN = l,53500*o01745329600 

RO a -,0297222D040C0SIGI+«0l0277700«PCnS(G+2,D0*WW) 

I -,0030555500«OCOS I 2,00*G+ 2,ri04ww | 

TA a -,00333300#DSIN(GI+,0163888D040SIN(GP) 

1 +,00500#OSINI2,00*WWI 

SG a -,030277700«OSINIG) + ,0102777O0*0SIN( G+2,00*WWl 
1 -,00305555D0*DSIN ( 2,004G+ 2,00*WW ) 

SG s I SG»,01745329600)/0SIN <YN) 

RO s RO*, 01745329600 
TA « TA*, 01745329600 
YN s VN + RO 
RO a OM + SG + DT 

C! a OCOS (R0»*0SIN <EPS1U*DSIN lYNl 
1 +OCOS tEPSIU*OCOS (YNI 
SI » 1.00 - CI**2 
SI a OSORT (sn 
SO a -OSIN (ROJ*OSIN IYN)/SI 
CO a 1,00 - S0**2 
CO a OSORT iCO) 

SO a -OSIN «ROI*OS!N lEPSIU/SI 

CO ■ «OS!W tR0I*$O*DCOS JEPSIU » DCOS IRO)*CO 

OL*ARKTNS<OeCOsSO» 

8 CONTINUE 

CA • OL ♦ CCR + TAJ - lOM + SGI 

SA ■ OSIN 4CA1 

CA a OCOS ICAI 

RO * OCOS «ROI*OS!N I gPSIU /( SI *COI 

EMUsU a CA*CO - SA*SO*CI 
ENU»2I a CA«$0 + SA*CO*C! 

a SA*SI 

EMUai ■ “SA*CO - CA*SO*CI 
ENI 2 t 21 * <*SA*SO + CA*CQ*CI 
ENUfll a CA*S 1 
EMCSai a SO*S! 

EMI 3.21 a “CO*SI 
EMI3.3I a Cl 
RETURN 
END 


162000 
162100 
162200 
162300 
162400 
162500 
162600 
162700 
162800 
162900 
163000 
163100 
163200 
163300 
163400 
163500 
163600 
163700 
163800 
163900 
1 64000 
164100 
164200 
164300 
164400 
164500 
164600 
164700 
164800 
164900 
165000 
165100 
165200 
165300 
165400 
1655(0 
165600 
165700 
165800 
165900 
166000 
166100 
166200 
166300 
166400 
I6o50(» 
166600 
166 700 
166800 
166900 
167000 
167100 
167200 
167300 
167400 
167500 
167600 
167700 
167800 
167900 
16R000 
168100 
16R200 
168300 
168400 
168500 
168600 


A'»S8 


SUBROUTINE MULT(A,H9C) 168700 

IMPLICIT ReflL*8< A-Hjfl-Z g* I 168800 

DIMENSION A(3,3Igfl(3g3IgC(3g3) 168900 

C DIMENSION A(3g3),Ri3g3lgCl 3g3) 169000 

DO 1 lrl,3 169100 

DO 1 J=l93 169200 

C( I g J) a 0,D0 169300 

on 1 K=lg3 169A00 

I C(IgJ) = C(IgJ) + AU g Kg J) 169500 

RETURN 169600 

END 169700 


A-f» 
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SUBROUTINE NUTAIT ITIME »nH ,CR * DT , EN, F PSI L ) 

IMPLICIT RE4L>«‘R ( A-HjO“Z» 4) 

OIMENSIUN EN(3o3) 

DIMENSION EN(3»3) 

007=0oD0 
D = TIME 
T = D/36525oOOO 
T2 = T*T 
T3 » T2«T 

OM= 12, H2790D0-oD 5295 3922D0*D+ o0020795D0*T+o0020fll00D0#T2+ 0OOOOO2 
ZOODO^TS 

CR a 6^0 37545200-H3, 176397DU‘i'D-o001l315 7Sn0*T-,00113015n0*T2 + o0000 
2019OO*T3 

GP a 208oB4399D0+, lll4040HD0»9‘D-o0l033400n0>«‘T~,01034300n0«T2-,0n00l‘ 
220000*13 

VL s 28O»0fll21DO+,98564734D0*O-»-,00030300nO*I T<-T2I 

G = 2ft2o0R053D0+, 4706840-04*0+, n0045S25no*T+, 0004«>750DO*T2+ ,000003 
200D0*T3 

OM a 0M*,017453296D0 

CR = CK*, 01743329600 

GP s GP*, 01745329600 

VL = VL*,01743329600 

G a G*,01745329600 

OEl = 25,58440000 * DCOSIOM) 

DE2 = -.251100 * OCOSl 2,D0*0M) 

DE3 a l,533600*OCnS( 2.00*VU 
DE4 a .066600* OCOSI 3,00* VL -G ) 

0E5 a -,025800 * OCOSf VL+G ) 

0B6 a - ,018300 * DCOSf 2,00*VL-nM) 

0E7 = -.006700 * OCOS( 2,00*GP - OM ) 

DE * OEl + 0E2 + 0E3 + DE4 + 0E5 + DF6 + 0E7 

001 a .245600 * DCUSl 2,00*CR) 

002 a ,050800 * DCOSI 2.00* CR - OM) 

003 a. ,036900 ♦ DCOStB.OO* CR - GP I 

004 a -.013900 * DCOSI CR+GPI 

005 a -.008600 * OCOSI CR-GP+OM ) 

006 a ,008300 * DCQSt CR-GP-UM \ 

008 a .006100 * DCOSI 3,00 * CR+GP - 2,00*VL I 

009 a .006400 * OCOSI 3,00*CR -GP -OM I 

00 a DOl + 002 + 003 + 004 + 005 + 006 + 007 + 008 + 009 
DTI a -{ 47.89270000 + .O402DO*T) * DSINIOM) 

0T2 a ,5800 * DSIN( 2,D0*0M ) 

0T3 » -3,536100000 • DSINI 2,00 • VL I 
DT4 a -.137800 * OSIN( 3.D0*VL -G I 
DT5 a .059400 * OSINI VL+G I 

0T6 a .034400 * OSSNl 2,00*VL - OWI 

0T7 a .012500 * DSINI 2,00*GP - OM J 

OTfl a .3500 • OSINI VL-G » 

0T9 a .012500 * 0SIN( 2,D0*VL - 2,00*GP I 

OT a OTl + DT2 + 0T3 + DT4 + f)T5 + DT6 + 0T7 * 0T8 + DT9 

DSl a -.565800 * nSINI 2,00*CR J 

DS2 a -.09500 ♦ OSINI 2,00*CR - OM > 

DS3 ■ -.072500 • OSIN( 3,D0*CR - GP J 

054 a ,031700 • OSIMI CR+GP I 

055 a .016100 ♦ OSINI CR -GP +0M 1 

DS6 a .015800 * OSINI CR -GP -OM ) 

0S7 a -.014400 * OSINI 3,D0*CR + GP - 2*00*VL » 

0S« a -.012200 ♦ DS!N« 3.00*CR - GP - OM T 

0S9 a ,187500 « DSINI CR - GP 1 

DSIO a ,007800 * OSlMf 2,00*CR “ 2,00*GP S 

OSLi a *041400 * OSIMI CR+GP - 2«00*VL I 

DSi2 a .016700 * OSIMI 2«00*CR - 2,00*Vt 1 

0S13 a -.008900 * OSINI 4,00«CR - 2,00*VL I 

OS a DS1+052+OS3+OS4+OSS+OS6+DS7+OS8+OS9+OSIO+OSU+OS12+OS13 

OE a ,174532960-5* lOE+OO) 

OT a ,1 74532960-5* «OT+OSJ 


169800 

169900 

170000 

170100 

1 70200 

170300 
1 70400 
170500 
170600 
170700 
170800 
170900 
17 1000 
171100 
171200 
171300 
171400 
171500 
171600 
171700 
171800 
171900 
172000 
172100 
172200 
1723O0 
172400 
172500 
172600 
172700 
172800 
172900 
173000 
173100 
173200 
173300 
173400 
173500 
173600 
173700 
173800 
173900 
174000 
174100 
174200 
174300 
174400 
174500 
174600 
174700 
174800 
174900 
175000 
175100 
175200 
175300 
1.75400 
175500 
175600 
175700 
175800 
175900 
1 76000 
176100 
176 ^'- 

f o300 





ea » 23o 

4^575fl7D0-*o0130«J4040U*T-,000000flRn0*T2-«‘,000000Sn0*T3 

1 76400 

EB » EB* 

*01745329600 

176500 

EPSIL w 

EB+DE 

1 76600 

EN4 i«n 

= i.ooc 

176700 

6M( l?2)* 

»DT*DCOS(EB) 

176B00 

EN( l93»x 

-0T*DS1N(EB> 

176900 

F.N(2j1 ) 

= *-EN( 1»2> 

177000 

EN( 2f2) 

= UODO 

177100 

Ef^4293» 

= -DE 

177200 

EN( 3^1 ) 

= -EN(l93) 

177300 

EM(3f2» 

a OE 

1 77400 

Ef^( 3»3) 

a 1»ODO 

177500 

RETURfS! 


177600 

END 


177700 
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SUBROUTINE ORBf^sOXjUtNOUT) 
implicit REAL*fie A-H^0“Z*$I 
DIMENSION CUOQO) 

DIMENSION K(3) 9D»(3UBI3) 

C THE FOLLOWING STATEMENT! SI HAVE BEEN MANUFACTUKED RVTHE TRANSLATOR TO 

C COMPENSATE FOR THE FACT THAT EQUlVALtNCF DOES NOT REORDER COMMON 

COMMON C 
COMMON C 

DIMENSION C(IOOO) 

DIMENSION XOIvDXUIfBn) 

eou 3 VALENCE tC (700) j SMA) , ! C( 701) ? FCC ) ? ( C ( 702 ) nO INC ) »(€( 703) »OMG) 

1* (C(704) 9 BEP) 9 (C(705) ®RCA) 9 ( C( 7 06 ) » C3 ) » ( C( 707) 9 THET) «( C(70R ) 9 PERV) 
2o(C(709),P)9 (C( 710) 90 b )9 (C( 711 ) 9 TPFR ) 

CALL CROSS OtflOXflf?) 

R2 * DOT(X*X) 

R=DS0RT(R2) 

V2 - OOnDXpDX) 

B2 « OOT(B,B) 

BB»DS 0 RTf B2) 

A*00T(X9DX)/U 
P * B2/U 
C3 = V2“2o0F>0*U/R 
SMA»“U/C3 

ECC=DSORT(DABS( l»0D0+C5«P/U ) ) 

0INC»ARKTNSnB09B4 3)9OSORT(B( 1 )>^« 2 '«'B( 2 2 ) ) 

OMGa ARKTNS(360»-B( 2) 98( 1 ) ) 

RCA*P/U,ODO+ECC) 

THTaARKTNS(3609 (P-R)9Bfl*A) 

BETa ARKTNS(3609X(2)*B( 1 ) -X( 1 ) #P ( 2 ) 9 X ( 3 )*BB) 

BEPa BET-THT 
IF(BEP) 29393 

2 BeP=BEP+6.2ft31«53D0 

3 CONTINUE 
RT0*57. 2957795100 
OINC*OINC#RTO 
OMGaOMG*RTD 
BEP=BEP*RTO 

PERVaDSORT ( DABS (C3-^2o ODO*U/RCA ) ) 

VIMP*PERV-OSORT(U/RCA) 

CTASa(P/R-l .ODO)/ECC 
IFtOABS!CTAS)»lT.l,DO) GO TO 200 
STASs0,00 

CTAS-DSIGNf 1 .OOOpCTAS) 

GO TO 201 

200 STASaDSORTt K0D0-CTAS*CTAS) 

STAS a OSIGNt STAS 9 A ) 

201 CONTINUE 

THE=ARKTNS« leOvCTASpSTAS) 

THET*THE*«TO 

CALL TC0NIC(U*ECC95MA;P9THE9TPER9FAC) 

TPER«TPER/fl6400»0D0 
IFISMAIiOalOfll 

10 WRITE « 6 f 7)SMA9ECCp0INC90MG9BEP9RCAeC39THET*PERV9 P 9 VIMP 9 TPFR 
7 F0RMATI5H SHA015*B»5H ECCDlSoBsSH lNCDl5.ft»5H LAN015.B9 

15H APFOlS.epgH RCAD15»8/5H C3D15*Bp5H THFTDlS.fipSH PERVDlS.flp 
2SH SLRDlgpApSH IHPV015 .SpSH TPFRD15.B) 

GO TO 12 

11 F1>A*U/0S0RT1U«SMA) 

F2-i.0D0“R/SMA 
S1NE»F1/ECC 
C0SE>F2/ECC 

E«ARKTNS ( 360« COSE « S ENE ) 

XNANa I E*ECC *S I NE } *R TO 
PERM « FAC * CI79) / 3600,00 
AFBK*$MA« i I . DO’t'ECC ) 

WRITE « 695)SMA9ECC«0INC»UMG«BFP9RCA9C39THFT9APBKp XMANpPERHpT 


177B00 
177900 
17B000 
17B100 
17B200 
17B300 
1 7B400 
17B500 
17B600 
178700 
178800 
178900 
179000 
179100 
179200 
179300 
1 79400 
179500 
179600 
179700 
179800 
179900 
180000 
180100 
180200 
180300 
1 80400 
180500 
180600 
180700 

180a00 

180900 
181000 
181100 
181200 
181300 
181400 
181500 
181600 
181700 
181800 
18 1900 
182000 
182100 
182200 
1P2300 
182400 
182500 
182600 
182700 
182800 
162900 
183000 
163100 
183200 
183300 
183400 
183500 
183600 
163700 
183800 
183900 
104000 
184100 
184200 
184300 


A-4S 



IPER ia*400 

5 FnRMAT(5H SMAOl5.8j5H eCCDl5«,8,5H iNCniSoF^SH LAND15o8» 18A500 

15H APF0l5.»«fl5H ftCAD15*B/&H C3ni5*Rp5H THETrn8,fl,5H APncUlSoR* 184600 

25H MAN015o8p5H PFRHD15.ft?&H TPFR015»8) 184700 

12 CONTINUE 184800 

RETURN 184900 

END 105000 


A-49 



SUBROUTINF ORB2X()(,VX»E’ U) 

IMPLICIT REAL*fl( A-H,0-Z 
DIMEMSICJN X(3)»VX(3l,eU6)jA<:is 21 
C EL CONSISTS OF INPUT ORBITAL FLFMFnTS 

C m SEMI-MAJOR AXIS (+ IF ELLIPTIC? -FOR HYPERBOLIC) 

C (2) ECCENTRICITY (LESS THAN 1» ELLIPTIC? GREATER HYPFRBOLTC) 

C (3) TRUE ANOMALY (OEr7RFES) 

C **#IF C379 CHANGED TO i*?' MEAN ANOMALY IS USED FOR INPtfT 

C379 = 0?D0 

C (A) LONGITUDE OF THE ASCENDING NODE (DEGREES) 

C (5) INCLINATION (DEGREES) 

C (6) ARGUMENT OF PFRIAPSIS (DEGREES) 

DATA- UR/ *01745329300/ 

CnM*>0C0S(DR“«‘EL(4)» 

SOM=DSIN(DR>)‘EL (4) ) 

C IN=DCQS( DR*EL (5) ) 

SNI = DSIN(DR*FL (5i ) 

IF(C379 )10?10920 

20 XM=FL(3)*DR 
ECC=EL(2) 

EsXM+eCC*OSIN (XM)+(ECC**2/2»D0)*DSIN { 2oD0*KM)+( ECC** 3/24 *D0 ) * ( 9 . 
lf)0>»0SIN ( 3*D0*XM)-3oD0*DSlN ( XM) ) -H ECC#*4/ 192*00 )«( 64 oDO*DS IN (4oD 
«0*XM)-32oD0*0SIN (2oOO«XMl) 

CTA:=(OCUS (E)-ECC)/( loDO'ECC^DCnS (E)) 

STA=DS(JRT ( 1,D0-ECC#*2)#DSIN ( E ) /( 1 * DO-FCC*DCOS (EM 
TAA=AKKTNS( 180?CTA?STA) 

GO TO 21 

10 CTA=OCOS (0R>>EL(3I) 

STAa0SlN(DR#EL(3)) 

TAA=EL(3)«DR 
21 TBB=TAA+EU6)’>DR 
CBA=DCOS(TfiR) 

SBA=0SIN(TBB) 

A( l?UaCOM*CBA-( SOM'^'CINJ'XSBA 
A(2»l )aSOM*CBA+(COM*ClN)*SBA 
A( 3? 1 )irSNI«SBA 

A( l»2)*-COM#S8A--(SOM»CIN)*CBA 
A( 2?2 )=-SnM*SBA-MCOM*CIN)*CBA 
A( 3,2)*SNI>«'CBA 
P*-EL( 1 )*( loOOO-FU2)’*'*2) 

R=P/( loODO+EL (2 >*CTA) 

C R IS THE RADIUS VECTOR MAGNITUDE 

VRsEL(2 )*STA>»DSORT{U/P ) 

C VR IS THE RADIAL COMPONENT OF VELOCITY 

VT=DSQRT(U#(2o0D0/R-1.0D0/EL( 1) )-VR*VR) 

C VT IS THE TANGENTIAL COMPONENT OF VELOCITY 
on I I«l?3 

X(I)sR*A(I?U ^ 

VX( n-VR*A(U 1)+VT*A( I 9 2) 

1 CONTINUE 
RETURN 
END 


185100 
185200 
185300 
185400 
185500 
185600 
185700 
105ROO 
185900 
186000 
186100 
186200 
106300 
I 86400 
186500 
186600 
186700 
1 86R00 
186900 
187 000 
187100 
187200 
187300 
187400 
187500 
187600 
187700 
187800 
187900 
188000 
188100 
188200 
188300 
188400 
188500 
188600 
188700 
188800 
188900 
189000 
189100 
189200 
189300 
189400 
189500 
109600 
109700 
189800 
189900 
190000 
190100 
190200 


A<-44 



SUbROUTIME OUKXpVflO^ERRORTNffLff H) 
IMPLICIT REAL*fif A-H9O-Z (.$ > 

[DIMENSION AI39 3I,EN(3»3)*EMN(39 3»9AR( 
DIMENSION 0ATE(393)9SDATE(3»3I sOATETI 
1ADATET(3*3) »DATEEC( 39 3 WDATECK 3931 9 X 
DIMENSION Y(20»9Dr20l9ERR0R( 30) 
DIMENSION RAT6(3)9PERT(3)9 PERTTI3) 
DIMENSION DACC(3) vCACCI 3) tOACC(3) 
DIMENSION RH39B)9V1(39B) 


DIMENSION XTEHP(3)9 VTEMP(3) 
DIMENSION DUMI393) 
COMMON/DISTRB/DACC 9 CACC 9 OACC 
COMMON /ROT/ DATE 9 SDATE 9 DATET 9 S DATBT 9 
C0MM0N/TIM6/TW.TF 

COMMON C(1000)9lC(5D)9Z(3)9VV(3)9N0R9 
COMMON/TUP/TLli 

COMMON/PRINT/ I C 5 O 9 ! SEL 9 MARDAT 9 ICENEC 
DATA ST/-1*/9ST1/-1./ 


393 ) 

393) 9 SDATFT(393) 9A0ATE(393), 
X( 3)9 '(V(3) 


ADATF* ADATET sDATEFCffDATECT 
!NP 

9 IEAREC 9 ISUNEC 9 ICDATE9IEARFX 


EQUIVALENCE ( Cl 40 1 ) 9 R 1 Cl 9 1 ) > 9 ( C ( 425 ) 9 V 1 CI 9 1 )) 

3 TIM sCC 593) + (X-t-C (506) )/fl64009D0 
CALL TFRAC(C(592)9TIM 9 TN 9 TF) 

IF lCa22) oNEo OoDO ) GO TO 3000 
STl = TW 
ST *= TF 

IFCC(594)-tC(506)^X)oGT.0o00) GU TO 50 
3000 IF(TW«EOoSTl„ANn,DABS(TF-ST) qLE*. 0000200) GO TO 50 
DO 1002 JJ=l93 


1002 PERTC JJ)=«OACC( JJ)-CACC( JJ) 
STi=TW 


ST=TF 

WR1TEC699) 

9 FORMAT! IHO) 

CALL DAT0UTCTW9TF90AYM9FDATE90) 

TP*TW<>TF 

TlUsTP 

CALL ROTEOCTPsA) 

CALL NUTAIT(TP9nM9CR9DT9EN9EPSI U 

CALL MULT(EN9A90ATEJ 

DUMClpU-loDO 

0 UMCl 92 )=O*D 0 

DUM4 UBI^O.DO 

OUMC29l)aO»DO 

DimC29 2 [*DC0S(EPSfU 

DUM(2,3)>EDStN(EPSIL t 

DUM(39l)*0-D0 

DUMC392)«-DUMC293) 

0UM(393)i: 0UMC292) 

CALL MULTCDUM9DATE90ATEECI 
DO BOO I»l93 
DO 800 J«U3 

800. DATETC J« n«DATEn 9 JI 


C 

100 IFCICDATE.EOoO) GO TO 30 
CALL MlCRUlpNORUOATEpXX J 
CALL MH vmvNORIpDATEfXVI 
WR1TE(6»48) CCN0R4-60} 

48 FORMATUHOpAStlXpSCENTEREDSSXp’TRUE EARTH EQUATOR AND EQUINOX OF 
IDATE * I 

CALL aUTXCXX,XV«NOUT» 

CALL ORBCXXpXVtCCNORIvNOUTI 
30 IFUC50.EQ«D«AND*C(777}«NE.O«ODS GO TO 31 
WR!TEC6»I0I CCN0R<«^60I 

10 FORMAT! IHO«A691X98HCENT£REO»3X951HCENTRAL BODY MEAN EARTH EQUATOR 
lANO EQUINOX OF 19501 
CALL OUTXCRlUtNORIpVl! i«N0R)9N0UT» 

CALL QRBfRIU»NnR»»VintN0ftl9UNnRl9MnyTI 


190300 
190400 
190500 
190600 
190700 
190R00 
190900 
191000 
19 1100 
191200 
191300 
191400 
191500 
191600 
191700 
191B00 
19 1900 
192000 
192100 
192200 
192300 
192400 
192500 
192600 
192700 
192R0O 
192900 
193000 
193100 
193200 
193300 
193400 
193500 
193600 
193700 
193B00 
193900 
194000 
194100 
194200 
194300 
194400 
194500 
194600 
194700 
194R00 
194900 
195000 
195100 
195200 
195300 
195400 
195500 
195600 
195700 
195B00 
195900 
196000 
196100 
196200 
196300 
196400 
196500 
196600 
196700 
196800 



A-45 


o o 


31 IF( ieARFXoEO*O.ORoWOKoNE«l) GO TO 
ITIMEsTP 

T2=> UP- 1 T I ME H<86A00<, DO 

TlsITIME 

DA-EN42» U 

CALL GHA(T2»Tl,fiHAiNj,DA»nMEGAl 
GHAN=GHAN*C(75) 

CALL MURin,NOR),DATE9XTFMPl 
SINT-DSIN(GHAN) 


200 iVf^TOO 

197000 

197100 

197200 

197300 

197A00 

I97*i00 

197A00 

197700 


COST=DCOS(GHAN) 

1)= XTEMPl n*COST+XTEMP( 2)*SIWT 
KK(2)-“XTEMP( U0SrNT+KTf:MP(2)*CnST 
XKOJaXTEMPO) 

KLONGI = ARKTNS{ ISO® XX ( 1 WXK 1 2 ) )*C<76) 

C 

200 CONTINUE 

iPf ISEL»E0.0onR.NnR,NE<.2l GU TO 300 
CALL MNA(TPf0M9CR9 0TflEPSlL»R0pG9GPjWWp EMN) 

CALL MULTI EMNpOATEpSDATE) 

00 BOl I=l»3 
DO 801 J = U3 

801 SOATETU, nsSDATEd 9 J) 

CALL MURlUpNOR IpSDATEpXX} 

CALL MUVlfUNORlpSDATEpXVI 
CALL MH PERT 9 SDATF 9 PERTT) 

WR|TE(69500) CCNOR-e-GO) 

500 FORMATUHOpAB, 1 X 9 ®FIXED TRUE MOON EQUATOR OF DATEM 
CALL 0UTX(XX9XV9N0UT) 

CALL 0RB(XX9XV9CIN0R) vNOUTI 
300 CONTINUE 

IF(MARUATeEQ.O«OR»NOR,NE.5) GO TO 1001 
CALL ARES(DAVMfFDATE9AR I 
CALL MULTI AR 9 A 9 ADATEI 
DO 802 1=1,3 
DO 802 J=l,3 

802 AOATETIJ 9 I )=AOATE|| , J1 
CALL MltRU I,NOR I 9 ADATE 9 XX) 

CALL MlIVU UNORUAOATEbXV) 

CALL MIIPERT9ADATF9PERTTI 
WR!TEI6,700} ClNOR+601 

700 FORMAT! 1H0,A6, IX 935HCENTERED MARS MEAN EQUATOR OF DATE I 
600 CALL OUTXIXX9XV9NOUTI 

CALL ORBIXX9XV,C|NOR),NOUT) 


1001 IOR»NOR 
600 CONTINUE 

lFn0R«E0.3,AN0.!5UNEC*E0.0) 60 TO 32 
CALL MUVm, I0R1»0ATEEC,XV| 

C 

MRITEI 699 OO) CH0R4>60» 

900 FORMAT! IHO 9 A 69 1 X 9 “CENTERED TRUE ECLIPTIC AND EOIHNOX OF DATE d 
CALL OUTXlXXfXVtNOUTt 
C 

1FU0R.E0«N0R} CALL ORBIXXfXVtClNnRUNOUTI 
32 CONTINUE 

IF! EOR«EO«3*AND.NOR*N£*3I 60 TO 16 

1FUOR*EO.U 60 TO &000 

IOR«| 

60 TO 600 

1000 IFUOR.EO«38 60 TO 16 
|F!N0it«£O«3$ 60 TO 16 
EOR-3 
60 TO 600 
16 CONTINUE 

IF«N0fl«NE«l«OR*tEARFX«Ee*0l 60 TO 33 


197800 
197900 
198000 
198100 
198300 
198600 
198500 
198600 
198700 
19R800 
1 98900 
199000 
199100 
199200 
199300 
199600 
199500 
199600 
199700 
199800 
199900 
200000 
200100 
200200 
200300 
20 O 6 OO 
200500 
200600 
200700 
200800 
200900 
201000 
20U00 
201200 
201300 
201600 
201500 
201600 
201900 
202100 
202200 
202300 
202600 
2025O0 
202600 
202700 
202800 
202900 
203000 
203100 
203200 
203300 
203600 
203300 
203600 
203700 
203800 


o o o 


WRITE(6,3001) «LATl9XL0W6I ?U3900 

33 CONTINUE P04100 

WftITE16»2nCi46n 20A200 

21 FORMAT! 20H0WUMB6R OF STEPS IS D20,8) 20-!>300 

204400 

204*500 

204600 

50 RETURN 204700 

end 204A00 
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SUBROUTINE OUTX U , VX ,NOUT ) ?04900 

IMPLICIT REAL«ft( A-H,0“Z*4) 205000 

DIMENSION XOIoVXOI pCf6) 205100 

OIMENSION X(3 I 9VX(3)fC{6) 205200 

CALL RVOUTIXjVXffC) 205300 

WRITE I 6fl5)K9VX,C 205A00 

5 FORMATISN XDl5,8e5H YOl5ofli^5H ZD15*ft?5H DXD15ofl» 205500 

1 5H DVD15,8j5H DZD15U/5H R015»flfr5H DEC015o8,5H RAD15»R, 205AO0 

2 5H V015.895M PTH015*8,5H AZD15o«» 205700 

RETURN 205A0O 

END 205900 




SUBROUTINE OVRLAY 

206000 


IMPLICIT REAL*R( A-HsO-Z 9$ ) ■ 

206100 


COMMON SnOOO)vIC(5( 

206200 


DIMENSION 6UFP{6»9lN0j6) 

206300 

16 

E0RMAT(6n39D1298) t 

206600 

996 

FORMATl IHO) 

206500 

995 

F0RMAT(6(5X9l3,D1798I i 

206600 

996 

FORMAT! « 1 S60X » 13H0VERLAY INPUT) 

2 06700 


K--1 

206ROO 

15 

CONTINUE 

206900 


REA0(5fll69END*20) (iNOdljBUFPt n»I«l 96 ) 

207000 


DO 10 N=l 96 

207100 


N«N 

207200 


GO TO 111912912, 12)«N 

207300 

11 

CONTINUE 

207600 


IF(IND(NI )139l3,lT 

207500 

17 

IFlK)l8,199l9 

207600 

18 

WRITE(6,996> 

207700 


K^O 

207ROO 

19 

CONTINUE 

207900 


WRlTE16,995f 1 I NO ( J) , BUFPIJ ) , Js 1 ,4) 

208000 


GO TO 16 

208100 

12 

IFUNDIN)) 15,15,16 

208200 

16 

CONTINUE 

208300 


J= INDIN) 

208600 


S( J)=BUFP(N) 

208500 

10 

CONTINUE 

208600 


GO TO 15 

208700 

13 

CONTINUE 

208800 


WRITE(69996) 

208900 


RETURN 

209000 

20 

STOP 

209100 


END 

209200 





c 

c 

c 

c 

c 

c 

C*>f 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c** 
c ♦ 

c 

c 

c 

c 

c 

c 


SIJHROUTINF RtAORl JEDjTSEC* lERRJ 
IMPLICIT REAL^Rf A-HjO-Z J 

READE, read, IWTFRPDLATF 9 TKAWSLATE JPL EPHEMERIS 

JoE* EKELUMD* MESA SCIENTIFIC CORP», 1965 SEPT 15 
C*L« LAWSONt JPL 9 1966 MAR 17 

READ JPL EPHEMERIS AT THE JULIAN EPHEMfiRlS DATF 
GIVEN BV UE0>e-TSEC/86«00oD0) 

ITEMS COMMUNICATED THROUGH THE CALLIWR SEOUENCF ** 

JED REFERENCE JULIAN EPHEME'<IS DATEo 

TSEC SECONDS OF EPHRMEKIS T'Mp PAST JFD« 

ANY COMBINATION OF VAlJFS OF JED TSEC 

IS ACCEPTABLE AS LUNG AS (JFD+TSEC/ R6400 cDO ) 

IS WITHIN THE RANGE OF THE FPHFMFRIS TAPE 
BEING USFO* HOWEVER TO OBTAIN THE 
FINEST POSSIBLE RESOLUTION IN INTERPOLATION 
THE NUMBER JED MUST BE AN EXACT MACHINE 
NUMBER o FOR EXAMPLE JED COULD BE A DATE ENDING 
WITH .0 OR -5 , 
lERR ERROR FLAG 

0=NO ERROR 

1 = { JE0+TSeC/t>o400.D0) LESS THAN FIRST DATE 
ON TAPE 

2»{ JE0+TSeC/86^00«D0) GREATER THAN LAST DATF 
ON TAPE 

3»S0ME IREOm IS NOT 0,1» OR 2 
4sICENT IS NOT IN THE RANGE I THRU 11 
5«ICW IS NOT 1«29 or 3 

THE FOLLOWING ITEMS ARE INPUT THROUGH COMMON 
COMMON BLOCK CETBLl * 

AO A,U. EXPRESSED IN DESIRFD OUTPUT UNITS 

RE EQUATORIAL RADIUS OF BARTH IN DESIRFD OUTPUT UNITS 

RE IS USED TO SCALF THE LUNAR FPHFMFRIS 
TPD DESIRED N(JM6ER OF TIME UNITS PFR DAY 

C0MM0N/CETBL2/lCWfl KENT » IRE0I13I 

COMMON BLOCK CETRL2 * 


C 

C 

C 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ICW 


ICENTR 


IREOI i 


FLAG INDICATING STATUS OF COMMON BLOCKS REC2 ANo CETBL3 

1 MEANS NEITHER BLOCKS CONTAIN VALID DATA 

2 MEANS ROTH BLOCKS CONTAIN VALID DATA 

3 MEANS REC2 IS VALID» CETBL3 IS NOT 
USER MUST SET ICW»1 BEFORE INITIAL CALL 
SPECIFIES CENTRAL BUDY FOR COORDINATE 


TRANSLATION AS FOLLOWS 

1 MERC 5 JUP 9 PLUTO 

2 VENUS 6 SAT 10 SUN 

3 EARTH 7 URANUS 11 MQON 

A MARS 6 NEP 

IREOUl SPECIFIES OUTPUT DESIRFD FOR 
BODY NO. J. 

IREOUl-0 NO OtjTPUT 


1 POSITION 

2 POSITION ANO VELOCITY 
J RUNS FROM 1 TO U AS FOLLOWS.. 

I MERC 5 jyp 9 PLUTO 


2 VENUS 

3 EARTH 

4 MARS 


b SAT 10 SUN 

7 URANUS n MOON 
fi NEP 12 ERTH-MN-BARYCENTFR 
13 NUTATION 


C * COMMON BLOCK CETBL3 * 

C TABS 629 OOU^E PREC* UORO 6UFFEK TO ACCOMODATE J*D. AND EPHEMERIS. 
C NUTAT 204 SINCCE PREC. WORD BUFPEft TO RCXOMODATE Nt^TAT ION DATA. 

C CKSUM 1 S.P. WORD FOR CHECKSUM. 

C*« THE FOLLOW INC ITEMS ARE OUTPUT THRO«lGH COMMON «* 


209300 
209400 
209500 
209600 
209700 
209800 
209900 
210000 
210100 
210200 
210300 
210400 
210500 
210600 
210700 
riOBOO 
210900 
211000 
211100 
211200 
211300 
211400 
2U500 
211600 
211700 
211ftOO 
211900 
212000 
212100 
212200 
212300 
212400 
212500 
212600 
212.700 
212RO0 
212900 
213000 
213100 
213200 
213300 
213400 
213500 
213600 
21^700 
213B00 
213900 
214000 
714100 
214200 
214300 
214400 
2145O0 
214600 
214700 
214R00 
2 14900 
215000 
215100 
215200 
215300 
215400 
215500 
215600 
215700 
215R00 


A-IO 


on n on non oooooo on o o oo ooooooooonoono 


# CIIMMON bLUCK CETPI. -ft ^ 

TABOUK » ) PLANFTARV AND LUNAR nUTPUT* SCALED AND 
TRANSLATED WITH RESPECT TO CENTFKo 

TABOUT MoJ) CONTAINS (lUTPlJT FOR 
BODY NOo Jo U oLEo J oLFo l?\ 

THE INDEX I IDENTIFIES COMPONENTS AS FOLLOWSoo 
1=X ?=Y 3*Z 

A=XDOT 5«YDUT 6 =ZDOT 

NUTt^ ) NUTATION OUTPUT 

NUT(1)=DELTA LONGITUDE 
NUTI2)=DELTA ObLiOUlTY 
NUTOIsTIME DERIVATIVE OF NUTIU 
NUTlA)aTIME DERIVATIVE OF NUTI2) 

COMMON /CETBL5/ BIVECTI 69 131 


COMMON BLOCK CETBL5 * 

BlveCTI j 1 WORKING ARRAY, CONTENTS ARF INTERPOLATED 

AND SCALED BUT NOT TRANSLATED, 1ST INDEX RUNS 
OVER X 9 Y 9 Z? XDOT 9 YDOTpZDOT AS IN TABDIIT 
BUT 2ND INDEX IS DIFFERENT AS FOLLOWS., 

BODIES I THRU 9 ARE HELIOCENTRIC, 

1 MERC 5 JUP 9 PLUTD 

2 VENUS 6 SAT 10 MOON REL TO EARTH 

3 ERTHMN 7 URANUS 11 F.RTHMN REL TO EARTH 

A MARS « NEP 12 FRTHMN REL TO MOON 

13 SEE A092+ 


THE COMMON BLOCK »CETBL9» IS FUR COMMUNICATION 
BETWEEN R0EP2 AND GETR2, 

COMMON /CETBLl/ AUpREpTPDpEMRAT 

COMMON /CETBL3/ T AB3( 8 29 ) 9 NUTAT ( 20A ) 9 CKSUM 
COMMON/CETBLA/TABOUT ( 69 12t vNUTlAl 
C0MM0N/CETBL9/jDl9TDAV» JDIF 9 1 ERR 1 
LOGICAL WFLAG 

INTEGER KREOI 121,MCENT( 11)9M1(2019JRFO( U»9 IPOSI U I 9 1 VFL ( li 1 
REAL NUTAT eSTP( 111 

DOUBLE PRECISION AU9R6pEMRAT9 TPDp TAG 39 B IVECT, TABOUT tNUT 

DOUBLE PRECISION JOl 9 TOAV 9 JOI F ,JEi) 

DOUBLE PRECISION TSECaRATp FACe U ( 2» 3 1 9 C 9 TFMP 

FAC=l/BftAOO 

DATA FAC/1.157A07A0TA07A07AD-5/ 


DATA STP/' 2,98*A.92A,5/ 

DATA KRE0/59592»6%59A939 1/ 

DATA MCENT/159l39096A15f iOfS/ 

DATA MI/II9O9IO9 HpOp 129l09 0e 12 p 0, 3p 3« 3? O 9 O9 3 9 1 3» 13 «0 «0/ 
DATA 

* IPOS /O29O929 lA69200e 25A9308»362»AU«470952A9 1/ 

♦ piVEL /A7 9 119, 1739 227»281^335 9 389, AA3 9 A97* 677 9 103/ 

DATA U/640oD0/ 


A99 CONTINUE 
JOUJED 
TDAY=TSEC*FAC 

IERRI«1ERR 
CAU GETTAP 

IFUERRl ,NE9 01 GO TO SOOO 

IFHCENT .GE. 1 «ANO. ICENT .IE. Ill GO TO 10 
IERRi«4 
GO TO 3000 
10 CONTINUE 


SET JREUUTO CONTROL INTERPOLATION 


P159DO 
2160D0 
216100 
216200 
216300 
216A00 
216S00 
216600 
216700 
216800 
2 16900 
217000 
217100 
217200 
217300 
217A00 
2176 00 
217600 
217700 
217BOO 
217900 
21POOO 
21B100 
21R200 
21R3O0 
21RA00 
21B600 
21R600 
21R700 
21RR00 
2 18900 
219000 
219100 
219200 
219300 
219A00 
219600 
219600 
219700 
219R0O 
219900 
220000 
220100 
220200 
220300 
220A00 
220600 
220600 
220700 
22 OR 00 
220900 
221000 
221100 
221200 
221300 
221400 
221600 
221600 
221700 
221ROO 
221900 
222000 
222100 
222200 
222300 
222406 


i 


A-SI 


u o o o u u 


DO 20 l»lvlO 

IF(JREg<I) ,GE. 0 oANO, IR6CH U „LE, 2) GD TO 20 
IERR1=3 
GO TO 5000 
20 JRE0( 1 )-lRE0(.I ) 

C RARYCENTFR FLAG 

JREO(3}-IREOU2» 

MA«PL=JREO( 1) 

00 2A I-2ao 

24 MAXPL=MAXO(MA«PL» JREOn > ) 

MAXEM-MAX0(IRE0(3) «IREO( 11 U 
MAXALL=HAK0rWA8PL«MAXEM) 
irfiCeWToFOsBoORaCENToEOa 11) GO TO ?« 

C CENTER IS NOT EARTH OR NOON 

C 10aM00N93=FRTHMN 

JREO( 10)sHAXEM 
JREOI 3)«MAX0( JREOO ) ,MAXEN) 

JREOI ICENT)*MAXALL 
GO TO 32 


CENTER IS EARTH OR MOON 


222500 

2226O0 

222700 

222ROO 

222900 

223000 

223100 

223200 

223300 

223400 

223500 

223600 

223700 

223200 

223900 

224000 

224100 

224200 

224300 

224400 

224500 


10=M0DN9 a^ERTHNN 

28 JREOI 1U)=MAXALL 
JRE043)=MAXPL 
32 JREOI 11)»IRE0( 13) 

LUNARsJREOnO) 

I8ARY=JRE0(3)43 


JREQI )I5 NOW SET 


SAVE*0- 
6001 CONTINUE 

DO 240 IBODYsiai 
IFUREOUBODY)) 240«240f40 
40 1F(STPUB00Y).E0,SAVE) GO TO 165 
.SAVE*STPilBOOV) 

160 TENP-JDIF/SAVE 
KK-TEMP 

U( leU=*TEMP-FLOATUK) 

IFIUt i«l m61»165»161 

161 CONTINUE 
U(2a)«l»D0-Unf 1) 

00 163 IU«U2 

unu«3)-ui !u»i )*unua> 

U4 IU92)«lunU93)-l«001/6.00 
163 UI IUv3)«fUnu»3)-4.D0)/20«D0 
165 IFUBODY-IOI 169a67t220 
167 C*RE 


172 

200 


203 


GO TO 172 
C«AU 

IGETl»aP0SUB00Y)4KK*9 
ICi-i . 

CONTINUE 

IFluntlU203»20l 
IGET2iilGETl46 
DO 204 !GET«I6ETI 
eiVECTIICI«!BOOVl 
* C«IUS2 

uu 

U«2 


INTERPOLATE IBOOV« 


4 

4 

4 

4 

4 


204 

GO TO 209 


aGET2»3 

«U4«TAB3UGET. i* 
92B*«TAB3UGET4 1»4 
»3»4 TAB3IIGET4 2BH4 
vU4«TAB3nGET4 9B4 
^2»*ITA69UGET4|0»4 
93i* TABBUGET^Uin 


l92««9*cl.O 


224600 

224700 

224800 

224900 

225000 

225100 

225200 

225300 

225400 

225500 

225600 

225700 

225800 

225900 

226000 

226100 

226200 

226300 

226400 

226500 

226600 

226700 

226800 

226900 

227000 

227100 

227200 

227300 

227400 

227500 

227600 

227700 

227800 

227900 

228000 

228100 

228200 

228300 

228400 

228500 

228600 

220700 

S2BB00 

228900 

229000 


n o o o non oooo non 


201 lC2=ICl+2 

DO 202 I*ICUIC2 

BIVECTUfl IBOOV fjsC*7AB3« IGETl ) 

202 IGETl»IGETl+3 
205 CONTINUE 

JRE0< Ib0DY)*JRE0(lB0DYJ-l 
IF( JREOnBODVH 240f240«207 
207 IGETl-g VEL(IB0DYJ^KK*9 
ICls4 
C=C/TPD 
GO TO 200 

INTERPOLATE IBODYsll NUTATION 

220 C*1«D0 

IGET1»IP0SI IB 00 VH'KK *6 
ICl^l 

222 IGET2=iGETl+3 
225 1 F(U( Iein 228 « 226«228 
228 DO 230 I GET«I GET U IGET2, 3 
NUTtlClJs 

* C*4U| 2v U*<NUTATUGET I <• 

* U« 292 )*JNUTATnGET+l> + 

Ut293H« NUTATnGET<‘2m + 

» UIl?n*(NUTATnGET-*- 6 K 

* UU,2)*«NUTATnGET+7» + 

* UUff3t* MUTAT J IGET-frH H J J 
230 ICl«ICI+l 

GO TO 232 

22A 00 227 IGET#IGETI, IGET2?3 
NUTJ lCU«C*NUTATnGET) 

227 IC1*ICI+1 
232 CONTINUE 

JREOJ IBODYJ*JREOnBC10Y)-l 
iF( jREonaoDYn 2^0*240,236 
236 C»C/TPt) 

IGET1«IVEUIB00VKKK*5 
IC1*3 
GO TO 222 
240 CONTINUE 

INTERPOLATION IS FINISHED 
RESULTS ARE IN BlVECTI 9 UNO NUT« } 
tEST MOON REOUFST 
IFILUNAR » 4020*402094010 

N0TE9*EMRAT«EARTH MASS/MCION MASS 

SET BIVECTI ,IU«ERTHMN CENTERED AT EARTH 

SET BiVECt* 9i2I«ERTHMN CENTERED AT MOON 

4010 RAT«1.00/IEMRAT->>1.D0) 

IMAX-LUNAR * 3 
. DO 4016 !«1*!MAX 
BIVECTH9iU<>RAT>t‘RIVECTn9lOI 
4016 BlVECTUtiai— EMRAT4BIVECTII, m 

4020 NFLAG-. FALSE. 

KCENTsMCENTUCENTI 

BEGIN TRANSLATION LOOP 

00 4108 lB0DY*l,i2 
IFUREOUBOOTH 4 108 *4108, 4024 
4024 |MAX«lRE0nB00Y$*3 

KASE>KCENt4KREOI IBODVI 
KI-MHKASEI 

60 TO I4O82t4O4O94O32»404A«40S2* 

* 4032«4068, 4040*404894056 « 


229100 

229200 

229300 

229400 

229500 

229600 

229700 

229R00 

229900 

230000 

230100 

230200 

230300 

230400 

230500 

230600 

230700 

230800 

230900 

231000 

231100 

231200 

231300 

231400 

231500 

231600 

231700 

231ROO 

231900 

232000 

232100 

232200 

232300 

232400 

232500 

232600 

232700 

232800 

232900 

233000 

233100 

233200 

233300 

233400 

233500 

233600 

233700 

233R00 

233900 

234000 

234100 

234200 

234300 

234400 

234500 

234600 

234700 

234800 

234900 

235000 

235100 

235200 

235300 

235400 

235500 

235600 


A*5S 


o r> o o cl n r> o no n o o non n o no 


♦ 40^2#407A^4076»<^040o40^89 

* 4084940aB94088f4064«4U80 )«KASF 

KASe«15 

402B Kl-IBODV 

KASE»l»3969ll 

4032 DO 403& l«lf IMAX 

4036 TABOUTJ I » I BODY ) «BI VECT (I p K1 ) 

GO TO 4108 

KASE«298pl4 

4040 DO 4044 I«lpIMAX 
4044 TABOUTU fllB0DY>«0oU0 
GO TO 4108 

KASe»4p9 

4048 K2«3 

GO TO 4100 

KASEaS 

4052 L2-11 

GO TO 4060 

KASEalO 

4056 L2»12 
4060 KlalBOOY 
K2«13 
GO TO 4092 

KASE«19 

4064 KlalCENT 

KASE«7 

4068 DO 4072 I*lpIHAX 

4072 TABOUTH 9 IBOOYta-Bl VECTJ I ,K1 » 

GO TO 4108 

KASE«12pl3 

4076 K2-KASE-1 
GO TO 4100 

KASE-20 

4080 KialBODV 

KASE«16 

4084 K2a|CENT 

GO TO 4100 

KASEal7»lfl 

4088 L2-ICEMT 
K2«KA5E-6 

4092 IffWFLAG) GO TO 4100 
HFLAG«»TRUE. 

BIVECTf «13)1S AN AUXILIARY VECTOR 

NEEDED NHEN KASE-5* 10« 17? 18, 

FOR KASE«05 BIVECH p13}«EARTH CENTERED AT SUN 

FOR KASE-10 BIVECTt pl3)aN0nN CENTERED AT SUN 

FOR KA$E-I7pl8 BIVECH «ri}"ERTHMN CENTERED AT iCENT 

DO 4096 I»i»IBARY 

. 4096 BIVeCmp]l3»-BIVECTn93»-°B!VECTn9L2} 

4i00 DO 4i04 falflMX 

4104 TAB0UT4 1 9lBQDY»>BlVECTU 9 KU-BI VECTH «K2» 

4i0a CONTINUE 
5000 lERRalERAl 
RETURN 
ENC 


^35700 
235800 
235900 
236000 
236100 
236200 
236300 
236400 
236500 
236600 
236700 
236800 
236900 
237000 
237100 
237200 
237300 
237400 
237500 
237600 
237700 
237800 
237900 
238000 
238100 
238200 
238300 
238400 
238500 
238600 
238700 
238800 
238900 
239000 
239100 
239200 
239300 
239400 
239500 
239600 
239700 
239800 
239900 
240000 
240100 
240200 
240300 
24 0400 
240500 
240600 
240700 
240800 
240900 
241000 
24U00 
241200 
241300 
241400 


i 


i 



i 

) 


>J 


A«ti 


SUBRGUTiwe ROTEO^TlMGpA) 
IMPLICIT REAL* 8 « A-«90-Z*»» 
OIMENSIOM A(3«3» 

T « TIMB/36525.nO 
T2 ■ T*T 
T3 ■ T2*T 


AIUU 

A(i«2» 

At2«l) 

AUvS) 

A43»U 

A(2«2» 

A(2v3( 

A(3«2) 

AI3,3I 

RETURN 

END 


1*00 - ,0002969700*12 - • 0000001 300* T3 
-,022369fl8D0*T - ,00000676D0*T2 + , 00000221 Dl)*^ 3 
-At 1,21 

-,0097171100*T ,00000207DO*T2 + ,0000009600*73 

-At 1,3J 

1.00 - ,0002497600*72 - .0000001500*73 
-.0001085900*72 - ,0000000300*73 
A(2»3( 

I, DO - .0000472100*72 + ,0000000200*73 


241500 

241600 

241700 

241800 

241900 

242000 

242100 

242200 

242300 

242400 

242500 

242600 

242700 

242800 

242900 

243000 

243100 


A-W 



SUBROUTINE RVOUT(A»B9C) 243200 
implicit REAL* 8 < A-H 9 O-Z 9 * » 243300 
DIMENSION A(3)f B(3l9C16) 243400 
C IS RAOIUSs IAT 9 tONff velocity, GAMMA, AZIMUTH 243500 
CALL SPERlAjC) 243600 
C(4)=FN0RMrBI 243700 
C(5)«90oD0-AD0T(A,B> 243B00 
VE*B(2KA(U-B( n«A(2) 243900 
VN*CUKB(3J-A(3)*D0T( A,B)/C( 1) 244000 
C( 6) =57*2957795 iOO*ARKTNS( IBO bVNjVEI 244100 
RETURN 244200 
END 244300 



SUBROUTINE SHADOWIRUV U » STEP) 

IMPLICIT REALMS! A-HpO-Z,S) * 

COMMON C( lOOOh ICI50P»Z<3)»VVI3 fpNOR 
COMMON /SHAO/SHAOK 
DIMENSION RlOpRIpVltBpfl) 

DIMENSION Kl(3t 
DIMENSION R(3pR)«V(3p8) 

DATA IPFLAG/l/» lUFL AG/ 1/9 IMF LAG/ 1 /p IMilFtr,/!/, FLAG/O* /,F LAGl/O */ 9 
lFLAG2/0<./pFLAG3/0e/ 

DATA IVFLAG/l/v IVUFLG/l/pPFLAG/Op/pPUFLAG/Oo/ 

6 FORMATUHOplSHBODY CENTER IS »A 6 ) 

IF(STEPoEO»OoDO)RETURN 

GO TO (^IpAlpSOS e400p400»4001 9 NOR 
41 CONTINUE 
DO 4 I=lfl3 

4 Xin)aftin»U-«l(!93} 

SEVA=AD0TI XlpftU 1* IJ ) 

RFS=FNOftMIRUlp3J) 

RFE»FNORM( R14 U U1 

DELSaDARSIN(65O0,D0/RFEf*57o2957795D0 
!F< SEVA-90.D0J BOO^SOOpBOI 

301 ZRS=6o965D5 
ZRE=6500,D0 

ZPSIG=DAR$INnZRS+ZREf/RFS)*57,2957795D0 

FSsSEVA^DELS+ZPSIG-lSO.DO 

IF(FS) 300*3029302 

302 FLAG*leDO 

IFURFLAGoWea) GO TO 1000 

K=1 

N=1 

IPFLAG=2 

KK=1 

WRITE(6»5) 

5 FORMAT aH0960HPR0RE ENTERED EARTH PENUMBRA. TIME AND ORBIT ARE 6 IV 
lEN BELOM) 

SHAOKal.ODO 
GO TO 20 
1000 CONTINUE 

2SIGMA*DARSINHZRS“2RE I /RFS)*57 .2957795D0 
FS»DELS-ZSIGMA+SEVA-180e00 
!FIFS)3039B04«304 
304 FLAG1«1.00 

imUFLAG.NE.U GO TO 310 

K«i 

IUFLA6*2- 

N«2 

KK«2 

WRITEI6«7I 

7 FORMATtlHO^BTHPROBE ENTERED EART^ UMBRA. TIME AND ORBIT ARE GIVEN 
IBELOWI 

SMADK«UODO 
GO TO 20 
300 FLA&«0»00 

303 FLAGlaO.OO 
310 00 30 I»l93 

30 Xlin«RUIt2»-^in93B 
SEVA*A00TUI»R1« i»2H 
IFI SEVA-90. DO I 2000« 2000*901 
90& ZRS-6.965D3 
ZRE«il73S.D0 

miG*DAIISlNUZRS4^ZRE»/RFS}937«299779900 
DELN»DAftS INI € 8 22 1 /FNOAM IR & 1 1 * 2 m«57 . 295779SD0 
FS«SEVA-FDEtN«>ZRSI6-’iaO«DO 
lflF$$ 2000* 902*902 
902 FLAG2«l*00 

tFHMFlAO.NE.U GO TO 3000 


244400 
244500 
244600 
244700 
244800 
244900 
245000 
245100 
245200 
245300 
24 54 00 
245500 
245600 
245700 
245800 
245900 
246000 
246100 
246200 
246300 
246400 
246500 
246600 
246700 
246800 
246900 
247000 
247100 
247200 
247300 
247400 
247500 
247600 
247700 
247800 
247900 
248000 
248100 
248200 
248300 
248400 
248500 
248600 
248700 
248800 
248900 
249000 
249100 
249200 
249300 
249400 
249500 
249600 
249700 
249800 
249900 
250000 
250100 
250200 
250300 
250400 
250500 
250600 
250700 
250800 
250900 





K-2 251000 

IMFLAG»2 25UOO 

N®3 251200 

KKal 251300 

WRITE(6,8) 251400 

8 FORMAT! 1HO*59HPROBE ENTERED MOUN PENIlMSRAe TIME AND ORBIT ARE GIVE 251500 

IN BELOW) 251600 

SHAOK*KODO 251700 

60 TO 20 . ' 251R00 

3000 CONTINUE 251900 

ZSIGMA=0ARSIN( ( ZRS-ZRE ) /RF S 1*57 o 295779500 252000 

FS*DELM-ZSIGMA+SeVA-lBOoOO 252100 

IF! FS12001*904j904 252200 

90^ FLAG3*K00 252300 

IF( IMUFLG.NE*U go to 2010 252400 

K*2 252500 

!MUFL6s2 252600 

N*4 252700 

KK*2 252BOO 

WRITE(6»9) 252900 

9 FORMAT ( IH0»56HFR0RE ENTERED MOON UMBRA, TIME AND ORBIT ARE GIVEN H 253000 

lELOWl 253100 

SHAOKslsODO 253200 

GO TO 20 253300 

2000 FLAG2»0,D0 253400 

2001 FLAG3«OoOO 253500 

2010 CONTINUE 253600 

IFl IPFLAGoNE*2,0R,FLAG,NE,0oD01 GO TO 11 253700 

MRITE(6»12) 253800 

12 FORMAT! 1H0*61HPR0BE HAS LEFT EARTH PFNUMHRA, TIME AND ORBIT ARE GI 253900 

IVEN BELOW) 254000 

SHA0K>0,0D0 254100 

K»1 254200 

N»5 254300 

KK=l 254400 

IPFLAG*! 254500 

GO TO 20 254600 

U IFnUFLA^..NE,2.0R,FLAGl,NE,0.n0) GO TO 14 254700 

WRITE46?. ) 254800 

13 FORMAT! Ir. -8HPR0BE HAS LEFT EARTH UMBRA, TIME AND ORBIT ARE GIVEN 254900 

I BELOW) 255000 

SHADK*0*ODO 255100 

RK*2 255200 

N*6 255300 

K»l 255400 

IUFLAG-1 255500 

GO TO 20 255600 

14 !FUMFLAG,NE,26OR*FLAG2,NE,O0OO) GO TO 15 255700 

MRITE«^*16) 255800 

16 FORMAT UH0»60HPR0eE HAS LEFT MOON PENUMBRA, TIME AND ORBIT ARE GIV 255900 

lEN BELOW) 256000 

SHADK«0«000 256100 

256200 

<<■2 256300 

N*7 256400 

IMFLAG*! 256500 

GO TO 20 256600 

15 !FnMUFL6«NE«2.0R,FLAG3«NE,0,D0) 60 TO 505 256700 

Wft|TEI6«lfi) 256800 

18 FORMAT UHO»57HPROBE HAS LEFT MOON UMBRA. TIME AND ORBIT ARE GIVEN 256900 

IBELOW) 257000 

SHAOR«0»000 257100 

*^•2 257200 

M-1 257300 

*<^■2 257400 

IMOFL6«l 257500 


A«M 


n n o 


20 CALL ITER(NORpRl«Vl«FS«K«KKpDTp A 9 V) 
TIMsC( 593 )44 X+C<5O6)+0Tf/86400o00 
CAlt TFRAC(C4592) 9 TIM 9 TT 9 TTF) 

CALL DATOUKTTtTTFpOATFpFDATEsO) 


GO TO U000931093000»2010p 11« lAp 15?505 9 405 » 502, 503 p505 ) 9 N 

400 DO 401 !Ri«3 

401 Xin )«R1( IpNORl-Rin ,3» 

SEVA^ADOTUlpRlllvNOR) } 

RFS=FNORM«RHl93H 

RFE = FN0RM4RHl»N0Rn 

DELS*DARSIN«C<N0R+20J/RFE 1*57,295779500 
IFI SEVA-90. 001500, 500*402 

402 ZRS»6.965D5 
ZRE«C(N0R+201 

ZPSIG*OARSINHZRS-»-ZREl/RFS 1*57. 295779500 

FS»SEVA+0ELS+ZPSIG“180*00 

IF4FS1500,403,403 

403 PFLAG-i.DO 

IFi IVFLAG.NE.U GO TO 405 

K»3 

N»9 

IVFLAG*2 

KK«I 

MRlTE{6,404)C(N0R-»-601 

404 FORMAT* 1H0,26HPR0BE ENTERED PENUMBRA OF ,A6,30HTIME AND ORBIT ARE 
IGIVEN BELOW) 

SHAOK-l.ODO 
GO TO 20 

405 CONTINUE 

ZS I6MA»0ARS IN n ZRS-ZRE 1/RFS 1*57. 295779500 

FS<=DELS-ZSIGMA-^SEVA-180.D0 

|F«F51501,404,406 

406 PUFLAG>1.00 
IFUVUFLG.NE.il GO TO 502 
Ks 3 

IVUFLG-2 

N»iO 

KK*2 

WRITE46,4071 C4N0R+601 

407 F0RMAmH0,23HPR0BE ENTERED UMBRA OF ,A6,30HT1ME AND ORBIT ARE GIV 
lEN BELDMS' 

5HAOK»i.ODO 
GO TO 20 

500 PFLAG-O.DO 

501 PUFLAGbO.OO 

502 IFS IVFLAG.NE.2.0R.PFLAG.NE. 0.001 GO TO 503 
WRITES6,5041 CIN0R*601 

504 FORMAT UH0927HPR0BE MAS LEFT PENUMBRA OF ,A6»30HTIME AND ORBIT ARE 
I GIVEN BELOW! 

SHADR>0.000 

K«3 

N-ll 

KK«i 

lVFLAG-1 

GO TO 20 

503 1F« IVUFi6.NE«2.0R*PUFLA6.NE.O.OO» GO TO 505 
WR|TE«G»506! C«NQIt460» 

506 FORMAT II HO »24MPR0BE HAS LEFT UMBRA OF *A6«30MfriNE AMO ORBIT ARE G1 
IVEM BiLOIII 
SHAORvO.OOO 
KK»2 

M»12 

Rtt3 


257600 

257700 

257BOO 

257900 

25R000 

25R100 

25R200 

25R300 

25R400 

25R500 

25R600 

25R700 

25RR00 

25R900 

259000 

259100 

259200 

259300 

259400 

259500 

259600 

259700 

259R00 

269900 

260000 

260100 

260200 

260300 

260400 

260500 

260600 

260700 

260600 

260900 

261000 

261100 

261200 

261300 

261400 

261500 

261600 

261700 

261800 

261900 

262000 

262100 

262200 

262300 

262400 

262500 

262600 

262700 

262800 

262900 

263000 

263100 

263200 

263300 

263400 

263500 

263600 

263TO0 

263800 

263900 

264000 

264100 


I 

\* 




I 




IVUPLG-I 26420a 

GO TO 20 264300 

505 CONTINUE 264400 

RETURN 264500 

END 264600 


SIJrtROUTIMe SHIFTPfNrjRpRjVjTWsTFvRR^VV) 

IMPLICIT RFAL»Xfl ( fi-M,n-Z9*l • 

DIMENSION C( 1000) p I C< 50 ) 9 K I 3 ) * Vf 3 ) »«R ( 3? « ) 9 VV ( A) 9 Pn t ?2 ) p Vt ( 22 I * 
IRK 392),Vl(392)9N(2)9SPAR(20p2) 9 TK 2p2) 

DIMENSION XlSIpWn) 

COMMON C 9 IC 

F.OU I VALENCE (C(AA3 ) 9 R 1 ( I 9 1> ) 9 ( C( 6B 1 ) 9 V 1< 1 9 1 > ) 9! CI695) 9 TI ( 1 9 !) ), 
l(IC( 1 ) 9 ItiJ) 9 nC(2) 9 IB 2 ) 9 ( IC( 12 I pNt 1 Mp ( C( 71), SPMSn) 9 ( C13 ) 9 OS) 
2,<C(801)9P0( IMp (Ci«23)9Vem )v(C{ 5611 9 SPAR ( 1,1) ) 

39iU6991 pSWNET) 

IC(2)=0,n0 
IPfSWNbTKOOp IOO 9 18 

100 nis«i.oio 
no 1 1 = 1,3 
xm=RU) 

1 w( n=v( I ) 

NM1«N0R-1 
IF<N0R-6)4f 3,2 

2 KKK=-l 
NM1 = 2 
GO TO 7 

3 NMl = 6 
^ CONTINUE 
Kt<K = 0 

I F I ItJ I “6 ) 5 9 5 96 

5 IP{ 182-6)9,9,6 

6 KKK =1 

7 CONTINUE 
DO 8 1=1,2 
K= ! 

TT=U0INT(TW-T1( 1,1 M+TF-TK 2, I ) )*SPMSO 

CALL STEPDINd 1 , TT ,T A, R U 1 9 1 ) 9 V 1 f 1 , 1) , OS ,RK ( 1 9 K ) 9 VVI 1 , K) , 1 , 
ISPARdpin 

8 CONTINUE 

9 CONTINUE 

CALL INTRHTHpTFpNMlpPn, 1 ,VE,D 1 S) 

IF{KKK)1D912,12 

10 NOB=NUR 

DO 11 1=1,3 
X( n=X(I )+RR( I pNOR) 

11 M( n=Nt n^vvn pNOR) 

12 CONTINUE 
00 13 J»l,5 
on 13 1=1,3 
L«3*J-3+I 

RR« IpJl»X( n-PO(U 

13 vvn9ji»wm-vEtu 
on 16 1=1,3 

RR n »6 i*x m -POU+ le » 

16 VVI l,6)=wm-VEH-»'lfU 
lF|KKKJl5,l7pl5 

15 DO 16 l»lp2 
K=I-f6 

DO 16 J=l,3 

RRf J,R)«RRIJ»31-RRI J 9 IU 

16 VVt J,KI»VVU93)-VVt J,K} 

17 CONTINUE 
GO TO 35 

iR TT=UTV«-TUl»£H->-TF-Tll2t IH^SPHSO 

CALL STEPOlM,TT,TA,Rm,U,ViU, lUCISURRIi, U,WU,D*1, 
isPAftfian 

TT« n T»#-TU 1 *2 1 S +TF-T 1« 2, 2 n*SPNSO 
1F« ICa I *IC«2 1-3 619,19 9 25 
19 M»2 

CALL STEH)«N,TT,TA,RU 1*2» , ViU , 21 tU U ,RR I 1,2 B vVVIi ,2 » «1 , 


26670D 
26^800 
266900 
265000 
265100 
2652O0 
2 65300 
265600 
265500 
265600 
265700 
265800 
265900 
266000 
266100 
266200 
266300 
266600 
266500 
266600 
266700 
266800 
266900 
267000 
267100 
267200 
267300 
267600 
267500 
267600 
267700 
267800 
267900 
268000 
268100 
268200 
268300 
268600 
268500 
268600 
268700 
2688 00 
268900 
269000 
269100 
269 2 00 
269300 
269600 
269500 
269600 
269700 
269800 
269900 
270000 
270100 
270200 
270300 
270600 
270500 
'270600 
270700 
270600 
270900 
271000 
27UOO 
271200 




lSPAR(l,2n 
fjn zk 1 = 1,3 
IF(N(JR-;^)20,21 ,22 

20 RR( 1 ,3> = -RR < I , I ) 

VV( I ,3J=-VV( I ,1) 

RR< I #n=0oD0 
vvn 9 1 )=n,f)0 

GO TO 23 

21 RRn,3»=:“R(<( I,U-RR(I,2) 
vv( ip3)=-vvn»ii"vvn,2» 

RR(1 ,2»=0o00 

VV( 1 ,2)=0.D0 
GO TO 23 

22 RR( I,2) = RR( l,2)+RR( 1,1) 

VVn,2)*VV( I,2J^VV< I ,1) 

RR( I ,3)^0a[)0 

vvn ,3)sOoDO 

23 DO 24 J=l,3 
RRn,j)=R(n-RRn p j> 

VV( !,J)«V< I )-VV(I, J) 

24 cnwTjwue 
GO TO 35» 

25 M=2 

IF1 ICa)-3)27,27,26 

26 KalCtU 
GO TO 28 

27 K=IC<2) 

28 CONTiNUe 

CALL STEPDlMBTT,TA,Rm,2) ,VH 1, 2 ) ,Cf 3 ) , RR U , tU , VV ( I , K) , 1 , 
ISPARd, 2n 
DO 34 1*1,3 
IF(N0R-3)29,3l,32 

29 RRU ,3)*-RRf I ,1) 
vvn ,3) = -VVU, 1 ) 

IF(K,E0,3) GO TO 30 

RR4 I ,K) = RR( E ,K)-RR( I , 1 I 
VVn,KJ*VVn,K)-WJ 1,1) 

30 RR< I ,1)»0,D0 
VVU,U»OoOO 
GO TO 33 

31 RRUsSl^OoDO 
VV4 1 ,3)*0.00 
GO TO 33 

32 RRII,U*-RRU,K)+RRn, IJ 
VV«I,U»“VV« S»K)+VVn,l) 

RRS S ,3)»-RRn ,K) 

vvi i,3)*-vvn,K) 

RRI !,KI»0.00 

VVJ ! ,K)«0»D0 

33 RR4 1 )-RRn« n 

vv(uu*vf n-vvn» u 

RR« I ,3)«RI I )*RRU ,3) 
vvn*3)«vn)-vvn,3i 

IPSK«E0.1.0R«K.E0*3» GO TO 34 
RRH )<-RRU ,K) 

vvn«)o>vm-wu«K) 

34 C0(4TJNU£ 

39 CONTINUE 

RETURN 

END 


271300 
271400 
271300 
271600 
271700 
271 POO 
271900 
272000 
272100 
272200 
272300 
272400 
272500 
272600 
272700 
27 2 POO 
272900 
273000 
273100 
27 3 2 00 
273300 
273400 
273500 
273600 
273700 
273800 
273900 
27^000 
27^100 
27 A 200 
274300 
274400 
274500 
274600 

274700 
274 POO 
274O00 
275 (»00 
275100 
275200 
275300 
2754f*0 
275500 
2756f»0 
2 75700 
275RO0 
^75900 
276000 
276100 
276200 
2 76300 
276400 
276500 
276600 
276700 
2768JJO 
276900 
277000 
277100 
277200 




SUBRnUTlNC SOLAR (RtSACC) 

277300 

IMPLICIT RFAL ♦ fl lA-Hvd-ZI 

277400 

COMMON C(IUOO) 

277500 

DIMENSION KI3USACU3) 9RX(3) 

277400 

P*C(774> 

277700 

A=C1770I 

277800 

eM*CI77U 

277900 

GAMMA*C« 772) 

27 B 000 

SIGMA=C(773) 

278100 

Cl = P'«‘A>»‘noO + 2 oO/3oO«GAMMA + l»0/3o0*GAMMA=«‘SiaMA)/EM 

278200 

DO 11 J>1,3 

278300 

RXU)»RU) / le^9‘»D8 

278400 

CONTINUE 

278500 

RAl)s(RXm*#2 + R«(2)«<' 2 +RX(3 )*«2)«*0*8 

278400 

C2-Cl/RAIJ*«3 

278700 

00 12 K«U3 

278800 

SACC(K)=C2*RX(K)«1oOO-5 

278900 

C 0 N T ! N U E 

27 9 f >00 

RETURN 

2791(»() 

END 

279200 


AM 



SUBROIjTlWfr SPERIX^Y) 

IMRLrCIT REAL^RI A-HflO-Z»S) 

C CnWVERTS CARTESIAN TO SPHERICAL 

oiM^hisinAj xf3)»vm 

C OlMENSlOH X(3»fV(3) 

RTr)=«»»72957795D+02 
A*KU J*Xm+Xl2)*X<2J 
Y( n*OSORT(A+KI3I*X(3J ) 
Y<2 )bDATAN(X( 3J /DSORTl An*KTD 
YJ 3}«ARKTHSUftOs Xmj^XfZM^RTO 
RETURN 
END 


279300 

279A00 

RETURNS In OEO 279500 

27 9ft no 

279700 

279ROO 

279900 

2BOOOO 

2B0100 

2B0200 

280300 

280400 





onooon o o o n n e*> n o on n 


SUBROUTINE STEPOfN « T « TA,R » V» U«RR 9 V SB) 
IMPLICIT REALMS (A-H90-Z 9 $) . 

CB IS A STORAGE VECTOR WHOSE ELEMENTS 


m INITIAL TRUE ANOMALV ( 11 ) 

(2) SEMI-LATUS RECTUM ( 12 ) 

(3) ECCENTRICITY (13) 

(4) RAOIUS-SOUAREO (14) 

(5) RADIUS MAGNITUr>E (15) 

( 6 ) SEMI-MAJOR AXIS (A) (16) 

(7) ABSOLUTE VALUE OF ( A ) (17) 

( MEAN MOTION ( 16) 

(9) TIME COEFFICIENT R/ABA 


( 19) 

(lO)TIME COEEFICIENT RV/S0RT(UA) (20) 

DIMENSION R(3)9V(3),HVEC( 3) «KR( 3) 9 VV ( 3) 9 E( 4 ) 9 CB( 20) 
DIMENSION SB(20) 

dimension R(3)9V(3)9HVEC(3)9RK( 3) 9 VV( 3) 9 E(4 ) 9 CB( 20 ) 
DIMENSION SB(20) 


ARE 

ANGULAR MOMENTUM SUUARED H? 
ANGULAR MOMENTUM MAGNITUDE H 
MU TIMES RADIUS MAGNITUDE UR 
R DOT V 
H?-U» 

h*rdutv 

U/(R*H2) 

RV/H 

TIME FROM PFRIAPSIS ( -o- OR -) 
SWITCH I ELLIP 9 2 HVPER 


N«N 

GO TO (l97)9N 

1 CONTINUE 
PI».314159265D+01 
TPI«*62«31fl5300+01 
CB(4)>cD0T(R9R! 

CB( 5)*OSORT(C6(4) ) 

V2*D0TIV9VS 
CBf i4)«D0T(R9V) 

VMnOSQRT ( VP ) 

CB(6)«CB(5)/I2.0D0-CB(5)>»V2/U) 
CB(7)>DABS(CB(6) } 

CB(8»«DS0RT(U/CB(7H/CB(7) 

CB(9|«CB(5)/CB|7) 

CB( 14}/0SDRT(U*CH(7H 

CALL CROSSIRfVfHVEC) 

CBI lU-OOTtHViCfHVECI 

CB(i2»>DS0RT(Cfimn 

CBf 13)«U*C6I5I 

CB(2B«CBmi/U 

EC2- 1 •DOO-CB 1 2 ) /cm 6 ) 

IR€C2»6E.I-l«0-6).AN0.EC2>LE.D«n0) GO TO 200 
CBf 3I«DSQRT(EC2) 

GO TO 201 

200 CBf3)-0*DO 

201 CONTINUE 
CBU5B«CBUU-CBf 131 
CBa6»-CBn2I«CBf 14) 

CBf 17B«U/ICBi5)*CmiU ) 

CBf i8|BCBa4)/CSf 121 

CBf U«AftKTNSf iSOvCBf 15 ) tCBf IBH 

TA2*CBf U 

lEfCBf6H3«392 

2 CONTINUE 

CBf 20) « I, 000 
EI»CBf B)4T 
GO TO 4 

3 CONTINUE 

CBf 20 I »2 *000 
E&«O.OfX» 

4 CONTINUE 

CALL TCCMICfU«C6f3B«CBf6)«CBf 2l9TA2«T(^p9EAC» 

N«N 

GO TO f6»5l9N 

5 CONTINUE 


PH0500 
280600 
280700 
280RO0 
280900 
281000 
281100 
28 1 200 
281300 
281400 
281500 
281600 
281700 
281800 
281900 
282000 
282100 
282200 
282300 
282400 
282500 
282600 
202700 
282800 
282900 
283000 
283100 
283200 
283300 
283400 
283500 
283600 
283700 
283800 
283900 
284000 
284100 
284200 
284300 
284400 
284500 
284600 
284700 
284800 
284900 
285 000 
285100 
285200 
265300 
285400 
285500 
285600 
285700 
285800 
265900 
286000 
286100 
286200 
286300 
286400 
2 86500 
286600 
266700 
286800 
286900 
207000 


A-M 



T»TFP-C6< 19)+T4P>*'P0 


2«7100 


GO TO U 


287200 

6 

CONTINUE 


287300 


CH( 19)*TFP 


287400 


N-2 


287500 


DO 60 l»U2n 


207600 

60 

SBf I > 


287700 


GO TO 700 


207800 

7 

CONTINUE 


287900 


00 70 l*l?20 


288000 

70 

connstif I ) 


288100 

700 

CONTINUE 


288200 


P0«TPI/CBt8j 


288300 


IF4M)8«a,ll 


288400 

H 

CONTINUE 


288500 


TA IS THE INCRENENTAL TRUE 

ANOHALY 

208600 


CT«DCOS(TA» 


288700 


ST«0SIN{ TAI 


288800 


EE*CT-ST*CBnflI 


288900 


GE*ST*COUI/CB« 12) 


289000 


R20ftl»CBmj/(CBt m+CBI 15KCT-CB4 16>#ST» 

289100 


EP«R20R1*EF 


289200 


GE*R20Ri*GE 


289300 


EFO*CB417J*«Cm lAI-CBI lA |»CT-CB 1 12J »ST S 

289400 


GEDa(Ctia5H-CH« 13J*CTI/CB4 lU 

289500 


TA2«CBU ntA 


289600 


TAP»I0INT«TA2/TPn 


289700 


TA2*TA2-TP|*TAP 


289800 


lFnOlNT4TA2/PI n 


289900 

9 

TAP*TAP-1.0D0 


290000 


GO TO 4 


290100 

io 

TAP-TAP+KOOO 


290200 


GO TO 4 


290300 

u 

CONTINUE 


290400 


K«CEH20I 


290500 


ENOT*Ui8 8f*T 


290600 


CAU. G0TQRIK«EM0T«Cb49nF* 

EU 

290700 


EE»“F«2J/CBfl9m.OOO 


290800 


GE«-Fm/CB«e»+T 


290900 


R0A«FI2I'KB«9I«F44I-«-CB9 !0I«F43I 

291000 


eFO«-Ctsm*CB««»*F43l/4CBl 5J*R0AS 

291100 


GEO»-F i 2 J /ROA+ 1 • ODO 


291200 


TAsARKTNSUBOvCRUI^EF-fCBI 14I*GE«CBI i2I*GEl 

291300 


IF«K-miO, 110*12 


291400 

no 

CONTINUE 


291500 


TAF»iOENm/POI 


;^91600 


|FITA«mH*112»U2 


291700 

in 

TAP«TAF40S IGN U . 000* T 1 


29 1800 

112 

TA«1A4TAP«TP| 


291900 

12 

CONTINUE 


292000 


00 13 l»I»3 


292100 


n«EF*ft 1 1 1 


292200 


vv u i -EFo*R m -»GEO«v 1 n 


292300 

13 

CONTINUE 


2924 00 


RETURN 


292500 


END 


292600 


SUenuUTINE TCONIC«U»ECsA,SU9TA2oT,FACJ ?92700 

IMPLICIT REAL^fif A“«90“E9* » r92ftO0 

TANGt OOOOFL ) *DS IN (OOOOFL I /DCOSIOOOOFL I 292900 

AB w DA0SIAJ 293000 

FAC»AB*DSORT( AB/UI 293100 

ecA*aoODo-ECs/a.ooo-«-ec» 293200 

ABF * DSORTC OABStFCA) ) 293300 

TH6*TANGU5nO'9'TA2) 293400 

IrlABE-.OOOOSOOm? 11912 293500 

12 CONTINUE 293600 

ECA*2.000*DATAN( ARE<‘THEJ 293700 

IFtA)149H9l3 293ROO 

13 T«FAC*JECA-FC*nSIN(ECAI} 293900 

GO TO 16 294000 

14 ANG««7B539A 16300 + ECA/2.D0 294100 

T»FAC4«EC*TANG!ECA>“DLnGtTANG4 ANGH > 294200 

GO TO 16 294300 

U FAC*DS0RT(SLP**3/«n*2#000/n U0n0+ECI**2J 294400 

ECl«eCA*THE**2 294500 

T*FAC*4THE-»-THE4434n l,0n0-2-0D0*FCA» /3.0n0-« ?*Or>0-3»000*ECA )*ECl/5 294600 

1 ,000+ f 3*000-4. 0004ECAJ4EC 1*42/7. 000-(4*ODO-5»OnO*FCAf*EC 14*3/9*000 294700 

2H 294 POO 

16 CONTINUE 294900 

RETURN 295000 

END 295100 


A-i? 


SlJhUfJUTING TERM(K 9 Y,DflF 1 
1‘^PLICIT REAL*R(d-H»0-Zp$) 

COMWDN U 1000) »IC(50) 9 Z (3) fVV( 3) 9N0R 
CDMMQN/APPER/ FPftOLD jOLO^XMAX 
Jl^^ENSIUM RH3«R>9Y11398) 
orwCNSlO^ Y(20)9D(20) 

CUMWON/STEPS/ STEPsXOLOn 

EOIJIVALEWCB < C « 40U p R U le 1 ) ) 9 (C( A26 ) 9 v/l ( 1 , 1 ) ) 
O&Ta XOLD/0./ 

STEP = X-XOLO 
XDLfJsX 
XriLOf) = XOUD 
V( 7)-X+Cf 506) 

C(4A1) s C(^ 6 U + loD0 
P=C(b94)-<C(506)+X) 

IF(C( 776) »E0,0*D0) GO TO 5 

FPR=C« 776)*«90*D0'-AD0T(RU I 9 NOK tpVll UWOR) ) ) 

IF (FPR^GT. OoDO 9 AN 0 * X oNEq O-iOO) nLD=FPK 
IF(FPROLDoEO.O*OOoANDoFPR*GT.OonO) GO TO 5 
IF (X ,E0. 0.00 .OR. OLO .EQ. 0.00) GO TO 5 
FPRf)LD=l.DO 
F = FPR 

5 CONTINUE 

IF(C(1000))4000»10910 

10 IFIOABSIX+CISOfiH.LE.XMAX) GU TO AOOO 
XMAX=OAtiS(X+C( S06) ) 

IF( XNAX-CI 1000) *86400. 00 >-*90009 1 19 11 

11 CONTBWUE 

CALL SHAD0WIR19V19XQLD9STEP) 

4()(jn CONTINUE 

IFIPNORMIRUlpNOR) ).GT*C(N0R+20) ) GO TO 2 
WRITE(6>3) CJNOR+60J 
3 FnRMATIlHOpHHPROBE HAS IMPACTED pA6) 

F^O.OO 

Cl 122)*10.00 
CI462)x0.00 
RETURN 
2 CONTINUE 
I0Ra3 

IFIFNflRMlRll I 9 1) }.LT,C( U ) ) lOKxI 

IF(FNORM(RHl92S).LToCt 12) ) I0R*2 

IFIFnoKMIRU l9^) I.LT.CI 14 ) ) I0R«4 

iFIFNDRHJKltlfS ) ).LT,C« 15) > I0Rx5 

lF(FNnRMfRm*6) ).LT*CI 16 ) ) I0R*6 

IF( lORoEO.NOR ) GO TO 1 

-MORxinR 

FxO.OO 

CU22)o2.DO 

RETURN 

1 IFrC(i22).EO.O.D0.AN0.F.GT*O.OO) F-OeOO 
RETURN 
FND 


295200 

295300 

295400 

295500 

295600 

295700 

295ROO 

295900 

296000 

296100 

296200 

296300 

296400 

296500 

296600 

296700 

296800 

296900 

297000 

297100 

297200 

297300 

297400 

297500 

297600 

297700 

297800 

297900 

298000 

29R100 

298200 

298300 

29F4O0 

298500 

298600 

298700 

298800 

298900 

299000 

299100 

299200 

299300 

299400 

299500 

299600 

299700 

299800 

299900 

300000 

300100 

300200 

300300 


A-M 



SU6R0UT!ME TFR AC 1 T 1 » T2 » T W »TF ) 

300400 


IMPLICIT REAL#B( A-HyO-ZfS). 

3 OOS 00 


TW=*T1 

30t)ft00 


TF«T2 

SL'07 00 

3 

CONTmUE 

3008CU 


lFfTFJU292 

3 00*? 00 

1 

CONTINUE 

301000 


TW«TW-KODO 

301100 


TF=TF+l,ODO 

301 200 


GO TO 3 

301300 

2 

i ONTINUE 

301400 


TDUM=I0INT(TF» 

301300 


TW-TW TOUM 

301A00 


TF=TF-TOUM 

301700 


RETURN 

301R00 


END 1 

301«?00 





o o o o o n o 


^.iBRUUTlNg TIMFC(T1,T29T3,T^» 

IMPLICIT REAL^flf A-H90-Z9S) 

n IS AN INPUT IN FORM ( YEARS FROM 1900) (MONTH OF YEAR) 
(DAY OF MONTH 

WRITTEN AS 6501*12 FOR JAN„ 12, 1965 

T2 IS AN INPUT IN FORM (HOUR OF HAYMMINUTE OF HOUR) 

(SFCONO OF MINUTE) 

WRITTEN AS 1301*3032 FOR IPM, 1 MINUTF 9 30,32 SFCONO 

T3 IS WHOLE DAYS FROM 1950 OUTPUT 

T4=0,ODO 

TA=IDINT(T1) 

TAl=TA/l00o0D0+,01D0 
TYU IDINT(TAl) 

TY= TY1-60,0D0 
IY=TY 

IF UY) 15,15*16 

15 CONTINUE 

WRITE ( 6,17) 

17 FORMAT (37H EXIT FROM TIMEC, DATE 1960 DR MEFORE) 

CALL EXIT 

16 CONTINUE 
ND=3652 
K“ 1 

UO 1 1=1, lY 
KK=K-I 

IF (KK) 2,2,3 

2 CONTINUE 
NO=NO + 366 
K=K+4 

GO TO 1 

3 CONTINUE 
NO=NO +365 

1 CONTINUE 

TV2sTA“TVIi«100o0DO+. iOO 
!V2=IDINT(TY2 ) 

IF (IY2) 10,10,11 
U CONTINUE 

IF a2-lY2) 10,12,12 
10 CONTINUE 

WRITE ( 6,13) 

13 FORMATI37H INPUT TIME IS WRUNG, FXIT FROM TIMEC) 

CALL EXIT 
12 CONTINUE 
00 4 1*1, IV2 
1*1 

C JAN FEB MAR APRIL MAY JUNE JULY AUG SEPT OCT NOV 

GO TO ( 4 ,6 ,7,6 , A ,6 , 6 ,6 , B ,6 ,B ) , I 

6 CONTINUE 
NO*ND + 31 
GO TO 4 

7 CONTINUE 
!F«KK-l»9sl4,9 

14 CONTINUE 
ND»NO 4>29 
GO TO 4 
9 CONTINUE 
ND»NO +2B 
GO TO 4 
B CONTINUE 
NO*NO 4^30 

4 CONTINUE 

T8l*m-TA»*lOOeOOO+eiOO • 

N0«N04 JOINT nan 

ND*NO-l 

T3«N0 

TA«I0INTIT2B 


302000 
302ino 
302200 
302300 
302400 
302500 
302600 
302700 
302R00 
3 02900 
303000 
303100 
303200 
303300 
303400 
303500 
303600 
303700 
303ft00 
303900 
304000 
304100 
304200 
3 04300 
3044 00 
304500 
304600 
304700 
304RO0 
30^900 
305000 
305 100 
3052('(' 
305300 
305400 
3 0*! 5 00 
305600 
30*^700 
30*i R 00 
305<»00 
306000 
306100 
30620t» 
3 06'-i00 
306400 
31 6500 
306600 
306700 
306R<'0 
3 06900 
307000 
307100 
307200 
307300 
307400 
307500 
307600 
307700 
307A00 
307900 
3UAOOO 
30ASO0 
30A21H) 
30A300 
30P4«N) 
30A500 




T4-< TiJ-TA) /R64o0n0 
T8*ta/ 1005000+, inn 
TBlsIDlNTI TB) 

TB2=TA-TBl<‘100eOn0+, ino 
TB3sIDINT(TB2» 

T-«f=T4+T81/24o00O+TB3/2A,0D0/60,0n0 

RETURN 

END 


^lOP^nu 
30R700 
30Hfln0 
3O«PO0 
309000 
309100 
309200 
3 093 00 


▲-n 



SUBROOTIftte TRAJJTlOflTlFsRsVjNF > 
implicit REAL*8{ A-HsO-Z,* > 

COMMON CU000>9lCf50KZ{3)9VV(3.UN0R 
COMMON/TIM6/TW jTF 

DIMENSION RI(3,(n,VU3»e)«R(3)« V(3) 9Y(20)«D(20» 

EQUIVALENCE (C (401 1,RU U I H p (CI425 WVlf 1? 1 H 

J*2 

N»6 

L^O 

C(46U«0«DO 
C(506)«Q.D0 
C( 590J*T1D 
C(59U»TIF 

C( 5941 «{C( 5901 -C( 592 H-C( 591) -C( 593 H4R6400. DO 
IF(C(594)«E0.0.00) RETURN 
Y(9)«-1,D0 
M*C ( 763 ) 

G«DSIGNU20.00»C(594)) 

Y(8)*C(764) 

2 DO 1 1-1«6 
I Y(n*O.DO 

XoO.OO 

CALL FN0L2<J»N«G,L»M9NEvXfV909DER!V» TERM, OUT) 
lF(Ca22)«EO«L.DO) GO TO 3 
IF(C(122).E0«2.00) GO TO 6 
1F(C(122) «EO*lO«OOt RETURN 
IF(C(594)•(C(506)'^XI.LE.0«00) GO TO 3 
4 CONTINUE 

C(506)«C(506)-«-X 
Cri22)«leOO 
DO 5 I«i,3 
C( !*499)«Rin,N0R) 

5 C4 l-»-502)>VUI vNOR) 

GO TO 2 

6 CONTINUE 
ca 22 )» 1.00 

WRITE(6,7) C(NOR4.60) 

7 FORMAT 120H0e0DV CENTER 1$ NOW A6) 

DO a 1-1,3 

C(H499)-RU IfNOR) 

S Cn4502)«VlU «NOR) 

C( 506)-C(506>-«'X 
GO TO 2 

3 CONTINUE 

DO 13 1-1,3 
R(‘n»R&n,NOR) 

13 vm-vu I, NOR) 

Ca22)-1.D0 

CaS92)-TW 

CI593S-TF 

RETURN 

END 


309400 

309500 

309600 

309700 

309800 

309900 

310000 

310100 

310200 

310300 

310400 

310500 

310600 

310700 

310800 

310900 

311000 

311100 

311200 

311300 

311400 

311500 

311600 

311700 

311800 

311900 

312000 

312100 

312200 

312300 

312400 

312500 

312600 

312700 

312800 

312900 

313000 

313100 

313200 

313300 

313400 

313500 

313600 

313700 

313800 

313900 

314000 

314100 

314200 

314300 

314400 

314500 








Case 1 i@ an example of the input and ou^ut formats for printout at apoapsis, 


INPUT 


C array 


Variable 

Semimajor axis of 535 GO km« 
Eccentricity is .86915888 
True anomaly is 180 deg. 

Ei^t ascension of node is 45 deg. 
Inclination to reference plane is 30 deg. 
Argument of periapsis is 45 deg. 

Epoch date is March 13s 1971 
Epoch time is O*' 0“ 0 * 

Input reference coordinate system is 
Earth mean equator and equinox of 1950, 


po 

time of 4 
1 . 
po 


in<ii 


I ) mTm 




first bloek of output shows that shadow times have not been requested and 
central body obiateness of the Earth, Moon or Mars, whichever is the centra! 
, is requested. However in this case, no oblateness is calculated since Venus 
e central boc^. 

next block of output is the epoch time and the input transformed into the 
A Earth equator and equinox of 1950, January 1 at 0 ^ , which is the working 
dinate system of TRIP. 

comes the epoch input transformed to those output options specified. In 
case, the initial position of the satellite was at apoapsls. This is followed 
:milar outout at the reouested times. The outout definitions are as follower 




THET 

APOG 

MAN 

P^Rtf 

TPER 


True anomaly in degrees 

Apoapsls in km 

Mean anomaly In degrees 

Period in hours 

Time to periapsis in days 


In the final block is the number of Integration st^ taken, ^ce printout at 

uires restartin 

grator, this is the number of steps taken from the previous output time. 





('W'-KLjiy lhu>|lT 


750 

41.53 sounooD 


751 

Ha 

Oil 

752 

r,l4i;ii(ii)(>liO (ii 

7*1^ 

ll.4R(>(Mlll|l|lO 

(It* 

75 * 

U.3UOOUOOOI) 


755 


(V 

754 

I'.TH'^riOOl^ U4 

TS7 

D.l) 


75*1 

U. 710317000 

Ikfr 

740 

M.4I 


741 

I'.AoiuxiDiion til 


II. HllltllllMlIlH 


744 

U.U 


744 

. f>a4rirtMkiir;ni) 

01 

T75 

1 , 10OMtO(H)l< 0 ) 

77h 

II, limill'IK'H') 

lil 

777 

4J,1UI110H.’'I0 

r\ 

0 

- 


0 

0.0 

1‘ 

U.O 



tS MIT Ht^lNIFSTPII 

HAMhIMjCS UF HAkTH I INCLUDliU 

HAMr.|l«ICS OF hlHW IS |MrLH»Fl) 

MADMIwICS OF MARS IS lMr.Lm>FII 

SHLAM KADIaT ll>M PHFSSIWt I& NUT iNCLtitIFIJ 
INPUT IN NkAN to AMO FO UF 19S0 

■>*ANCH 13f i9Tl« OHMS* OMH* U> OSFf. J|I4,|AM ijaT*. 4»S|>Oli(U.flli 

)l-U(AA*)»7ll>in fW V-0.9SNOli!fcNO (IF 7-IU OR »S 4itW»fl?ni7TU tKl l>V 0,*TAAT»l •••'“IM 'l/-l!.^ tn 

H U< KHMrUUUCIt) Ofr OFC^OaTOTOAeiAn O? R A-n«9A10AS9nfl U7 W MthA|AAHl70 On PTH-il.NPOnPAAIirVOA A/ (i. 1 1 ^ 7(1 (fcFti ii^ 


(4ARCH I3f 197li OHftSt (VttNf ()• USFC JIiIIa>m l>aTfc /AA ) ii^^ • S(ii>('immiii 

vemtis CFNTEhtO CEMTRAL Winv HfeAN fcA»Th FIIKATOk ANO FCHlIMtllt IIP 191>fi 

V-U*hA 9 HTlAln fl* V- 0 > 91 30 lift 90 01 i-n.lblSISAAO 01 . OB ll.AOft?fU?TO 00 IlV O.AlAATIl'iri-Ul Hi’' 0 .>:, 1 lH 7 li|i 1 ii Oil 

R (jiinouuuuon no opc-o>2o7n4Bitri o? ra-o.raioairro o? v o.aiihahipo on PTH-o.^jimpAAort-oA */ (i.ii/’/MT'-Fp 

SNA 0.93300U00D 05 PCC 0«8fc9l5B<Htn no |Mf. 0.100000010 O? LAM O.AICKIIHKtlU O? APP Il.AIOl.OlHi 1 n U7 xlia ii. #,yNNy«A»ti oA 

C3-U>fv07{H579n 01 THPT'KIt lAOlKHHHin ftj APno O.IOOOOOOOfl UA »AW O, I PlriWUlOlO 01 PPmO o.1TFOlPV.Mn u? TuF«-o.Vhnsami?h ho 

VtNOS CFMTtOtn TRUE ECLIPTIC AMO EOUlNOA OF OATfe 

X-U.AlTBA|77n 04 V-U>9qf><iq0{l30 01 7 O.AtoHl (iPOVI 04 OK U.00R443410 OO 0V-O.4n41?14V>-ul I1/-0 . ,!chhI S 7*.li oo 

B 0, lOlKtUODUn Ch OEC 0.2t"33517n 01 R A-0*O354673T0 (l? V 0.h1|R4Hl1U 00 PTN-n. 1709 u 410-U4 m 0.ilvm<‘l'‘ Hi 

SNA U,5310UU(Km 05 FCC a.BA9l5ltHbn on IMC 0.707197440 U? LAM 0.915415950 07 APF 0.157401050 01 XLA ii. ^inmOoiK'*' 04 

C 3 -U, 4 nT 04579 n fll TPFT- 0 , IRnOOUUUO 01 APOr, 0 . JuniHJOOOn 04 PAM O.lflOnnuOlil (I 1 PFhM 0 , 17 A 9 p 141 o U/ 1PFk-.H.7iS9549(>1Ii I'll 

SUM CEMTEKEO TRUE ECLIPTIC AMO EOUTMnX Of OATb 

4-0.45021)22411 OR 7-0,984204570 0« f. 0.122527270 U7 04 O. 32212205|> 07 OV-0. 147 10| Tin u2 >V-o. ^ 2 iamihaii hJ 

R 0.10B5V947n 09 OEC 0.644452940 00 R A-0. 1 14492810 01 V 11.355 |0354l| 02 PTH 0,244347440-01 ai ii.934419i< I ii 02 

MUMbER OE STbPS IS 0. lOOOOvUUIl 01 


MARCH I4t 1971, 15HRS, 54MIN, iH, 22SFC J|>LIAU OATt 2*4li)25.i'T*»17*>2#. 

VENUS CEnTEREO CEMTRBL H1R1V MEAN EARTH RnUATlB AMT) EnolMtlV Of iqSo 

4-0*644 T66U 31) 04 V-0.9 11 1 HtS20 05 7-0.153815440 05 OS 0.40740092U 00 PV 0.4407*17550-0 1 02-0. 23«i1 Mol >’ IN) 

B U. 10Ul)15T4r) 06 OEC -0.2071 74440 f )2 fl A-0.94074Q410 02 V ().451I149]U OO PfH I). 1 42 104 550- H a 2 0 , 112/>5>>ii') 01 

SNA 0.535O0359O 05 ECC 0,H6944052O OO 1MC 0 . 1 o 02 )n| 8 (l P 2 LAN 0,450421410 02 APE 0.449954910 U 2 RCA 0 .h 9 K 49 TK 8 n 14 

C3-U.AOTU4172n 01 THET 0,lft00nut>00 01 nPfir. 0.1U0O|5740 04 Mam U. 18000000,0 01 perm 0,178987420 0 ? TPF9 l). 7 B 9 *»«V| 2 '* iiii 

VENUS CEMTEHtn TRUE ECLIPTIC AMI) EUUIM04 OF nATb 

H-U ,4 1 2 744 SOI) 04 V-0,9q71R9670 05 7 0.44417715)) 04 l)S U .407440.950 IIO 1)7-0.480394440-01 0j-U.i£2HH9ii>i)>li iiO 

R U.l(R)0lS74n 04 DEC 0*247154710 01 84-0*935142410 O? V 0.45I114V1U (Ml PTh-0,847427I)40-uh a 2 il.UU40*7«i' ()1 

SMA 0.535003591) 05 ECC 0.R4944U520 00 n»r. 0,207497500 02 LAM ti,9i550«71u 0 ? APF 0.152447180 01 MCA ll.6V)*49T8«l" i '4 

03-0,407041 770 01 ThFT 0,lHOOflUUUO 01 APflfi 0.1000157411 04 MAM I'. IROOOOOOU 1)1 PERM 0,178987420 i)2 TPFh 0 , 7B955T1 2I< On 

SUN CEHTEKEO TRUE ECLIPTIC AMO EUUlUnU OR OATt 

K-0.40442T2U) OR V-0. lOOTFBOSO 09 7 0 * 99 Ta 55 llO 04 OK U,32H184R7U 0? DV- 0 , 1 11171 15n ti7 l)2-().27V19«i«'i i'| 

R V.10R43UU7D 04 DEC 0.499711520 OO HA- 0 . 111481350 ui V Ij,355I|R7B5U 02 PTH 0.5R715|440-U1 a2 0.93449174H ))2 

NUM 9 ER OF STEPS IS O.SRSOOUUUO 01 


MARCH 14, 1971, 3HRS, 4Wl|N. ia.373SEC JULIAN UAlb 7441024. 45NT7T4T 

VENUS CENTERED CEUTRAL OOOV NbAN EARTH ECHIArnR AbO EOUtWIK (IF 1950 

8-0,45937(1710 04 ¥- 0,911205540 05 7 - 0.154079490 05 OX U. 4044921 IU 00 DV 0 * 445 OA 9 |)On-Ul 0 /- 0 * 7 iW 2 N 1 8 1 r> Oli 

H 0 * 104 ) 029451 ) 04 l)FC-n* 2 nT 504 ISO 02 RA- 0 . 9404 l 451 () 02 V Ii. 45 fl 490510 DO PTH-O, I 22589940^)9 47 < 1 . 1 1 228 la 7 '> li-) 

SNA U.SSSOOftRAD 05 ECC D*I694RR7S0 00 IbC, n. 1 u(l 422 tin U 7 LAN U.450R1h05U 02 APE U, 4 A« 9514 R 0 U7 RCa D.bVTlTA^I'' |I4 

CS-0.607USHU20 01 TUFT 0. 1 ROODUlHfO 01 pPMO 0 . 1 U 0029 M 0 04 MAM U.lAODOUOnU 01 PFRM 0,17MU9D)*90 U2 TMFk ( 1 . 7H95441>*)' iMi 

VpIPiS CEhTEbEU true ecliptic AMD t'AltMOB of OATt 

K-D. 40734 B50II) 04 V-D, 997171590 05 7 0,4b41('2N40 04 OB 0.40494VR4U ftO IIV-0.4TM 47410-01 07-0, JJM977 D*u* |i i 

R 0.100029630 06 DEC 0.265926040 01 6 4-0.934844090 U2 V U.65045U51U OO PTH-0.0792T425O-0A A7 U.11U41A47)) i'1 

SNA U. 5354)06040 05 ECC 0.669696730 0(1 |WT. 0.2079794)40 U2 LAN U.91519300I) O? 8 PF 0.15249)870 u1 RCa Ii. 497) T422i> dm 

CS-U.6U7|)S602() ni THET 0. IRMRhtUUI) 01 aPOO 0.104)079650 1)6 NAN U«lROnnU4tOO 4)1 PERM 0,1TP«44M9h u2 TPFH ll. 769564150 (h) 

stm CEHTCREO TRUE ECLIPTIC AN|) feCtulNni) op oaTe 

S-^.16226SS60 4)0 V-O*10244|p)6O 09 7 |),667809T8n 06 Ob U.351R2H6U 02 ITV-4i.liai«94?n 07 <i/-ii. 22 Hl 45 | 5 r> u] 

R U«14)R659660 09 DEC 0*957t349uD 00 R A-f). IU94T5290 01 « U.955D6384t) 02 PTM 0.90TI6A540-UI A7 l)*N36k7927o u 2 

MRiOER OE STEPS IS 4).S«104)(N|UA 01 


HAKN U, 1971, 2SHRS, 59H1n* 59.9995PC JULIAN iaTF 2F41tii( (.M manRhi* 

VENUS CEH 7 EHKD CfMTRM. ROOV MfeM EA8fH FOl)iTrs« ano AQUlNim HE 19541 

R-0. 196943710 4»5 V-0. l2436Sl5n 05 7 4),7n7745*»[M> 04 D»-4i. 194477 <90 4.1 UV-0,4l2141l5o t)l "i(-t),»«677Fl5li il' 

R 0*2949MOM) 05 OEC 0«74t549«7wn ni OA-O* 147712064 1)1 V I..AA4A5415M OD P7*4 li.41A1|]2«n U2 Aj |i. tl92nitu70 ii 1 

&4R 0.5S494SB06 M ECC 0,669*42220 00 INC r),5u04991tn 02 LAN 0,441U1S7*b) 4)2 APE 0.A4997O140 ^i2 RLa U,69hlA24*U' U4 

C 3-0.6072216141 01 THE? 0. 120*47250 01 fePOO U, lil4)RD7490 u6 mjo, 0.11TA74«7lt 02 PERU 0.17H9t9)i2H u2 TPE<4 ll.•,147^1M■lH.ll^ 

VEWTS CINTERiHl T*UE ECLIPTIC **!» EOUIMAR UE n«T6 

H-0.B95496540 4)» 7-0,1054533*1) fW £ 4V,T*41«52»)) It* l*R-4i,i*21l642U «S l<7«4i,*14A(l74Tn ul U( <i.A2M*(4«4r 1,1 

R 0,a34«6<H|gn 03 OK 0.1*0364710 02 RA- 0 . 151661544 Ul V Uo444455|1U Of NTH |1„5*>A1tg 29u ,•> pj ••.UA5A97AAn 

S*6 0,534*65110 05 PCC •o**9*4222* 0* 1*C 0«2U*lH»*>*n U7 LaN ti,9154f»772« l>2 Ape **aR“75I))9AD 01 >*4.A u,A94l42«(*r ,i« 

C3-0,607Sil6Oa 01 THET 0,120*47254 *1 8P4R *, 1U4H*I|7(UR) U6 NA4 M,|17424RAU U2 PERM 0,1788)9010 ll? |84 m II, s 1577184*0-4) I 

SUM C«4T««64 TMlt PCLBPTIC *W k<JUt*Mt)l V DATE 

H-«. 3364* 1140 m 7-#» 10SB936PO 69 2 4,92ll4l49(» *6 »» 6,31S15991B 07 »7-*i, 1*241 14570 U7 Di-UaLd t'l 

* O.B**596*gO 09 OfiC »a2WllRl639P) (9) *A-OaTl»*B6t49n 05 V Ub34«‘*5592u 02 PTm 4i.7«a 99AN20 4>| »/ li, v24«44N2ID 07 

« H |M|>| t6 OO 6Tf»5 B5 Oa06l«*«IMA *9 





CASE 2 Periapsis Printout* 

The only difference here from ease 1 is the input array c(776), which now re 
quests output at periapsis. Notice that the initial position was at periapsis. 


^Pcfiapiis pfiMoMt is a fans of osoct past sod le^dtos teststdog ihm iotcginnofo 
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J«WiT 


73(P 

u.b-jsoooooo 

OS 

TSl 

ll.kA'I ISkHU'I 

00 

TSZ 

1 |0 


TSS 

ii.9S0ri(ifiiKin 

it/ 

739 

0 . PIUIHIOOOOO 

(1? 

ISS 

fl, AS04Hi('(Wlfl 

(V 

TS6 

'',Tlo3M(>0(l 

09 

7ST 

0.0 


TSR 

U.Tl03lT(Ui:i 

(19 

760 

ll.O 


761 

0.9(l1i|)|>ll00(l 

Ul 

763 

0 . 1 ii(i,ii-iiti(iii 

06 

7h9 

0.0 


766 

0. 4ti(i(iuOii(,n 

01 

7TS 

('. l(1(>()(ltM)(|(l 

Ul 

776 

-11. 1(10110(10(1(1 

01 

777 

0 . 100 101 OOO 

Oft 

0 

0.0 


0 

0,0 


0 

(1,11 



is >«IJT KFOIIESTf:*) 

HAH<4rM|CS IIF FAKTH 

MAHHrMics n»: ^fin.» iMCLimfiti 

kA^MiMTCS tte (9ANS IS iwr.Ldtifi) 

MH.&H 4AI11AT1IIN PbIfSSlIME IS ^(IT InCUIH^O 
|nPI)T in HlrAh pa APMt PO ItF 144(1 

«*«Ch n, 1471 1 OHAS, «n|M, u. USF.* JI<t|AN (lATt aAAllWi.MinOOOfm 

X 0.4Aft41UH9» US V (1.«44|nHtlHn n* r n.?A7A«TlBa 0* Ol i1¥~n,A27nl 404(1 UH oz 0.3il4iA4l 4(1 «1 

« (>A OFC II. iiOTllAltl^n n? UA 0,A»ifli»S4*ITn o? w U.OXIPAUIOI) in PTM ».A(l49ATd40-«A U 0. (>779249411 


nAMCh n, 1971t IIHHSr UhIMi 0. OSFC .IHtlAM iMTfc 299lU2ii.41l(H’IHMHl 

ypNIIS CFk 7 PhP 0 Cf 4 T 9 AL HIIOV OfeApI EAHTH FOIIoTOtt F(M‘ iNliB IW= 14 S(* 

» U, 9 AH 910 *t 2 |l P 3 V O,A'i'<ll(MiH 00 09 / 11,29 79 h 74 K(I Ii 9 0 X-ll,HAnP 4 K <>20 0 | |IV-(I,A 74 h 1 4 * 14(1 UO r>Z 0 .J 2929414 H lU 

H U,A<I 999999 II 09 OfC 0 . 207091 ( 1^11 02 KA (I, RSn<> 44 qTll >12 V 0 . 441 24 U 10 'i 01 27 ^ ( 1.0040470 40 -U 9 AZ 0 . 97792444(1 (»2 

S'lA U.SlSOiHiOIMl ns FCC O.AAOlSHHbll OO |Mr. 0 . 4 u(|(>]i(H>|n u> |.AN |j. 49 llOOoOl|] 02 A 2 F 0 . 44 ( 110001(1 U 2 hCA 0 .(i 94 <« 4999 II 09 

C 3 - 0 .A 07097 } 9 iI O] T 0 F 7 0 , 4 > ppor, 0 . IwOOIHUHiO IM |., 40 S I lH 90 o -|9 PFKH 0 . 47 A 9 B 4 Aii(l 02 IPFk O.U 

VFMiS CFNTEItkII 7«tOF FCLIPTIC AMO feOllIxOH OF liATk 

X Ui 447 S 03 M 4 ij» OS V (l. 997 M 9 Ss 7 n (>9 7 - 0, 427714940 04 nn- 0 .n 9424 ]l 4 flU (II OV 0 , 992 |P 799 ri (JO O/ (l,i 29 H 7994 n ll| 

K 0.999999990 09 OFC-U. 29044^)40 Ol UA O.P(> 94479 <»n 02 V Ii. 94 | 29 ul 90 01 P 7 H ( 1 . 79724472 ( 1-09 AZ (i. 0992 Ai>H<in ii 2 

S 9 A U«S 4 SUOOa 9 tl ns fCC O.H 99 | 9 HH 90 (HI |MC 0 . 2 U 749794 O 02 1 AM 0.944949940 02 APF 0 . 44290409(1 04 9 Ca 0.699999900 09 

r. 4 - 0 . 907 O 9 *i 79 O m 79 FT 0.0 APOi; 0 . lUOOOOOlO 09 kAW O. 390 O 0 UO 0 U (<4 PFHO 0 . 47 P 9 M 4940 U 2 TPFh O.U 

SON CF^lTEHlin TPIM: FCLIPTIC AMD tOtllAinK OF OATt 

X-0.9S0I3M40 OP V-0.9n714979n OP 7 0.122029490 07 OX 0.27«4|90Tli 02 OV-0. 1 49A944P(i U2 0 / 0.12910249(1 01 

R O. IOR9996IIII 09 (IFr. 0.649902990 00 P A-0. 1 194 1 4(1(91 04 V (1.299909460 02 PT9 0,999999040 Ul AZ 0.979144490 02 

JIIHUFH OF S7KPS IS (I. mOOOUOOO 01 


hARCH 19. 1971, 1)HM$, SSMIM, lb,49MSFr. JliUAM 0A7k 2991U29.0I0HI4944 

VFMIIS CFM7FHFO CFWTHAL IOP07 OtAn EARTH FOIIATOM APO F^IIJMIIX OF 19S0 

n 0.969299100 04 V 0 . 6 S 1645 T 60 09 7 n. 29709149'1 iJ 9 OX -((. 669749010 Ol nv-ll, 9411 ?A 490 UO HZ 0.329674440 (It 

P. II. 99 H 971 S 90 09 OFC 0.207179740 02 “A 0 . 01 , 9799(910 02 V 0.942429440 Ol PTH- 0 , 794944440-07 AZ 0 . 6 T 7 T 2 A 920 ii 2 

S«tA O. ‘> 49996691 ) 09 FCt 0 .fl 99993 b 60 (Hi IMT, 0 . 4 oO 2204 rHl U 7 LAD 0 . 94092209 U 02 APF 0,999942940 U 2 »CA ■■. 69 H 4 T 1490 o 9 

C 4 -OtAOTUSUSR» «1 THft- 0 . 11299364 ( 1-114 APOp 0.1010119940 OA PAM 0.460000000 (i 4 PFttH 0 . 4769797(10 (|7 TPFH-O. 1 TU 92 A 690-09 

VFHKS CFMlEAkO TRUE FCLIPTIC AMO 6 (t«)*inH nP OATt 

X 0 . 92796927(1 09 V 0 . 69649 H 62 O 09 7 - 0 . 42449791(1 HI nH-U.H 70 t 064 Au 01 (IV 0.647990190 W 1 111 O.jUTTSOOTO Ol 

4 0,696971390 09 OFC- 0. 267 142690 01 6 a 0 .P 694114410 U? V II.V 3732944 U 01 PTM- 0 . 4 A 77 14940-07 AI 0.6939341690 02 

S «6 0 . 44499309(1 oS FCC 0 ,P 6 R 9 *S 37 o (Ml INC ( 1 . 2 O 799 O 290 07 LAd* 11.944403240 07 APF 11 . 442*97310 (J 4 PCA 0.696971490 09 

C 3 - 0 . 60 70409 VI Ol TMFT-O. 162443660-04 upor. 0.100014960 06 «AN 0.360000000 04 PppH 0 . 47 M 9 T 9290 U 2 TPEM-O. 2960 O 24 TO- 41 M 

SII 6 CEMTEhtrr 7 HIIF EClIPTIC AMO EmUMOK »F OATb 

8-0.9063931611 no V-n. 100622110 06 2 0.09244049(1 06 OH 0.234299661) 02 |TV-0. 124R04460 U2 Ol 0.123990010 111 

A 0.106426661) 00 OFC (1.997303990 HO RA-O.l IZOOOITO ()4 » 0.26TIO4TRU (1^ PTM II. 614776900 1)1 61 0.0730990*0 02 

P;l*HFit Of STEPS tS 0 . 216000(100 04 


mAHCH lA, 1971, 3MRS, 30HM, 33. 39SFC JHLIAm UA7E 299l(l26.044l|0997 

DEMOS CEMTEAEO CEMTHAi OUOV Mb60 cPIfH FOIIATiX 660 Fnt)|Ht|f (if 193 ;’ 

a O.4396902W1 03 V 0.630367230 O* 2 0 , 296 -, 69760 09 08-U .P70402 *20 0| l)V-n.64PVI17T0 00 01 0.33092744(1 Ol 

« 0.697124*6I> (19 OFC (1. 207306360 (‘2 PA o,A5**479020 07 V (1,943292270 01 PTH (1.44427137(1-19 61 0,6773961911 «2 

S«A 0.43499219H 03 FCC O.P69696STn OO IMC 0.40IM.2T490 02 1*9 b.94liPlli99u 02 OFF 0.999644270 02 6 CA 41.697139660.09 

C 3-0,9(17099 7 HI Ol 1^T-(I,9T9'/6733I1-('A npIWi 0. l00P2?|*n (>6 990 (1,360000000 04 PFMM 0.47*974240 02 7PFR-0. 1U1A9990(W0>I 

VFjMIS CFMTFPHi TPhF ECLIPTIC *90 t‘HII<«lK rs OATfe 

H (1.923376(1211 03 V 0,643(111396(1 (9 r-«l.424d,3602(' 04 OH-0,*T«P72*2i) Ol IIV 0,6*4327660 uO HI U.92P330A40 111 

(4 0.697129060 09 oPT. -0. 2639232VO Ol 96 0.663199270 u2 V U.933792?2u Ml P1(H-0.P6f(H|ln4n-0M Al 0*693699910 02 

SriA 0.939992100 03 FCC 0.069699300 00 |MC 0.2U79P24T0 02 lih« 0,‘3444BB644) 02 APP 0,4429«2n7n U4 AC* 0.697129*60 09 

CB-O,607u3466tl 01 tH(l7-0, l44l327l9n-09 Apor, O, lu002TL*( 06 060 U.36000V0OU 04 PEMH 0.47*974300 02 T*Eit-0.2O67 199*0-09 

SON CEHTCHKO 7*UP FCll*7IC AN«I Eoulooa n* OAV£ 

8-0.462192010 0* V- 0 . 1024436341 09 7 0.6626(19790 06 08 t>,24066*A*0 02 MV-O. II 197492*1 U2 Ol 0.123107*10 01 

0 0.1«'.33ni90 m OEC 0.399T2S6OO (Hi H A-o, m»96n72(l 04 V- 0.263317170 02 PTm 0,94799*370 01 Al b.«7364774n l >2 

MOpiF* « 1 * S7PPS IS 0.2I200UUOO 03 


•lAftCM I6« 1971, 24MKS, 396B4*. 39«999S*r. JilUAM MTE 299tUIt7,34MW10n«HI 

DE(M»S CEMlCNtO CE4*7*AL 6007 NfcA* E6*7H POiUTOR 6*01 F9MBM08 d* 1930 

H-0,5999«696n 06 7-0.92*793600 09 E«Aa469l99740 04 08 1i.6136«97*U 00 1(7 0,1 7**796341 (*I 112-0.1140109114) 041 

* O.99fl0O29AA 03 9M-0.214SM739n 4i2 96-«l.«24332690 O? V 0.6A4fh6|3«j 4*1 PTH-0. 14)3341 1300 07 *1 0.111607621 4)3 

S4U t».33313atf90 113 KC 0.«697734»?n 04* IMC fl*40M33|9n »7 B .*0 0.930*994*tt 07 AeW «.6*9«A124A U2 9C* D.6964I93NAA 06 

C9-4|.««6«7*6Bn 4*1 TM*T-0.l7*?(9667n 04 *|H)R f).lwW6124>0 06 hAM o.1«1449A‘0a 04 pp*«4 n««a»t6916n 442 TPFtt-0, 760122*40 (to 

9PMIS CC4>7HI6» 9M4P FCLgP39C *4441 6«l)tlM4R HP noTE 

*-«,492973«M)W 06 7-6.996732310 4)3 2 4*. 464H99*4*> (I* t»8 U.6t4p26l70 OO 147 n.61«91994l47'«B 01-9.2366124(241 OO 

« 0.«9mie2964l 03 HOC O.SO*3«S9M 08 *6 -«>d 9B970742O 02 9 0.664036860 041 07M-41, 1033*834)41 02 AS 0.1EW7B4B24i 0» 

(4.333861090 4)9 *«X 9.M9773676 OO 3«<C 1.20*8324)30 02 L*9 l,.w4923B339 02 «p* o.4323B«6*4« 04 3C« 0«696693*T4> 06 

C4MI,4A6*7«6in 4)1 70*7-0.87*396690 (*4 APno )).B«s4KI6 12416 1)6 664* U.16B«49*70 n3 **h)4 M.470866166 02 7P«40-0. 760122*441 09 

84*^ CffOTIMM) TOUT HUP78C «3» 6404 BKW 0* 0676 

l-«i«33l(3S9064t 00 9-0.194272?^ 1« 2 4).3a74(UI4Mf» 06 O* O.VM323a«IO *S 149-0. M363*80» 02 OS-«o2293?123*) 08 

* 0.84)60 7 6606 O* )I*C 0.272746331) 491 M*-*. lo38467in 04 (4 0.3366607*0 02 «t9-4)B84>62*8261 04) 62 0.937076681 02 

ni 6«*M W S762S !S 0»7l|l9*%M494a 62 
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MK le^ac exact pdAioiit tmd pdatowt <»r |»€fiipsis m saM« fan. TiUP 
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SMA U. 536001120 05 FCC U,PA9AHh 240 391 |9C n.3u04249Rn 317 I AN 3'.46|i4|266U 07 ApF 0.443*9 3 / 1 All u? K.A 0.69T1693Rit 04 

C3*0.63ITU4453I) 0| TNFT-0. 14R7R662n 03 pPOR 0. |l|3>n7H63n 06 MAM U,323HK444m iI 3 PFHH 3I.37"9A4790 U7 TPFp-o. 15n4l724h no 

VENUS CEMTErEO TRUF FCLIPTIC AN3J bOUINnX OF DATE 

X 0,2190611411 05 V-0.444756450 06 7-3'. T293T I "60 04 OK H.7093056TU-0I UV 0.767316760 ul l'Z-0.n3>.3yi>3 fyk-o] 

R (1.901021500 3)5 nFC-0.637IKV909rt 0| 9 4-0.636693360 317 V 0.26266139U 311 PThhi. 631366671)0 1*7 a2 O. IUNI07h1o >'* 

SNA 0.635001120 05 ECC d.R69AH»24n IH> |NT. 0.2u7979A|0 07 LAN li.93637740u 07 APF 0 . 36749 il«i 50 u3 hCA 0 . 6 Vrl 693 no 114 

C 3 -(1.607044920 0| TNET-0. 14"2R66zn 03 APOi; O, 1003124630 06 MAN ll.323R"444)J 03 PFhH 31.37A9H4T90 07 TPFn-O. |5 h417740 00 

SIM CENTEMEO true FCLIPTIC AND tuUiNnh OF RATE 

1^0.366429750 OR V-0. 107224 lOO 3*9 } O.AqXTSTNND 316 IIH U.37T94T«6|) 132 MV-3l,936777l>Nn ut IIZ-U.214l«i9| 3)> 3>| 

•I 0. 100595310 09 OFC 0,360756090 on ■6-0,109720640 313 V 11.34149570.1 3>2 Pfu-ll, 3A QP I ll«|n Ul AZ l),935T7RAHn 11 / 

4I1MWH OF STEPS IS 0.9A60000UO 03 


MARCH 16. 1971. 23HRS. 54M|H, 69.999SFC J'R.IAN imTE 7441u< 3 y.6fPNP»P>0 

vends CEMTEREO CENTRAL WIOT HEAM EARTH FO|IATOH AMO FflillMlV OF 1950 

»-0, 3949949010 04 V-0.920T4M2n 3*6 7-0. 6«31 64740 (16 OK U.*134«474u on w 0.17AR7696N un nz-«*. lnUlll)44)i 3P> 

A U,96M(mi375p 05 nFC-0o?|43RT420 02 . RA-0.924353160 U2 V 0.6633i65f*2U OO PfH-0.eN667RMH* U2 aZ 0.11|6M264n |i« 

SM U. 535151420 05 FCC 0.A69775940 00 fRC 0.63|fl65614n 02 IAm U.A6(MH»96Ri 1 02 AIlF 0.449N4IP30 02 "CA 1). b96MU6A4i. ||4 

C3-«,AQAtl76U40 ni THFT-e. 1TP39E730 03 APtm 0.1003*61320 06 RAN 0.191339363^ 03 PFRH li.37«|A4A|n up T MFct-i*. I40l2463n (>» 

VEMIIS CEMTEREO TRUE ECLIPTIC AN|» ENUloNlB '*F MATE 

X-U.3429R3740 04 v-fi,99«i7534)Tn 06 7 0.549102670 334 OR 0«6|32233V6w Oft |tv Q,4|«R7«T3n-3i) i>Z-3i.294MI2693) 3R* 

fl 0.99Rr'(l37Rni 1*5 OFC 0.200504400 Ol RA-0.9197aT760 07 V 0.463065N2U 03* PT04-O.)O667«6rn 07 AZ 0 . 1 l 0 Tl 6 | 7 i> |i3 

SNA U.5391 51420 05 FCC 0.H697TS94O (*A BMt 0.20RI 323*61) 02 Lam 0.9352%1A6U 02 Am 3P 0.3626I9A90 u3 «CA 0.69M>N6AiN> t)4 

C5-U.404R7404I) 01 7*3FT-3). 170914730 3*3 Apf«A 0,10003)1320 06 NRN 0,141359360 03 PFMN 31,374)44610 02 1 PF«H». 740)24630 (it) 

SIM CeN7EREO T41K FCUPTIC AMO EOUlMriii Of OiA*- • 

R-«. 9943)9940 M V-0, 105272740 1*9 / 3*. 6173011673) 3*6 OK 0.33A623R3Nt o2 |.9-n, |««664220 o7 ♦H-0,22N5722«** »l 

R 0.104474440 M OPC 0*2727M55fl P3) BA-3>,10«)547|n 03 V U.3644POV6w 3*2 PTuhI, 104274469 3M) AZ ».9!l7liV443«0 |i 2 

HIMOfR nf STEPS as l*.5563K6HrtK* (.3 
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CASE 5 First and Last Point Printout 

This case represents the most efficient type of run but also the one that gives 
the least information^ being useful only when the final point is all that is de- 
sired* Note that the total steps ts^en will be the same as for case 4. The 
differences in input from case 4 are: 

(1) c(763) Is again set to 100 000* 

(2) c(764) is set to 0. 
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I 4 il II l( Hi !■ I • * uf\ 
I > 0 14* I Mill' III' 


^HAUrltl tS Mill Khi.')IF-STH’ 

HAH"<II|U|I'.S IJf “iiHTtl rS I'*'CL1IIJt-ll 

H6M*-lln tr.S (Ih l*llll^ IS |Mrt.NllHl 

naf4llniics UF HAHS IS 

bllLAk *<aIj1a1 IIII*' Fki-SSHI>i- Is ki|iT InH.LtliH-ii 
IMFIJT IN Hhatii l-ll Anil) Fy i|4- lUSll 

’'AKf.H IJ, i*Y7l, rwRS, (iHlN, 0, uSfi. ,1IILI/If‘‘ miTi: i,Siii'i'''i|iii 

X u.4>bH<9|(jHj;n OA V n.hsx inHHnii n& i yAt-vx!*ii ii« ii*-ii,hahh*>a«;'i> (O i v-(i,is;'1hi snsn iji) ’ i>t II, ic'ii'A xi I'l 

H ll,ft<f*jg*)'Y*J‘YI) 04 riFC. (>« /OY04H1 ^lt H? hA n, A4, a!>HX 07 h ii? ii, 4^ I ^4u| n ij ii] hT*' n, *4 i'Xo t7 h hi»-iiis ti, sV /v;'X 4 F|i 


'•AKCH 13» 1*171* (IHMS, (IH|M, 0, l(SF(; JI'LIa*' ijJiTt ,i4sli)iti,t>|i|'l>0llli(i 

VfNHS CE«T6XtM central KOMY ►-twii tAHTO FhhaThh XnI> Fill|nn)r IIK |4S*> 


X 03 V r;.h*i3lriHi*Mn ii4 / c. ?<,74 «7'<«n ,i^ 

H 04 OFC (i. jioTOAH i ^0 i.f -te. ii,f'*jj.N3iij70 OY 

bi*iA O.SS^OOOUOO 04 FCr 0,MA«»ls«()n0 00 | >>'r o, Xi,OOOo(* 1 1’ 07 

C 3"U, 007044790 01 TnFT 0,0 aHOO 0, 1*j|iIiOOOO'* OA 

WFnUS CFMFMtll ruilE FClIFTlC ANIi tOolAliX OF OfTe 

X 0,43? 4U3HOI1 03 ¥ ti,A*J7>i*Y3*»yi' o* 7 -o, 77 1 X AXn nx 

H 0,6*1 *4999991) 04 nFr.-0,?6(*?31>340 Hi wa I), f o4»*m 74 iji i „/ 

SMA 0, SBSOOlJUAO 04 FCf. 0,069191(1**40 00 1*0. (*, ?u7X97o*ii (iV 
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SON CFNTFOEO TBOF FCLIPTIC AN|i HUOlMiX OF 0^77. 
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V 0,h31?4u19o <lt PTO il.M S*|FTKX'l-oF a/ 0,((7 f¥;»XaFli ,1,1 

l_ar* 0 ,4* OI lOo** 1 1 J "7 a¥F |l, 4 SOl lOOO 1 O li7 vCa O, |jV*J*7¥¥tl¥l' ijA 

•■0* li,4*JSll|'aiH)-l6 PFAO 0, 37 m9U 3F|I|' (O' I F>K 0,0 


im -II. nft*./ X (4aij o| HV ii.fiu?! HjiaAi* o.i o, 3,!Fk7*if30 hi 

¥ o .¥? 1 , 1401 '>0 "1 Hi** l■.7'JV?'■37?l'-^J6 417 o.b<74^6*JHUii i,;, 

tfli" l•.u3SFX*l'•4l* 0? hkF O. 3474(’Xj;Hn ox 7(.« O , (,*7'(**'490'H1 n'l 

Frt'i ('..(MHIOHoOii MX |JH(0 II, X7W9H XFSI* O'* I AFK 0,0 


HX I . .'>*4X I 4o7ii 07 iiY-ii, I Xgmi*iXim o7 0/ o, i/„l(iyxj, . , | 
u ii,,;6Fim4 F"ii 07 P|t( m, 6 mo',**fm xn yj *7 i,,,. V4l X3*.*,i, o/ 


MARCH 16, 1971* Z3HSS* 49H|9* 99 , *(>J*/5Fr'. JOLIi* '- M" 1 t 7441,'? 4. SOI- 

VFNOS CFNtFkEO CFNTRAL WIO¥ AFaN FAKTh FiiimTC'* n *'0 F''ii|.iinr ,|f 1 .imi 


X-<),394994Hnn 04 V'0,9 7464*10 MS 7-{>, ?ftX I >,4 7411 oS 

K 0.99R003T40 04 0Fr,-(t,ilX3K74?O O? a A*-(i,‘t/4 3 S3 1 F.i (iv 

SNA (I.S341414?0 04 FCC 0, H6977494i 1 OM 1 nF. 0,XuMS5sl4ii i ? 

r 3-0. 606 1*740411 Ol rHFT-0. 174X96730 (X aP'iS M.U)006|X?>1 OS 

VFNtlS CFNTFkFO TIKE FCLIPTIC AOIt Ff’U|**nx IIF OATf 

X-0. 1429837611 04 Y-(l.996 V4>07rt OS 7 0 ,X|jX lOPS 71 , o4 

H 0,9980113760 04 OFF 0, 2oH4046li0 M I *< A-o , o 19 707 76M o7 

SNA 0,634141420 04 FtC (t,h6(»774940 MO l*r. 0, ?i'A inPOH,) 07 

C3”O,606874<i4O 01 THFT-(I, 17H3967 30 Ii3 aki'ft o, |*,tiMM '*?.i os 

SON CFNTEkfcli THOF FCUPtlf ANll t'lHiNOX op osTE 

X-0. 338319940 08 Y-0. lti3PT27hO (**> 7 1 '. * I 7X0 86 7n 06 

M (1.IOH6T4460 09 Off. 0.277T3844O MO ■( A-0, U'aI , 3 « 7 lo i,x 

»I) 8 hFm IIF STEFS IS 0,53S0OOUi'O 03 


OX l.,6lX69676u (M* OV I.I,174W7 Sm60 um I ?-o.16i'11i fil' Ml 

¥ ii.h6Xn4M'?ij (til MT**-o,|M«S7*i6rni u7 ni < 1 , 1 1 1 6 «,' 6 a-, ' ■ ( 

1 A *1 1 I, 4 SO*>*IVX |'0 m7 6 HF II.4A6.W. I 7 x 11 li 7 ^(,rt M, 6 ** 6 fl¥ S ‘I *J ' - i-i 

AA'I M.HlIXXyXSij MX MFkm ii,37ai66S10 (i7 1 PF><*,tl. I'O'l 76 * Xi 1 


oX i:,6l X72366IJ MO OY 0,Q 1MH79Y3I i-O I O/-I 1 ,? a 4«< I 7SXl Mil 

u (i,66Xos^Ii7ii Mi jM4t,V"min 1,7 fti (I, 11,1 4 1 X I 7-1 ,m 

KA'8 I'.'sXS7S 1 ASiI ,'P *ipF *1, *4741 96 ,0 (jX ¥(;« M,696F9S''|*I 1 114 

I’ll'' O, [U1XX¥X4|J riX Jaah M,XV9166Sl'l 07 | ¥*■< -M, Y«*u I 7* **XI I till 

■*X M, -SX6S7X U,1|, 1,7 ,|V-M. 1,196447211 iM H?-,).f29.(727nM I | 

¥ M, *S44>4||74|* H7 MT**->I, 111,776860 MO A/ M , 9 4M 466 1 II 1 M7 
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itHAimu IS kF(mi*^STfiv 
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1 . 1 ifi.. 

ts 

iHCm!iii-o 

Far -0-1 ICS 

of 

••A*S 

IS 

IMCLIiOFO 


SKLAH «flrn*TlfW' P«FSS(I»F tS l^r.LI'^*•=l: 


IWPIIT |p4 mi-oF F» 
^ahCH in, 1>*TI, 
X-U. 60267711^0 
ft O,66O*>S106 |i 


AhO FO OF l<»6ri 
17MRS, 5M*«, 9**,V‘)9^«C 
06 V n.6 7?077itn im 
06 iiEr.-0.2^fl67i;*wn n? 


7-n,?67?2't060 Ii6 
Rfl 0.176621160 Oft 


JItLIAP iMTt 266l02l.>l?6Oonn 

OX“L'. 206061 JftO 01 0V-0,10?ftfl26l0 02 02 0.27S37I6MO III 

o 0.10HO60RO0 02 PTH 0.7l612ftft6'»-U6 02 0.736021600 02 


MARCH in, 197), ITHMS, iPlM, >i9.909<PT. 


JI^IAF* 00 Tt 2661 U2 1.21260000 


►AkTH CFHTFRf" 
X-0.6023ft««6|i 
k 0.660951090 
Kria 0.106236700 
L3-0, 36260197(1 


TBOF FaHTH FOIiaTIIB AMO FOOlNOr OP OATH 
06 V 0,46ftl6rtK6f1 Oft 7-<*. 26"ATA 6»t 06 
1*6 nFr.“O.26966U010 02 0.176790660 lift 

06 FCC 0.936591330 OO I*i(; 0.2K7671 «ho 02 
01 THFT U.n npor. O.201H66MJO 06 


OH>U. 1997671 60 01 
V O.10Mfi69»no 02 
LAM 0.229656610 Oft 
PAM 0.360000000 Oft 


oy-0, i«76fni9o 02 
PTH o,706tift>)aon-«6 
APF n.ft026ftn62f» Oft 
PFMH 0.oftnftft2t1ft0 02 


02 0.270917720 Ol 
A2 0.735966020 02 
9CA I1.AM1951090 06 
TPFh U.U 


FAkTm CPMTFKrO 
ft-O. 602677030 
q 0,66(i95l09l> 
S~iA 0,106236690 
(,3'-0,ftM2A02020 


CFWTPAL MtiOy PtA6 6 6FT6 toO*TfH Aftri OF I960 

06 ¥ 04672077360 Oft »-0, 267?2ft06O O* 206067990 Ol l>V-(l,ln2ft"?510 ii2 

06 OFC-o, 2ftH67i2UO ('2 P* 0.176621160 Oft W 0.L06m9R90 02 PP' «,716123ft5l’-06 

P6 Fcr. <».9ft6S91320 on Ii-r 0.2A67T69in «2 LAM 0.2?9Aft6690 Oft APF f1.ft0269AHTn Oft 

Ol TOFT-n.lF6l 756^0-l'6 Apor, 0.2ulft6ft67n «6 9AM 0.36IHtOUr>nO Oft PFKH 0,9ft0ft31A6n 02 


02 0,27937l6fli< 01 
A2 0.735H21600 02 
RCA 0,660961090 06 
TPhH-0.229MflH7«o_09 


FilftTO FtHFO LATITOOF lOFr, 1-0. 239660010 02 UI^M TIlOF lOFRi 0.1U67l<530 Oft 


Mi-hl-H IIF STtPii Is O, lOIHKK.01"O m . --- - 

PHfluF FiJTFMFO FARfh PFA»l«>rK6A. Tl*9 AMO OROIT APF r.lWFM RFLOlf I 
-Hi0.h 10, 1971, I7HRS, 160H, 23. lft6ft6C I 

P^O.tF r .f PRiiO FAhTH T1l*F A"**j OJHIT Piyso JFI.Ow / 

. 10, 1971, I7M9S, |6»|9, 17.22259r 1 

> TMtoW€T 

wkiloF OAS Ltn FARTh PFMOmhra, TI>F aF'» ORpIT Akt RlMFH FFLOV' \ 

-arch |0, 1971, '.7HM5, 16F|M, 23.2I66FC \ 

pKiIkF has left FAMTH OiPPA. TImF a,JiI h»ht 6RF RIVFM hFLOM j 
►.ARCH |0, 1971, I7MRS, Ifc l**, 17.301 6FC 

HftlHF FmIFHMI PAHTH PFfcMO-lPA. TI-|F A ,9 IVMIT A9F r.|VBM MFLOM 
iaNCH 16, |97|, IftHRft, RFI-i, 6P,6-i95Fr. 
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pftOhF HAS LFFT FARTH PF»06 «I«a. Tlr.f aMO IIUkIT AHt P|9 Fb- HFLOP 
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